CN108169316A - 基于支持向量机和树叶磁学的大气重金属污染评估方法 - Google Patents
基于支持向量机和树叶磁学的大气重金属污染评估方法 Download PDFInfo
- Publication number
- CN108169316A CN108169316A CN201711380579.9A CN201711380579A CN108169316A CN 108169316 A CN108169316 A CN 108169316A CN 201711380579 A CN201711380579 A CN 201711380579A CN 108169316 A CN108169316 A CN 108169316A
- Authority
- CN
- China
- Prior art keywords
- sample
- leaf
- heavy metal
- atmospheric
- sampling
- 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
- 229910001385 heavy metal Inorganic materials 0.000 title claims abstract description 65
- 230000005291 magnetic effect Effects 0.000 title claims abstract description 62
- 238000000034 method Methods 0.000 title claims abstract description 40
- 238000012706 support-vector machine Methods 0.000 title claims abstract description 31
- 230000005389 magnetism Effects 0.000 claims abstract description 23
- 239000013618 particulate matter Substances 0.000 claims abstract description 15
- 238000012360 testing method Methods 0.000 claims abstract description 12
- 238000011156 evaluation Methods 0.000 claims abstract description 10
- 238000004458 analytical method Methods 0.000 claims abstract description 9
- 238000010276 construction Methods 0.000 claims abstract description 4
- 238000005070 sampling Methods 0.000 claims description 47
- 239000003344 environmental pollutant Substances 0.000 claims description 17
- 229910052751 metal Inorganic materials 0.000 claims description 17
- 231100000719 pollutant Toxicity 0.000 claims description 17
- 238000012549 training Methods 0.000 claims description 17
- 230000029087 digestion Effects 0.000 claims description 12
- 239000000243 solution Substances 0.000 claims description 12
- 238000009826 distribution Methods 0.000 claims description 11
- 229910052742 iron Inorganic materials 0.000 claims description 11
- 238000005520 cutting process Methods 0.000 claims description 10
- 229910052745 lead Inorganic materials 0.000 claims description 10
- 238000011160 research Methods 0.000 claims description 10
- 229910052793 cadmium Inorganic materials 0.000 claims description 9
- 238000012795 verification Methods 0.000 claims description 9
- 229910052785 arsenic Inorganic materials 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 8
- 239000012528 membrane Substances 0.000 claims description 8
- 239000000919 ceramic Substances 0.000 claims description 7
- 229910052804 chromium Inorganic materials 0.000 claims description 7
- 238000001035 drying Methods 0.000 claims description 6
- 229910052748 manganese Inorganic materials 0.000 claims description 6
- 239000011573 trace mineral Substances 0.000 claims description 6
- 235000013619 trace mineral Nutrition 0.000 claims description 6
- 229910052787 antimony Inorganic materials 0.000 claims description 5
- 229910052802 copper Inorganic materials 0.000 claims description 5
- 230000005415 magnetization Effects 0.000 claims description 5
- 229910052720 vanadium Inorganic materials 0.000 claims description 5
- 239000004698 Polyethylene Substances 0.000 claims description 4
- 238000000605 extraction Methods 0.000 claims description 4
- 229910052759 nickel Inorganic materials 0.000 claims description 4
- 238000005457 optimization Methods 0.000 claims description 4
- -1 polyethylene Polymers 0.000 claims description 4
- 229920000573 polyethylene Polymers 0.000 claims description 4
- 229910052719 titanium Inorganic materials 0.000 claims description 4
- 239000013598 vector Substances 0.000 claims description 4
- 229910052725 zinc Inorganic materials 0.000 claims description 4
- 229910052782 aluminium Inorganic materials 0.000 claims description 3
- 239000008277 atmospheric particulate matter Substances 0.000 claims description 3
- 238000013210 evaluation model Methods 0.000 claims description 3
- 229910017604 nitric acid Inorganic materials 0.000 claims description 3
- 239000003755 preservative agent Substances 0.000 claims description 3
- 230000002335 preservative effect Effects 0.000 claims description 3
- 239000010453 quartz Substances 0.000 claims description 3
- 229920006395 saturated elastomer Polymers 0.000 claims description 3
- 238000007789 sealing Methods 0.000 claims description 3
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 claims description 3
- 238000005303 weighing Methods 0.000 claims description 3
- 239000012490 blank solution Substances 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 claims description 2
- 239000000356 contaminant Substances 0.000 claims description 2
- 238000012937 correction Methods 0.000 claims description 2
- 230000005347 demagnetization Effects 0.000 claims description 2
- 230000009977 dual effect Effects 0.000 claims description 2
- 230000003203 everyday effect Effects 0.000 claims description 2
- 239000012634 fragment Substances 0.000 claims description 2
- 238000012417 linear regression Methods 0.000 claims description 2
- 238000013507 mapping Methods 0.000 claims description 2
- 101100545275 Mus musculus Znf106 gene Proteins 0.000 claims 2
- 238000004368 synchrotron infrared microspectroscopy Methods 0.000 claims 2
- 230000001360 synchronised effect Effects 0.000 claims 1
- 238000004088 simulation Methods 0.000 abstract description 17
- 230000007613 environmental effect Effects 0.000 abstract description 8
- 230000002265 prevention Effects 0.000 abstract description 2
- 239000005427 atmospheric aerosol Substances 0.000 abstract 1
- 238000011002 quantification Methods 0.000 abstract 1
- XEEYBQQBJWHFJM-UHFFFAOYSA-N Iron Chemical compound [Fe] XEEYBQQBJWHFJM-UHFFFAOYSA-N 0.000 description 17
- 239000011133 lead Substances 0.000 description 13
- 238000012544 monitoring process Methods 0.000 description 11
- 235000019083 Osmanthus fragrans Nutrition 0.000 description 10
- 244000242564 Osmanthus fragrans Species 0.000 description 10
- 230000006870 function Effects 0.000 description 10
- 239000011651 chromium Substances 0.000 description 9
- 239000011572 manganese Substances 0.000 description 8
- 239000002184 metal Substances 0.000 description 8
- 239000010949 copper Substances 0.000 description 7
- 239000011701 zinc Substances 0.000 description 7
- 239000002245 particle Substances 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 description 4
- 239000011707 mineral Substances 0.000 description 4
- 235000010755 mineral Nutrition 0.000 description 4
- 238000003908 quality control method Methods 0.000 description 4
- 230000000875 corresponding effect Effects 0.000 description 3
- 230000005294 ferromagnetic effect Effects 0.000 description 3
- 238000009616 inductively coupled plasma Methods 0.000 description 3
- 238000011084 recovery Methods 0.000 description 3
- 241000196324 Embryophyta Species 0.000 description 2
- 241000282414 Homo sapiens Species 0.000 description 2
- UQSXHKLRYXJYBZ-UHFFFAOYSA-N Iron oxide Chemical compound [Fe]=O UQSXHKLRYXJYBZ-UHFFFAOYSA-N 0.000 description 2
- 229910000831 Steel Inorganic materials 0.000 description 2
- 238000005299 abrasion Methods 0.000 description 2
- RQNWIZPPADIBDY-UHFFFAOYSA-N arsenic atom Chemical compound [As] RQNWIZPPADIBDY-UHFFFAOYSA-N 0.000 description 2
- BDOSMKKIYDKNTQ-UHFFFAOYSA-N cadmium atom Chemical compound [Cd] BDOSMKKIYDKNTQ-UHFFFAOYSA-N 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 201000010099 disease Diseases 0.000 description 2
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 2
- 238000001095 inductively coupled plasma mass spectrometry Methods 0.000 description 2
- 238000002354 inductively-coupled plasma atomic emission spectroscopy Methods 0.000 description 2
- 239000006249 magnetic particle Substances 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- PXHVJJICTQNCMI-UHFFFAOYSA-N nickel Substances [Ni] PXHVJJICTQNCMI-UHFFFAOYSA-N 0.000 description 2
- 241000894007 species Species 0.000 description 2
- 239000010959 steel Substances 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- OYPRJOBELJOOCE-UHFFFAOYSA-N Calcium Chemical compound [Ca] OYPRJOBELJOOCE-UHFFFAOYSA-N 0.000 description 1
- 208000024172 Cardiovascular disease Diseases 0.000 description 1
- RYGMFSIKBFXOCR-UHFFFAOYSA-N Copper Chemical compound [Cu] RYGMFSIKBFXOCR-UHFFFAOYSA-N 0.000 description 1
- DGAQECJNVWCQMB-PUAWFVPOSA-M Ilexoside XXIX Chemical compound C[C@@H]1CC[C@@]2(CC[C@@]3(C(=CC[C@H]4[C@]3(CC[C@@H]5[C@@]4(CC[C@@H](C5(C)C)OS(=O)(=O)[O-])C)C)[C@@H]2[C@]1(C)O)C)C(=O)O[C@H]6[C@@H]([C@H]([C@@H]([C@H](O6)CO)O)O)O.[Na+] DGAQECJNVWCQMB-PUAWFVPOSA-M 0.000 description 1
- FYYHWMGAXLPEAU-UHFFFAOYSA-N Magnesium Chemical compound [Mg] FYYHWMGAXLPEAU-UHFFFAOYSA-N 0.000 description 1
- 238000010220 Pearson correlation analysis Methods 0.000 description 1
- ZLMJMSJWJFRBEC-UHFFFAOYSA-N Potassium Chemical compound [K] ZLMJMSJWJFRBEC-UHFFFAOYSA-N 0.000 description 1
- HCHKCACWOHOZIP-UHFFFAOYSA-N Zinc Chemical compound [Zn] HCHKCACWOHOZIP-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- XAGFODPZIPBFFR-UHFFFAOYSA-N aluminium Chemical compound [Al] XAGFODPZIPBFFR-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011575 calcium Substances 0.000 description 1
- 229910052791 calcium Inorganic materials 0.000 description 1
- 239000000969 carrier Substances 0.000 description 1
- 238000009614 chemical analysis method Methods 0.000 description 1
- JOPOVCBBYLSVDA-UHFFFAOYSA-N chromium(6+) Chemical compound [Cr+6] JOPOVCBBYLSVDA-UHFFFAOYSA-N 0.000 description 1
- 238000002485 combustion reaction Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 238000009792 diffusion process Methods 0.000 description 1
- 238000010790 dilution Methods 0.000 description 1
- 239000012895 dilution Substances 0.000 description 1
- 239000000428 dust Substances 0.000 description 1
- 239000010419 fine particle Substances 0.000 description 1
- 239000000446 fuel Substances 0.000 description 1
- 239000012535 impurity Substances 0.000 description 1
- 238000009776 industrial production Methods 0.000 description 1
- 235000014413 iron hydroxide Nutrition 0.000 description 1
- NCNCGGDMXMBVIA-UHFFFAOYSA-L iron(ii) hydroxide Chemical compound [OH-].[OH-].[Fe+2] NCNCGGDMXMBVIA-UHFFFAOYSA-L 0.000 description 1
- 239000011777 magnesium Substances 0.000 description 1
- 229910052749 magnesium Inorganic materials 0.000 description 1
- 230000005381 magnetic domain Effects 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 1
- 229910052753 mercury Inorganic materials 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000011591 potassium Substances 0.000 description 1
- 229910052700 potassium Inorganic materials 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 239000012925 reference material Substances 0.000 description 1
- 239000013558 reference substance Substances 0.000 description 1
- 208000023504 respiratory system disease Diseases 0.000 description 1
- 238000009991 scouring Methods 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 229910052708 sodium Inorganic materials 0.000 description 1
- 239000011734 sodium Substances 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
- G01N27/62—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating the ionisation of gases, e.g. aerosols; by investigating electric discharges, e.g. emission of cathode
-
- 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/71—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light thermally excited
- G01N21/73—Systems in which the material investigated is excited whereby it emits light or causes a change in wavelength of the incident light thermally excited using plasma burners or torches
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/12—Measuring magnetic properties of articles or specimens of solids or fluids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
- G01R33/00—Arrangements or instruments for measuring magnetic variables
- G01R33/12—Measuring magnetic properties of articles or specimens of solids or fluids
- G01R33/1215—Measuring magnetisation; Particular magnetometers therefor
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
Landscapes
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Life Sciences & Earth Sciences (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Condensed Matter Physics & Semiconductors (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Evolutionary Biology (AREA)
- Bioinformatics & Computational Biology (AREA)
- Electrochemistry (AREA)
- Operations Research (AREA)
- Probability & Statistics with Applications (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Algebra (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Plasma & Fusion (AREA)
- Other Investigation Or Analysis Of Materials By Electrical Means (AREA)
Abstract
本发明公开了一种基于支持向量机和树叶磁学的大气重金属污染评估方法,包括以下步骤:S1:大气颗粒物样品采集;S2:树叶样品采集;S3:大气颗粒物重金属提取及分析;S4:树叶样品磁学参数测试;S5:支持向量机磁学模型构建;S6:模型误差评估;S7:模型实际应用。总之,本发明所提供的一种基于支持向量机和树叶磁学的大气重金属污染评估方法,搭建环磁学与重金属污染定量化评价的桥梁,拓展利用环境磁学评价大气污染的方法和技术;所构建的模型可用于城市大气颗粒物重金属污染水平的快速模拟和预测,进而识别、把握城市大气重金属污染状况和变化趋势,为大气重金属污染综合防控提供重要参考。
Description
技术领域
本发明涉及城市大气颗粒物重金属监测和污染研究技术领域,具体是涉及一种基于支持向量机和树叶磁学的大气重金属污染评估方法。
背景技术
大气颗粒物来源广泛、结构复杂、比表面积大,易富集多种污染物如重金属。研究表明,大气颗粒物中重金属与人类的呼吸系统疾病和心血管疾病的发生具有重要联系。我国的环境空气质量标准(GB3095-2012)将大气汞、砷、六价铬和镉列为省级政府机构建议制定和实施监测的环境空气质量指标。传统的大气重金属污染监测方法是利用常规收集装置收集一定粒径大气颗粒物,化学提取并使用原子吸收光谱仪或电感耦合等离子质谱等方法检测颗粒物中重金属浓度,工作量较大、耗时且费用昂贵,因此,选择高效、便捷的方法研究大气颗粒物中重金属污染具有重要的现实意义和实践价值。
支持向量机(Support Vector Machines,SVM)是由Vapnik开发的新一代基于统计学习理论的学习机器,能解决小样本、非线性、高维数和局部极小点等实际问题,遵循结构风险最小化原则,可以成功应用于分类、归类和时间序列预测等领域。Libsvm是台湾大学林智仁等开发设计的一个简单、快速、有效的支持向量机模式识别与回归的软件程序,它提供了编译好的基于Windows操作系统的执行文件以及相关的软件程序源代码。
环境磁学的原理是通过测量环境介质中物质的磁学特征,从中提取环境信息。与传统的化学分析方法相比,磁学测试具有样品用量少、简便快捷、灵敏度高、费用低和非破坏性等特点。大气颗粒物中的金属元素主要有钠、钾、镁、钙、铝、铁等地壳元素以及铅、锌、砷、镉、铜等主要来自人类污染的元素。城市大气颗粒物中磁性铁氧化物和铁的氢氧化物占铁元素总质量的10-70%,主要来自工业磨损、含铁杂质燃料燃烧、钢铁加工和机动车零部件磨损等。这些磁性颗粒与大气重金属的共源性使得利用磁性特征描述或量化大气重金属成为可能。
目前针对区域性大气污染时空分布的研究多以模拟分析为主,如数值模拟和统计预报。数值模拟通常需较详尽的污染源信息及时空分辨率很高的气象预报模式等,在实际中往往难以获取。大气污染物浓度的统计预报是通过分析污染物变化规律,建立典型因子与污染物浓度之间的定量或半定量关系,根据这些关系作浓度预报,无需构建复杂的数理模型,比较简便、易于实现,且在借助合适的输入变量及模型下,常可以获得较高的模拟预测精度。近年来常用的支持向量机模型可以解决小样本、非线性、高维数和局部极小点等实际问题,是基于支持统计学习理论的学习机器,已经成功应用于建立NOX、CO、SO2、O3、PM的联系及污染物预测。尽管目前我国很多城市大气颗粒物已基本实现自动、在线、连续监测,但由于大气重金属监测分析工作量较大、研究不够,还缺乏数据积累,并未开展系统的监测和分析工作。
城市常绿树木分布广泛,叶片比表面积大,生长周期长,可以通过叶片表面的蜡质层吸附或通过气孔直接吸收大气中的颗粒物,尤其是更容易吸附或吸收粒径小于10μm的细颗粒物,是颗粒物天然的“接收器”。因此,可利用植物叶片作为载体来反演短时间内大气颗粒物中重金属的污染状况,为其污染水平、时空分布和来源等方面的研究提供重要信息和依据。空气监测站通常监测的是实时颗粒物浓度,而利用植物叶片作为载体可以很好地反演短时间内大气颗粒物污染状况。近年来,国内外已有利用树叶的磁性特征来表征大气颗粒物污染的研究,但目前还没有利用树叶磁学参数构建非线性数学模型来模拟城市大气颗粒物中重金属浓度的报道。
因此,本发明采用支持向量机的非线性数学方法,以树叶磁学参数、常规大气污染物和气象数据为基础,构建大气颗粒物中重金属磁学评估模型,模型可用于模拟预测城市大范围内大气颗粒物中重金属的时空分布,为我国城市大气中重金属污染的快速监测、评估和风险管理提供支持。
发明内容
本发明解决的技术问题是:本发明针对现有技术的不足,提供一种基于支持向量机技术和树叶磁学参数的大气颗粒物重金属污染模拟、评估方法。
本发明的技术方案是:一种基于支持向量机和树叶磁学的大气重金属污染评估方法,其特征在于,包括以下步骤:
S1:大气颗粒物样品采集
使用中流量采样装置收集大气颗粒物PM2.5,采样介质使用高纯石英膜,采样前后滤膜均恒温恒湿48h并称重,所述恒温恒湿条件为温度20-25℃,湿度40-50%,按一定条件进行采样;
S2:树叶样品采集
在距离大气颗粒物采样点1km范围内调查常青灌木或乔木分布情况,以1-3种树木作为研究对象,选择树龄和高度相似、生长状况良好的树上成熟、健康无病害、无虫斑的树叶,采集新生的树枝上最成熟的叶片,每4天进行一次与颗粒物的同步采样,采样时佩戴一次性手套,避免其他物品直接接触树叶样品,与大气颗粒物同步采样,周期为1年;样品放入自封袋带回实验室放入冰箱冷藏室保存,以备测试;
S3:大气颗粒物重金属提取及分析
将大气颗粒物采样膜用陶瓷剪刀剪下1/2,并剪成条状碎片,置于消解管中,加入HF、HNO3和HCl进行提取,使用电感耦合等离子体质谱仪ICP-MS测试消解液中微量元素,使用电感耦合等离子体发射光谱仪ICP-OES测定消解液中大量金属元素;取等面积空白滤膜,按样品处理方法制成空白溶液并检测元素浓度;使用美国国家标准与技术研究院标准参考物质NIST SRM 1648a(城市颗粒物)进行加标回收,作为分析质量控制;
S4:树叶样品磁学参数测试
将树叶样品放入烘箱烘干,烘干后剪碎,称取2g样品用保鲜膜包裹放在圆柱状聚乙烯样品盒中,压实待测;利用KY-3S卡帕桥磁化率仪测量样品的体积磁化率κ,经密度校正得到样品的质量磁化率χ;使用Molspin交变退磁仪对样品进行退磁后,使用JR-6A旋转磁力仪测量样品的非磁滞剩磁ARM,并计算出无磁滞磁化率χARM。使用ASC Scientific ModelIM-10-30对样品加1000mT的磁场进行脉冲磁化后,使用JR-6A旋转磁力仪测试样品的等温剩磁作为饱和等温剩磁SIRM。其中,质量磁化率χ是环境介质中铁磁性矿物含量的粗略量度;非磁滞剩磁ARM提供了铁磁晶粒的磁畴信息,对单畴晶粒的反映尤为敏感,可反映磁性颗粒物的粒径组成;饱和等温剩磁剩磁SIRM反映了磁性矿物在样品中的富集程度,与环境介质中铁磁性矿物的类型和含量有关;
S5:支持向量机磁学模型构建
分别以“气象因子+污染物+树叶磁学参数”,以及“树叶磁学参数”作为输入参数,所述气象因子包括:温度、相对湿度、气压和风速;所述污染物包括:PM2.5、SO2、NO2和O3浓度;所述磁学参数包括:树叶的χ、χARM和SIRM;以PM2.5中重金属浓度作为输出层,建立大气颗粒物重金属污染磁学评估模型;
S6:模型误差评估
使用相关系数,平均绝对误差和均方根误差来评价模型性能及预测效果;相关系数用于反映变量之间的相关关系密切程度,用于测量模型的拟合性能;平均绝对误差是所有单个模拟值与实测值偏差绝对值的平均值,用于评估模型的预测效果;均方根误差是误差平方和的平均值的平方根,用于判定该模型的可靠性,其计算公式分别如下:
S7:模型实际应用
实际应用中,模型训练完毕后,可针对任意区域任意时间段,采集不同位置用于建模的树种的叶片并分析其磁学参数,然后将“树叶磁学参数”或者“气象因子+污染物+树叶磁学参数”作为输入因子,用上述已训练好的模型预测该区域、该时段该位置对应大气颗粒物重金属的浓度,进而获得重金属浓度的时空分布情况。
进一步的,所述S1中,采样方法具体为:按季节采样,每次采样持续4天,每天固定采集8-12小时,每个样品共采集32-48个小时,避免大风和雨雪等极端天气,采样为期1年,采用上述采样方法采集的样品代表性强,能有效评估一年四季大气颗粒物重金属的含量,能使得评估的准确性更加精确。
进一步的,所述S3中,测试消解液中微量元素有As、Cd、Co、Cr、Cu、Mn、Ni、Pb、Sb、Ti、TI和V,但不仅限于这十二种;测定消解液中大量金属元素有Al、Fe和Zn,但不仅限于这三种。
进一步的,所述S4中,使用陶瓷剪刀进行剪碎,通过使用陶瓷剪刀剪碎叶片,可以防止金属剪刀在剪碎时夹带金属成分影响结果的准确性。
进一步的,所述S4中,脉冲磁化的磁场强度为1000mT。
进一步的,所述S5中,随机选择80%的数据作为训练样本,且包含其中的极大值和极小值,剩下的20%作为验证样本,采用随机选择80%作为训练样本与20%的验证样本对比,对比充分,更易对比得到的评估结果,且评估结果更客观准确。
进一步的,所述方法中支持向量机模拟模型建立具体过程为:在MATLAB R2013a环境下安装Libsvm-3.21,选取支持向量机中常用的径向核函数,RBF核函数,epsilon-SVR模型进行模拟,通过交叉试验的方法选取模型最佳核函数g值及惩罚因子c值,建立重金属浓度的磁学模型,在至少100次成功建模之中,选择训练和验证模型中相关系数最大的一次作为模拟模型。
进一步的,支持向量机用于回归统计的方法具体如下:
设训练样本为(xi,yi),(i=1,…,l),建立如下线性回归方程:
其中,为非线性映射函数,ω为法向量,b为偏置量;
进而训练目标可表示为:
且满足以下约束条件:
其中,ε为不敏感损失函数,ε≥0,表示允许的最大误差,其取值的大小影响支持向量数目;C≥0为惩罚系数,它控制对超出误差样本的惩罚程度;ξi,为松弛变量,指示样本超出允许误差的程度,分别对应计算值小于实际值及计算值大于实际值两种情况;
为得到以上最优化问题(式(5)),引入了拉格朗日乘子法,得其对偶优化问题为:
满足以下约束条件:
其中,αi和是拉格朗日乘子;
对式(9)求解,对应的预测输入样本为xi,得到如下预测值:
所述方法采用常用的径向基核函数RBF作为核函数,其表达式为:
其中,δ2是核函数的方差。
与现有技术相比,本发明的有益效果如下:
1、本发明所提供的一种基于支持向量机和树叶磁学的大气重金属污染评估方法,搭建了环磁学与重金属污染定量化评价的桥梁,拓展了利用环境磁学评价大气污染的方法和技术。
2、鉴于大气颗粒物重金属污染特征复杂、分析费时费力,对于其污染的快速监测具有一定紧迫性。本发明所构建的模型可用于城市大气颗粒物重金属污染水平的快速模拟和预测,进而识别、把握城市大气重金属污染状况和变化趋势,为大气重金属污染综合防控提供重要参考。
附图说明
图1是本发明的流程框架图;
图2是磁学模型对PM2.5中As、Cd、Cr、Fe和Pb的模拟值与实测值的对比,(a)-(e)是模型1的模拟结果,(f)-(j)是模型2的模拟结果;
具体实施方式
下面结合具体实施方式来对本发明进行更进一步详细的说明,以更好地体现本发明的优势。
一种基于支持向量机和树叶磁学的大气重金属污染评估方法,包括以下步骤:
S1:大气颗粒物样品采集
近年来,长三角地区大气颗粒物污染比较严重。南京是江苏省省会,长三角地区的第二大城市,我国重要的综合性工业生产基地,同时是华东地区重要的交通枢纽。本发明选择位于南京市北部的南京大学仙林校区作为研究区域,进行大气颗粒物和树叶样品的采集及模型构建。仙林大学城是江苏省发展高等教育产业的重点地区,因距离南京城北、江北工业区较近,其空气质量受到一定影响。
PM2.5采样设备为XY-2200型大气综合采样器(流量100L/min),采样器采气口距离地面1.5m,采样膜选用高纯石英膜,采样时间为2015年12月4日到11月30日,为期一年。每次采样持续4天,每天8小时,采样时间段为7:00-15:00,每个样品共采集32个小时,大风和雨雪等极端天气暂停采样,最终每个季节获得21个、共计获得PM2.5样品84个。同步记录采样点附近南京大学空气质量自动监测站的实时气象数据和大气污染物数据。
S2:树叶样品采集
桂花是南京市最常见、分布最广泛的常绿树种之一,其叶片对大气颗粒物吸附能力较强,且树木高度及叶片大小很适宜进行本发明的采样及分析。在距离PM2.5采样点约0.5km范围内采集桂花树叶,采样时间与频率和大气颗粒物相同,每隔4天采集一次,选取树龄为3年、高度为2m左右的3棵树,在每棵树的不同方向分别收集4片树叶,采样高度为1.5m。选择成熟、健康无病害、无虫斑的树叶进行采样。采样时佩戴一次性手套,用陶瓷剪刀剪取树叶,避免其他物品直接接触树叶样品。样品采集后放入聚乙烯袋带回实验室。采样时间为期1年,共获得84个树叶样本。
S3:大气颗粒物中重金属提取及分析
每张滤膜的1/2用于分析金属元素含量。将大气颗粒物采样膜用陶瓷剪刀剪下1/2,并剪成条状碎片,置于消解管中,加入HF、HNO3和HCl进行提取;As、Cd、Co、Cr、Cu、Mn、Ni、Pb、Sb、Ti、TI和V浓度使用电感耦合等离子体质谱(ICP-MS,Perkin Elmer SCIEX,Elan9000,Norway)测定,最低检出限的0.01μg/L;Al、Fe和Zn使用电感耦合等离子体发射光谱仪(ICP-OES,Perkin Elmer SCIEX,Optima 5300DV,Norway)测定,最低检出限为0.001mg/L。以美国国家标准与技术研究院标准参考物质:NIST SRM 1648a(城市颗粒物)为质控标样,加标回收作为分析质量控制,绝大部分元素加标回收率为106%,符合质控要求。
PM2.5中绝大部分金属元素含量在冬季的浓度最高,只有Co和Cr在春季最高;大部分金属元素含量冬季和秋季高于夏秋季。具体地,金属元素浓度季节分布情况如下(单位ng/m3):春季:Al(1297)>Fe(930.7)>Zn(199.6)>Ti(50.66)>Ni(46.65)>Mn(44.28)>Cr(39.32)>Pb(39.20)>Cu(34.35)>As(7.392)>V(5.505)>Sb(3.52)>Cd(1.308)>Co(0.496)>TI(0.411);夏季:Al(1177)>Fe(445.4)>Zn(196.8)>Ni(36.72)>Ti(25.89)>Cr(25.67)>Mn(23.83)>Cu(18.05)>Pb(14.69)>As(4.946)>V(4.249)>Sb(2.916)>Cd(1.276)>TI(0.372)>Co(0.284);秋季:Al(1186)>Fe(354.3)>Zn(214.4)>Ni(39.12)>Pb(29.51)>Mn(28.88)>Ti(28.78)>Cr(27.26)>Cu(17.53)>V(3.875)>As(3.848)>Sb(3.145)>Cd(1.477)>TI(0.279)>Co(0.223);冬季:Al(1550)>Fe(1056)>Zn(275.7)>Pb(73.77)>Mn(56.79)>Ti(51.59)>Ni(47.7)>Cu(43.28)>Cr(36.85)>As(8.754)>V(5.9)>Sb(4.742)>Cd(2.345)>TI(0.697)>Co(0.457)。
S4:树叶样品的磁学参数
将采集到的桂花树叶样品放入烘箱烘干,温度为55℃,烘干后用陶瓷剪刀剪碎,称取约2g样品用保鲜膜包裹置于用于磁学参数测试的圆柱状聚乙烯样品盒(体积为11.15cm3)中,压实待测;在875Hz的频率下,使用卡帕桥磁化率仪KLY-3测试低频磁化率χlf;使用Molspin交变退磁仪进行退磁后,将样品置于0.04mT的直流偏置场中和100mT的交变磁场峰值下获得非磁滞剩磁,再利用JR-6A旋转磁力仪测量样品的非磁滞剩磁(ARM),计算成无磁滞磁化率(χARM);使用ASC Scientific Model IM-10-30对样品加1000mT的磁场进行脉冲磁化后,利用JR-6A旋转磁力仪测量样品的等温剩磁(IRM1T),作为样品的饱和等温剩磁(SIRM)。每个季节树叶的磁学参数如表1所示,χlf:冬季>春季>秋季>夏季;χARM值冬季最高,夏季最低;SIRM:冬季>春季>秋季>夏季。
表1桂花树叶磁学参数特征
S5:支持向量机模型构建
将大气颗粒物中重金属浓度与树叶磁学参数、污染物浓度、气象数据进行皮尔逊相关性分析,结果显示,PM2.5中金属浓度和树叶的χlf、χARM、SIRM值显著正相关,说明树叶所吸附的颗粒物中的磁性矿物与大气重金属之间具有一定的共源性及内在联系。所有金属元素浓度与PM2.5浓度均呈现正相关,Pb与PM2.5浓度相关系数最高,为0.781,Co与PM2.5浓度相关系数最低,为0.229;绝大部分金属元素浓度与O3、温度、相对湿度、风速呈显著负相关,与SO2、NO2和大气压正相关。大气重金属与大气污染物之间具有较强的共源性如工业排放、交通污染和道路扬尘等。风速对大气污染物有一定稀释作用,一般来说,温度越高,大气压越低,空气对流运动越明显,大气颗粒物扩散速率越快,其浓度也同时有所下降;空气湿度对大气颗粒物最显著的影响是对其的冲刷作用。
在分析了重金属与输入因子之间的内在关系之后,以金属浓度作为输出层,将树叶的磁学参数(χlf、χARM和SIRM)、污染物浓度(PM2.5、SO2、NO2和O3)和气象数据(温度、相对湿度、压力和风速)作为输入层作为本案例的模型1;另外,为提高模型的可应用性,考虑到采样点如不能获取气象数据和污染物的情况,再直接将“树叶的磁学参数(χlf、χARM和SIRM)”作为输入因子构建本案例的模型2。随机选择80%的数据作为训练样本,剩下的20%作为验证样本,为了确保模型的可靠性,训练集中包含了不同季节数据的极值。应用MATLAB R2013a和libsvm-3.21建立支持向量机模型,输入的参数首先利用核函数的径向基函数映射到高维特征空间,然后在特征空间中进行线性回归。
S6:模型误差评估
模型的R、MAE、RMSE值如表2所示。模型1(输入因子为树叶的磁学参数、污染物浓度和气象数据)的训练和验证模型的R值均大于0.6,误差在可接受范围内,其中Fe、Mn和Pb训练和验证模型的R值均大于0.8。对于模型2(树叶的磁学参数),只有树叶磁学参数作为输入因子时,训练和验证模型的R值也均大于0.6,误差都在可接受范围内,且训练和验证模型As、Cd、Fe、Mn和Pb的R值均大于0.7。两类模型对Fe和Pb的模拟效果较好,误差也较低,而对V、Co和Sb的预测效果相对较弱。模型对不同金属元素的模拟结果与研究区域、采样时间及来源有关,具体原因仍需进一步验证。总体上,模型1和模型2均可对大气颗粒物重金属浓度进行较为理想的预测。以As、Cd、Cr、Fe和Pb为例,树叶磁学模型对金属元素浓度模拟值与监测值的比较如附图2所示。
表2磁学模型对PM2.5中重金属浓度模拟的相关性系数(R)、平均绝对误差(MAE)和均方根误差(RMSE)
S7:南京大气颗粒物重金属浓度分布预测
上述模型实际应用时,采用网格布点法在南京市划定5km×5km的网格,每个网格中选择2个共50个居民小区作为树叶采样点。采样前一周及当天均晴朗无雨,在每个小区较为开阔的地带选择树龄、高度相似及生长状况良好的桂花树3棵,于每棵树的不同方向采摘4片共计10片左右成熟、健康、无虫斑的树叶,放入自封袋带回实验室。将树叶烘干并测试其磁学参数χlf、χARM和SIRM。将三个磁学参数作为输入因子代入上述支持向量机模型2中,即可模拟出50个点位PM2.5中重金属的浓度值。结果显示,距南钢集团、金陵石化或扬子石化等地较近或距交通密集区较近的小区,树叶磁学参数值及模拟输出的PM2.5中重金属浓度均较高,而附近交通密度较低又无明显工业源影响的小区所模拟出的重金属浓度较低。本发明所构建的模型可用于模拟预测城市大范围内大气颗粒物中重金属的时空分布,从而为我国城市大气中重金属污染的快速监测和风险管理提供参考和支持。
最后应说明的是:以上实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述实施例所记载的技术方案进行修改,或者对其中部分技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明实施例技术方案的精神和范围。
Claims (8)
1.一种基于支持向量机和树叶磁学的大气重金属污染评估方法,其特征在于,包括以下步骤:
S1:大气颗粒物样品采集
使用中流量采样装置收集大气颗粒物PM2.5,采样介质使用高纯石英膜,采样前后滤膜均恒温恒湿48h并称重,所述恒温恒湿条件为温度20-25℃,湿度40-50%;按一定条件进行采样;
S2:树叶样品采集
在距离大气颗粒物采样点1km范围内调查常青灌木或乔木分布情况,以1-3种树木作为研究对象,采集叶片,每4天进行一次与颗粒物的同步采样,采样时佩戴一次性手套,避免其他物品直接接触树叶样品,与大气颗粒物同步采样,周期为1年;样品放入自封袋带回实验室放入冰箱冷藏室保存,以备测试;
S3:大气颗粒物重金属提取及分析
将大气颗粒物采样膜用陶瓷剪刀剪下1/2,并剪成条状碎片,置于消解管中,加入HF、HNO3和HCl进行提取;测试消解液中微量元素和测定消解液中大量金属元素;取等面积空白滤膜,按样品处理方法制成空白溶液并检测元素浓度;
S4:树叶样品磁学参数测试
将树叶样品放入烘箱烘干,烘干后剪碎,称取2g样品用保鲜膜包裹放在圆柱状聚乙烯样品盒中,压实待测;测量样品的体积磁化率κ,经密度校正得到样品的质量磁化率χ;进行退磁后,测量样品的非磁滞剩磁ARM,并计算出无磁滞磁化率χARM;对样品进行脉冲磁化后,测试样品的等温剩磁作为饱和等温剩磁SIRM;
S5:支持向量机磁学模型构建
分别以“气象因子+污染物+树叶磁学参数”,以及“树叶磁学参数”作为输入参数,所述气象因子包括:温度、相对湿度、气压和风速;所述污染物包括:PM2.5、SO2、NO2和O3浓度;所述磁学参数包括:树叶的χ、χARM和SIRM;以PM2.5中重金属浓度作为输出层,建立大气颗粒物重金属污染磁学评估模型;
S6:模型误差评估
使用相关系数,平均绝对误差和均方根误差来评价模型性能及预测效果;相关系数用于反映变量之间的相关关系密切程度,用于测量模型的拟合性能;平均绝对误差是所有单个模拟值与实测值偏差绝对值的平均值,用于评估模型的预测效果;均方根误差是误差平方和的平均值的平方根,用于判定该模型的可靠性,其计算公式分别如下:
S7:模型实际应用
实际应用中,模型训练完毕后,可针对任意区域任意时间段,采集不同位置用于建模的树种的叶片并分析其磁学参数,然后将“树叶磁学参数”或者“气象因子+污染物+树叶磁学参数”作为输入因子,用上述已训练好的模型预测该区域、该时段该位置对应大气颗粒物重金属的浓度,进而获得重金属浓度的时空分布情况。
2.根据权利要求1所述的一种基于支持向量机和树叶磁学的大气重金属污染评估方法,其特征在于,所述S1中,采样方法具体为:按季节采样,每次采样持续4天,每天固定采集8-12小时,每个样品共采集32-48个小时,避免大风和雨雪等极端天气,采样为期1年。
3.根据权利要求1所述的一种基于支持向量机和树叶磁学的大气重金属污染评估方法,其特征在于,所述S3中,测试消解液中微量元素有As、Cd、Co、Cr、Cu、Mn、Ni、Pb、Sb、Ti、TI和V;测定消解液中大量金属元素有Al、Fe和Zn。
4.根据权利要求1所述的一种基于支持向量机和树叶磁学的大气重金属污染评估方法,其特征在于,所述S3中,测试消解液中微量元素有十二种,所述微量元素分别是As、Cd、Co、Cr、Cu、Mn、Ni、Pb、Sb、Ti、TI和V。
5.根据权利要求1所述的一种基于支持向量机和树叶磁学的大气重金属污染评估方法,其特征在于,所述S4中,使用陶瓷剪刀进行剪碎。
6.根据权利要求1所述的一种基于支持向量机和树叶磁学的大气重金属污染评估方法,其特征在于,所述S4中,脉冲磁化的磁场强度为1000mT。
7.根据权利要求1所述的一种基于支持向量机和树叶磁学的大气重金属污染评估方法,其特征在于,所述S5中,随机选择80%的数据作为训练样本,且包含其中的极大值和极小值,剩下的20%作为验证样本。
8.根据权利要求1所述的一种基于支持向量机和树叶磁学的大气重金属污染评估方法,其特征在于,支持向量机用于回归统计的方法具体如下:
设训练样本为(xi,yi),(i=1,…,l),建立如下线性回归方程:
其中,为非线性映射函数,ω为法向量,b为偏置量;
进而训练目标可表示为:
且满足以下约束条件:
其中,ε为不敏感损失函数,ε≥0;C≥0为惩罚系数;ζi,为松弛变量;
为得到以上最优化问题(式(5)),引入了拉格朗日乘子法,得其对偶优化问题为:
满足以下约束条件:
其中,αi和是拉格朗日乘子;
对式(9)求解,对应的预测输入样本为xi,得到如下预测值:
所述方法采用常用的径向基核函数RBF作为核函数,其表达式为:
其中,δ2是核函数的方差。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711380579.9A CN108169316A (zh) | 2017-12-20 | 2017-12-20 | 基于支持向量机和树叶磁学的大气重金属污染评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711380579.9A CN108169316A (zh) | 2017-12-20 | 2017-12-20 | 基于支持向量机和树叶磁学的大气重金属污染评估方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108169316A true CN108169316A (zh) | 2018-06-15 |
Family
ID=62522801
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711380579.9A Pending CN108169316A (zh) | 2017-12-20 | 2017-12-20 | 基于支持向量机和树叶磁学的大气重金属污染评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108169316A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110333556A (zh) * | 2019-06-03 | 2019-10-15 | 深圳中兴网信科技有限公司 | 空气质量预测方法、装置、计算机设备及可读存储介质 |
CN110411968A (zh) * | 2019-09-03 | 2019-11-05 | 南京信息工程大学 | 一种大气铬元素浓度检测方法 |
CN111768861A (zh) * | 2019-11-27 | 2020-10-13 | 复旦大学 | 膳食重金属水平评估模型的构建方法、评估方法和系统 |
CN115796034A (zh) * | 2022-12-01 | 2023-03-14 | 信阳师范学院 | 一种基于机器学习和数值模式的道路扬尘贡献评估系统和方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6961719B1 (en) * | 2002-01-07 | 2005-11-01 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Hybrid neural network and support vector machine method for optimization |
CN102661889A (zh) * | 2012-04-27 | 2012-09-12 | 上海大学 | 大气颗粒物中分离磁性重金属元素的方法 |
CN103175713A (zh) * | 2013-03-01 | 2013-06-26 | 中国环境科学研究院 | 一种适用于大气干湿沉降中重金属分析的样品采集及提取方法 |
CN103913508A (zh) * | 2014-01-14 | 2014-07-09 | 常州千帆环保科技有限公司 | 一种快速检测土壤中重金属含量及空间分布的方法与装置 |
CN104198512A (zh) * | 2014-08-18 | 2014-12-10 | 北京农业质量标准与检测技术研究中心 | 基于支持向量机的x射线荧光光谱分析方法及装置 |
CN107121650A (zh) * | 2017-05-04 | 2017-09-01 | 南京大学 | 基于支持向量机技术的大气颗粒物重金属磁学评估方法 |
-
2017
- 2017-12-20 CN CN201711380579.9A patent/CN108169316A/zh active Pending
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6961719B1 (en) * | 2002-01-07 | 2005-11-01 | The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration | Hybrid neural network and support vector machine method for optimization |
CN102661889A (zh) * | 2012-04-27 | 2012-09-12 | 上海大学 | 大气颗粒物中分离磁性重金属元素的方法 |
CN103175713A (zh) * | 2013-03-01 | 2013-06-26 | 中国环境科学研究院 | 一种适用于大气干湿沉降中重金属分析的样品采集及提取方法 |
CN103913508A (zh) * | 2014-01-14 | 2014-07-09 | 常州千帆环保科技有限公司 | 一种快速检测土壤中重金属含量及空间分布的方法与装置 |
CN104198512A (zh) * | 2014-08-18 | 2014-12-10 | 北京农业质量标准与检测技术研究中心 | 基于支持向量机的x射线荧光光谱分析方法及装置 |
CN107121650A (zh) * | 2017-05-04 | 2017-09-01 | 南京大学 | 基于支持向量机技术的大气颗粒物重金属磁学评估方法 |
Non-Patent Citations (2)
Title |
---|
冷湘梓: "基于非线性数学方法的PM2.5中重金属浓度模拟", 《中国环境科学》 * |
王呈: "南京市树叶附尘对大气重金属污染的磁学响应", 《中国环境科学》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110333556A (zh) * | 2019-06-03 | 2019-10-15 | 深圳中兴网信科技有限公司 | 空气质量预测方法、装置、计算机设备及可读存储介质 |
CN110411968A (zh) * | 2019-09-03 | 2019-11-05 | 南京信息工程大学 | 一种大气铬元素浓度检测方法 |
CN111768861A (zh) * | 2019-11-27 | 2020-10-13 | 复旦大学 | 膳食重金属水平评估模型的构建方法、评估方法和系统 |
CN111768861B (zh) * | 2019-11-27 | 2024-05-17 | 复旦大学 | 膳食重金属水平评估模型的构建方法、评估方法和系统 |
CN115796034A (zh) * | 2022-12-01 | 2023-03-14 | 信阳师范学院 | 一种基于机器学习和数值模式的道路扬尘贡献评估系统和方法 |
CN115796034B (zh) * | 2022-12-01 | 2024-01-30 | 信阳师范学院 | 基于机器学习和数值模式的道路扬尘贡献评估系统和方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Chang et al. | First long-term and near real-time measurement of trace elements in China's urban atmosphere: temporal variability, source apportionment and precipitation effect | |
CN107300550A (zh) | 一种基于bp神经网络模型预测大气重金属浓度的方法 | |
Pan et al. | Wet and dry deposition of atmospheric nitrogen at ten sites in Northern China | |
CN108169316A (zh) | 基于支持向量机和树叶磁学的大气重金属污染评估方法 | |
Gao et al. | Insights into the chemical characterization and sources of PM2. 5 in Beijing at a 1-h time resolution | |
Bove et al. | An integrated PM2. 5 source apportionment study: positive matrix factorisation vs. the chemical transport model CAMx | |
Tan et al. | Source of atmospheric heavy metals in winter in Foshan, China | |
Choi et al. | Source apportionment of PM2. 5 at the coastal area in Korea | |
Kim et al. | The chemical composition of fine and coarse particles in relation with the Asian Dust events | |
Hedberg et al. | Source contributions to PM10 and arsenic concentrations in Central Chile using positive matrix factorization | |
Furusjö et al. | PM10 source characterization at urban and highway roadside locations | |
Khuzestani et al. | Quantification of the sources of long-range transport of PM2. 5 pollution in the Ordos region, Inner Mongolia, China | |
Engvall et al. | Changes in aerosol properties during spring-summer period in the Arctic troposphere | |
CN107121650A (zh) | 基于支持向量机技术的大气颗粒物重金属磁学评估方法 | |
Zhang et al. | Continuous monitoring, compositions analysis and the implication of regional transport for submicron and fine aerosols in Beijing, China | |
CN109142650A (zh) | 一种蔬菜中镉含量预测模型的建模方法及其应用 | |
CN103913508A (zh) | 一种快速检测土壤中重金属含量及空间分布的方法与装置 | |
Brignole et al. | Chemical and magnetic analyses on tree bark as an effective tool for biomonitoring: A case study in Lisbon (Portugal) | |
CN108072597A (zh) | 基于支持向量机和树叶磁学的大气颗粒物浓度监测方法 | |
Huang et al. | Study on the seasonal variation and source apportionment of PM10 in Harbin, China | |
CN110503348B (zh) | 一种基于位置匹配的个体空气污染暴露模拟测量方法 | |
Zabalza et al. | Concentration and sources of PM 10 and its constituents in Alsasua, Spain | |
Su et al. | Is biomagnetic leaf monitoring still an effective method for monitoring the heavy metal pollution of atmospheric particulate matter in clean cities? | |
Seto et al. | Precipitation chemistry in western Japan: Its relationship to meteorological parameters | |
CN110907318B (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: 20180615 |