CN108956881A - SDABB Source Apportionment, system, terminal device and computer readable storage medium - Google Patents
SDABB Source Apportionment, system, terminal device and computer readable storage medium Download PDFInfo
- Publication number
- CN108956881A CN108956881A CN201810778540.0A CN201810778540A CN108956881A CN 108956881 A CN108956881 A CN 108956881A CN 201810778540 A CN201810778540 A CN 201810778540A CN 108956881 A CN108956881 A CN 108956881A
- Authority
- CN
- China
- Prior art keywords
- source
- class
- component
- matrix
- sdabb
- 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
- 239000002245 particle Substances 0.000 claims abstract description 94
- 239000011159 matrix material Substances 0.000 claims abstract description 78
- 238000001228 spectrum Methods 0.000 claims abstract description 45
- 239000008277 atmospheric particulate matter Substances 0.000 claims abstract description 39
- 238000004458 analytical method Methods 0.000 claims abstract description 32
- 238000000556 factor analysis Methods 0.000 claims abstract description 29
- 238000000034 method Methods 0.000 claims abstract description 12
- 239000013618 particulate matter Substances 0.000 claims abstract description 12
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 11
- 238000004590 computer program Methods 0.000 claims description 17
- 238000011109 contamination Methods 0.000 claims description 10
- 239000005416 organic matter Substances 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 8
- 229910052710 silicon Inorganic materials 0.000 claims description 8
- 229910002651 NO3 Inorganic materials 0.000 claims description 7
- NHNBFGGVMKEFGY-UHFFFAOYSA-N Nitrate Chemical compound [O-][N+]([O-])=O NHNBFGGVMKEFGY-UHFFFAOYSA-N 0.000 claims description 7
- 229910052782 aluminium Inorganic materials 0.000 claims description 7
- 229910052791 calcium Inorganic materials 0.000 claims description 6
- 229910052742 iron Inorganic materials 0.000 claims description 6
- 229910052719 titanium Inorganic materials 0.000 claims description 6
- 239000004568 cement Substances 0.000 claims description 5
- 238000007405 data analysis Methods 0.000 claims description 5
- 238000013480 data collection Methods 0.000 claims description 5
- 239000002689 soil Substances 0.000 claims description 5
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 4
- QAOWNCQODCNURD-UHFFFAOYSA-L Sulfate Chemical compound [O-]S([O-])(=O)=O QAOWNCQODCNURD-UHFFFAOYSA-L 0.000 claims description 4
- 229910052793 cadmium Inorganic materials 0.000 claims description 4
- 229910052799 carbon Inorganic materials 0.000 claims description 4
- 229910052804 chromium Inorganic materials 0.000 claims description 4
- KNHUKKLJHYUCFP-UHFFFAOYSA-N clofibrate Chemical compound CCOC(=O)C(C)(C)OC1=CC=C(Cl)C=C1 KNHUKKLJHYUCFP-UHFFFAOYSA-N 0.000 claims description 4
- 238000002485 combustion reaction Methods 0.000 claims description 4
- 239000000470 constituent Substances 0.000 claims description 4
- 229910052802 copper Inorganic materials 0.000 claims description 4
- 238000010790 dilution Methods 0.000 claims description 4
- 239000012895 dilution Substances 0.000 claims description 4
- 239000000428 dust Substances 0.000 claims description 4
- 239000000284 extract Substances 0.000 claims description 4
- 230000004907 flux Effects 0.000 claims description 4
- 238000002290 gas chromatography-mass spectrometry Methods 0.000 claims description 4
- 238000001095 inductively coupled plasma mass spectrometry Methods 0.000 claims description 4
- 238000002354 inductively-coupled plasma atomic emission spectroscopy Methods 0.000 claims description 4
- 238000004255 ion exchange chromatography Methods 0.000 claims description 4
- 229910052745 lead Inorganic materials 0.000 claims description 4
- 229910052749 magnesium Inorganic materials 0.000 claims description 4
- 229910052748 manganese Inorganic materials 0.000 claims description 4
- 229910052753 mercury Inorganic materials 0.000 claims description 4
- 229910052759 nickel Inorganic materials 0.000 claims description 4
- 229910052700 potassium Inorganic materials 0.000 claims description 4
- 229910052708 sodium Inorganic materials 0.000 claims description 4
- 229910052720 vanadium Inorganic materials 0.000 claims description 4
- 238000004846 x-ray emission Methods 0.000 claims description 4
- 229910052725 zinc Inorganic materials 0.000 claims description 4
- 239000002028 Biomass Substances 0.000 claims description 3
- 125000005575 polycyclic aromatic hydrocarbon group Chemical group 0.000 claims description 3
- TWNIBLMWSKIRAT-VFUOTHLCSA-N levoglucosan Chemical compound O[C@@H]1[C@@H](O)[C@H](O)[C@H]2CO[C@@H]1O2 TWNIBLMWSKIRAT-VFUOTHLCSA-N 0.000 claims description 2
- 239000011236 particulate material Substances 0.000 abstract description 5
- 230000009286 beneficial effect Effects 0.000 abstract 1
- 239000000243 solution Substances 0.000 description 6
- 238000004364 calculation method Methods 0.000 description 3
- 229920001503 Glucan Polymers 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- FGIUAXJPYTZDNR-UHFFFAOYSA-N potassium nitrate Chemical compound [K+].[O-][N+]([O-])=O FGIUAXJPYTZDNR-UHFFFAOYSA-N 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Classifications
-
- 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/0004—Gaseous mixtures, e.g. polluted air
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N15/00—Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
- G01N15/02—Investigating particle size or size distribution
-
- 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/0004—Gaseous mixtures, e.g. polluted air
- G01N33/0009—General constructional details of gas analysers, e.g. portable test equipment
- G01N33/0062—General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method or the display, e.g. intermittent measurement or digital display
-
- 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/0004—Gaseous mixtures, e.g. polluted air
- G01N33/0009—General constructional details of gas analysers, e.g. portable test equipment
- G01N33/0062—General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method or the display, e.g. intermittent measurement or digital display
- G01N33/0068—General constructional details of gas analysers, e.g. portable test equipment concerning the measuring method or the display, e.g. intermittent measurement or digital display using a computer specifically programmed
-
- 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/16—Matrix or vector computation, e.g. matrix-matrix or matrix-vector multiplication, matrix factorization
-
- 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
- Y02A50/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE in human health protection, e.g. against extreme weather
- Y02A50/20—Air quality improvement or preservation, e.g. vehicle emission control or emission reduction by using catalytic converters
Landscapes
- Engineering & Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Biochemistry (AREA)
- Analytical Chemistry (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Medicinal Chemistry (AREA)
- Mathematical Physics (AREA)
- Food Science & Technology (AREA)
- Combustion & Propulsion (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Dispersion Chemistry (AREA)
- Analysing Materials By The Use Of Radiation (AREA)
Abstract
The present invention relates to SDABB Source Apportionment, system, terminal device and computer readable storage mediums, and wherein method includes the following steps, step 1: collection analysis Atmospheric Particulate Matter source class, identify the mark component in the class of Atmospheric Particulate Matter source；Step 2: the particle diameter distribution of mark component is calculated；Step 3: acquiring and analyzes the component in ambient particle object pollution sources class；Step 4: " multi-source spectrum matrix nesting parallel factor analysis " three-dimensional Factor Analysis Model is utilized to calculate ambient particle object concentration of component matrix, and the particle diameter distribution of the source class obtained in step 2 mark component is included in algorithm, source spectrum matrix and source contribution matrix are obtained, identify source class and quantifies the contribution of source class.Beneficial effects of the present invention: component particle diameter distribution will be identified in airborne particulate material resource class and is included in three-dimensional factorial analysis Source Apportionment, obtained result more meets the physical significance of particulate matter source particle diameter distribution, improves the accuracy of source resolution result.
Description
Technical field
The invention belongs to Particulate Pollution origin analysis fields more particularly to SDABB Source Apportionment, system, terminal to set
Standby and computer readable storage medium.
Background technique
Gray haze pollution is one of the main atmospheric problem that current China city and region are faced, and shows generation frequency
It is secondary it is high, at the big feature of haze region area.Particulate matter origin analysis is science, effectively carries out the basic and preceding of gray haze prevention and cure of pollution
It mentions, provides indispensable scientific basis to formulate urban atmosphere Particulate Pollution control way.
PMF method and ABB method are the particulate matter Source Apportionments delivered, and PMF method is suitable for single particles
The particle size of material resource parsing, atmospheric environment is complicated, and PMF does not consider relationship between partial size, is easy to appear and do not meet actual particle size
The result of distribution；ABB is not included in the particle diameter distribution index of source class mark component, cannot reflect particle swarm size distribution feelings comprehensively
Condition.
Summary of the invention
To solve the above problems, the present invention a kind of SDABB Source Apportionment, system, terminal device and computer are provided can
Storage medium is read, wherein SDABB (Abb factor analyis method based on size distribution of
Markers in source profies) it is based on the three-dimensional Factor minute for identifying component particle diameter distribution in airborne particulate material resource class
Analysis, can identify particulate matter primary pollution source class, and each source contribution of quantitative estimation, the application will identify in airborne particulate material resource class
Component particle diameter distribution is included in three-dimensional factorial analysis Source Apportionment, and obtained result more meets the object of particulate matter source particle diameter distribution
Meaning is managed, the accuracy of source resolution result is improved.
The technical solution SDABB Source Apportionment of invention, it is characterised in that the following steps are included:
Step 1: acquiring and analyzes Atmospheric Particulate Matter source class, obtains the number of components in the class of Atmospheric Particulate Matter source
According to；
Step 2: according to the component data of the Atmospheric Particulate Matter source class of acquisition, Atmospheric Particulate Matter source class is identified
In mark component and calculate mark component particle diameter distribution；
Step 3: acquiring and analyzes the component in ambient particle object pollution sources class；
Step 4: " multi-source spectrum matrix nesting parallel factor analysis " three-dimensional Factor Analysis Model is utilized to calculate ambient particle object
Concentration of component matrix, and the particle diameter distribution of the pollution sources class obtained in step 2 mark component is included in algorithm, obtain source spectrum
Matrix and source contribution matrix, identification pollution sources class and the contribution of quantitative contamination source class.
Further, Atmospheric particulates are obtained using the acquisition of Atmospheric particulates source sampling device in the step 1；
Preferably, Atmospheric particulates source sampling device is dilution tunnel sampler, ELPI sampler or settling flux sampler.
Further, pollution sources class includes city raised dust, coal-fired source, motor vehicle source, sulfate source, nitre in the step 1
Hydrochlorate source, biomass combustion source, soil source, building building cement source.
Further, component includes element class in the step 1: Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Fe,
Ni, Cu, Cd, Zn, Hg, Pb etc.；Ionic species: Cl-、NO3 -、SO4 2-、NH4 +；EC/OC class: EC, OC；Organic matter species: the organic matter leaves of pulse plants
Alkane polycyclic aromatic hydrocarbon, levoglucosan.
Preferably, element class is measured using ICP-MS, ICP-AES or XRF, and ionic species are measured using ion chromatography, EC/OC
Class is measured using carbon component analysis instrument, and organic species are measured with GC-MS.
Further, the particle diameter distribution that component is identified in the step 2 is mark component in PM2.5Middle content and the mark
Component is known in PM10The ratio of middle content.
" multi-source spectrum matrix nesting parallel factor " analyzing three-dimensional Factor Analysis Model is to more partial size three dimensional receptors in step 4
Data carry out factorial analysis and obtain identical sources tribute to the particulate matter receptor of each partial size during Synchronization Analysis three-dimensional data
Trend matrix is offered, each partial size extracts respective derived components spectrum matrix respectively.
Further, the multi-source spectrum matrix nesting parallel factor analysis " three-dimensional Factor Analysis Model:
In formula, xijkIt is the concentration of j-th of component in i-th of sample of kth different particle types object, unit μ g/m3；
aihIt is the contribution for simulating h-th of source class being calculated to i-th of sample, unit μ g/m3；
bjhkIt is the content of j-th of component in p-th of source constituents spectrum of kth different particle types object, unit g/g；
ejhkIt is the residual error of j-th of component in i-th of sample of kth different particle types object, unit μ g/m3。
Further, object concentration of component matrix is as follows up to formula in step 4: X1=c1AB1
X2=c2AB2
…
Xz=czABz
In formula, XzIt is the acceptor density data matrix of z-th of particles things；
A is the identical source contribution value matrix of each partial size；
BzIt is the derived components spectrum matrix of z-th of particles things；
czIt is scale factor.
SDABB source resolution system, it is characterised in that comprise the following modules:
Particulate count obtains atmosphere for acquiring and analyzing Atmospheric Particulate Matter source class according to acquisition and analysis module
The module of component data in grain object pollution sources class；
Component particle diameter distribution module is identified, for the component data according to the Atmospheric Particulate Matter source class of acquisition, identification
Mark component in the class of Atmospheric Particulate Matter source and calculate mark component particle diameter distribution module；
Ambient particle object data collection and analysis module, for acquiring and analyzing the component in ambient particle object pollution sources class
Module；
The identification of source class and source class contribute quantitative module, utilize " multi-source spectrum matrix nesting parallel factor analysis " three-dimensional Factor minute
It analyses model and calculates ambient particle object concentration of component matrix, and the pollution sources class obtained in step 2 is identified to the particle diameter distribution of component
It is included in algorithm, obtains the module of source spectrum matrix and source contribution matrix, identification pollution sources class and the contribution of quantitative contamination source class.
SDABB source resolution terminal device, including memory, processor and storage are in the memory and can be described
The computer program run on processor, it is characterised in that described in the processor is realized when executing the computer program
The step of SDABB Source Apportionment.
Store the computer readable storage medium of SDABB line source analysis program, the computer-readable recording medium storage
There is computer program, it is characterised in that the computer program realizes the step of the source the SDABB solution method when being executed by processor
Suddenly.
The medicine have the advantages that being included in three-dimensional factorial analysis source for component particle diameter distribution is identified in airborne particulate material resource class
Analytic method, obtained result more meet the physical significance of particulate matter source particle diameter distribution, improve the accuracy of source resolution result.
Detailed description of the invention
Fig. 1 is the flow chart of the SDABB Source Apportionment of the embodiment of the present invention 1 and 2.
Fig. 2 is the structural block diagram of the SDABB source resolution system of the embodiment of the present invention 1 and 2.
Fig. 3 is the PM that the analysis of the embodiment of the present invention 2 obtains2.5Source spectrum.
Fig. 4 is the PM that the analysis of the embodiment of the present invention 2 obtains10Source spectrum.
Fig. 5 is every provenance class mark component of the embodiment of the present invention 2 in PM2.5Middle content in PM10The ratio of middle content
Datagram.
Fig. 6 is the PM for acquiring and analyzing in the embodiment of the present invention 2 in ambient particle object2.5Component datagram.
Fig. 7 is the PM for acquiring and analyzing in the embodiment of the present invention 2 in ambient particle object10Component datagram.
Fig. 8 is that source spectrum matrix data figure is obtained in the embodiment of the present invention 2.
Fig. 9 is that source contribution matrix data figure is obtained in the embodiment of the present invention 2.
Figure 10 is the identification source class and PM of the embodiment of the present invention 22.5Quantitative source class contribution data figure.
Figure 11 is the identification source class and PM of the embodiment of the present invention 210Quantitative source class contribution data figure.
Specific embodiment
It explains with reference to the accompanying drawing to a specific embodiment of the invention.
The technical solution of invention: SDABB Source Apportionment, comprising the following steps:
Step 1: acquiring and analyzes acquisition and analyzes Atmospheric Particulate Matter source class, obtains Atmospheric Particulate Matter source class
In component data；
Step 2: according to the component data of the Atmospheric Particulate Matter source class of acquisition, Atmospheric Particulate Matter source class is identified
In mark component and calculate mark component particle diameter distribution；
Step 3: acquiring and analyzes the component in ambient particle object pollution sources class；
Step 4: " multi-source spectrum matrix nesting parallel factor analysis " three-dimensional Factor Analysis Model is utilized to calculate ambient particle object
Concentration of component matrix, and the particle diameter distribution of the pollution sources class obtained in step 2 mark component is included in algorithm, obtain source spectrum
Matrix and source contribution matrix, identification pollution sources class and the contribution of quantitative contamination source class.
Atmospheric particulates are obtained using the acquisition of Atmospheric particulates source sampling device in step 1；Atmospheric particulates source sampling device is
Dilution tunnel sampler, ELPI sampler or settling flux sampler.
Pollution sources class includes city raised dust, coal-fired source, motor vehicle source, sulfate source, sources of nitrate, biomass in step 1
Combustion Source, soil source, building building cement source.
Component includes element class in step 1: Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Fe, Ni, Cu, Cd, Zn, Hg,
Pb etc.；Ionic species: Cl-、NO3 -、SO4 2-、NH4 +；EC/OC class: EC, OC；Organic matter species: organic matter hopance polycyclic aromatic hydrocarbon, left-handed
Glucan.
Element class is measured using ICP-MS, ICP-AES or XRF, and ionic species are measured using ion chromatography, and EC/OC class utilizes
The measurement of carbon component analysis instrument, organic species are measured with GC-MS.
The particle diameter distribution that component is identified in step 2 is mark component in PM2.5Middle content and mark component are in PM10Middle content
Ratio.
" multi-source spectrum matrix nesting parallel factor " analyzing three-dimensional Factor Analysis Model is to more partial size three dimensional receptors in step 4
Data carry out factorial analysis and obtain identical sources tribute to the particulate matter receptor of each partial size during Synchronization Analysis three-dimensional data
Trend matrix is offered, each partial size extracts respective derived components spectrum matrix respectively.
Multi-source spectrum matrix nesting parallel factor analysis " three-dimensional Factor Analysis Model:
In formula, xijkIt is the concentration of j-th of component in i-th of sample of kth different particle types object, unit μ g/m3；
aihIt is the contribution for simulating h-th of source class being calculated to i-th of sample, unit μ g/m3；
bjhkIt is the content of j-th of component in p-th of source constituents spectrum of kth different particle types object, unit g/g；
ejhkIt is the residual error of j-th of component in i-th of sample of kth different particle types object, unit μ g/m3。
Object concentration of component matrix is as follows up to formula in step 3: X1=c1AB1
X2=c2AB2
…
Xz=czABz
In formula, XzIt is the acceptor density data matrix of z-th of particles things；
A is the identical source contribution value matrix of each partial size；
BzIt is the derived components spectrum matrix of z-th of particles things；
czIt is scale factor.
SDABB source resolution system, it is characterised in that comprise the following modules:
Particulate count obtains acquisition and analyzes Atmospheric Particulate Matter source according to acquisition and analysis module for acquiring and analyzing
Class obtains the module of the component data in the class of Atmospheric Particulate Matter source；
Component particle diameter distribution module is identified, for the component data according to the Atmospheric Particulate Matter source class of acquisition, identification
Mark component in the class of Atmospheric Particulate Matter source and calculate mark component particle diameter distribution module；
Ambient particle object data collection and analysis module, for acquiring and analyzing the component in ambient particle object pollution sources class
Module；
The identification of source class and source class contribute quantitative module, utilize " multi-source spectrum matrix nesting parallel factor analysis " three-dimensional Factor minute
It analyses model and calculates ambient particle object concentration of component matrix, and the pollution sources class obtained in step 2 is identified to the particle diameter distribution of component
It is included in algorithm, obtains the module of source spectrum matrix and source contribution matrix, identification pollution sources class and the contribution of quantitative contamination source class.
SDABB source resolution terminal device, including memory, processor and storage are in memory and can be on a processor
The step of computer program of operation, processor realizes SDABB Source Apportionment when executing computer program.
The computer readable storage medium of SDABB line source analysis program is stored, computer-readable recording medium storage has meter
Calculation machine program is realized when computer program is executed by processor such as the step of the solution method of the source claim SDABB.
Embodiment 1
The technical solution of invention: SDABB Source Apportionment, comprising the following steps:
Step 1: acquiring and analyzes and obtains acquisition and analyze Atmospheric Particulate Matter source class, obtains Atmospheric Particulate Matter source
Component data in class；
Step 2: according to the component data of the Atmospheric Particulate Matter source class of acquisition, Atmospheric Particulate Matter source class is identified
In mark component and calculate mark component particle diameter distribution；
Step 3: acquiring and analyzes the component in ambient particle object pollution sources class；
Step 4: " multi-source spectrum matrix nesting parallel factor analysis " three-dimensional Factor Analysis Model is utilized to calculate ambient particle object
Concentration of component matrix, and the particle diameter distribution of the pollution sources class obtained in step 2 mark component is included in algorithm, obtain source spectrum
Matrix and source contribution matrix, identification pollution sources class and the contribution of quantitative contamination source class.
SDABB source resolution system, it is characterised in that comprise the following modules:
Particulate count obtains acquisition and analyzes Atmospheric Particulate Matter source according to acquisition and analysis module for acquiring and analyzing
Class obtains the module of the component data in the class of Atmospheric Particulate Matter source；
Component particle diameter distribution module is identified, for the component data according to the Atmospheric Particulate Matter source class of acquisition, identification
Mark component in the class of Atmospheric Particulate Matter source and calculate mark component particle diameter distribution module；
Ambient particle object data collection and analysis module, for acquiring and analyzing the component in ambient particle object pollution sources class
Module；
The identification of source class and source class contribute quantitative module, utilize " multi-source spectrum matrix nesting parallel factor analysis " three-dimensional Factor minute
It analyses model and calculates ambient particle object concentration of component matrix, and the pollution sources class obtained in step 2 is identified to the particle diameter distribution of component
It is included in algorithm, obtains the module of source spectrum matrix and source contribution matrix, identification pollution sources class and the contribution of quantitative contamination source class.
SDABB source resolution terminal device, including memory, processor and storage are in memory and can be on a processor
The computer program of operation, processor realize the step such as claim 1-6SDABB Source Apportionment when executing computer program
Suddenly.
The computer readable storage medium of SDABB line source analysis program is stored, computer-readable recording medium storage has meter
Calculation machine program is realized when computer program is executed by processor such as the step of the solution method of the source claim 1-6SDABB.
Embodiment 2
SDABB Source Apportionment, comprising the following steps:
Step 1: acquiring and analyzes and obtains acquisition and analyze Atmospheric Particulate Matter source class, obtains Atmospheric Particulate Matter source
Component data in class；
Wherein Atmospheric particulates are obtained using the acquisition of Atmospheric particulates source sampling device in step 1, Atmospheric particulates source sampling
Device is dilution tunnel sampler, ELPI sampler or settling flux sampler.
Wherein pollution sources class includes city raised dust, coal-fired source, motor vehicle source, sulfate source, sources of nitrate, life in step 1
Substance combustion source, soil source, building building cement source.
Step 2: according to the component data of the Atmospheric Particulate Matter source class of acquisition, Atmospheric Particulate Matter source class is identified
In mark component and calculate mark component particle diameter distribution；
Component includes element class in step 2: Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Fe, Ni, Cu, Cd, Zn, Hg,
Pb etc.；Ionic species: Cl-、NO3 -、SO4 2-、NH4 +；EC/OC class: EC, OC；Organic matter species: organic matter hopance polycyclic aromatic hydrocarbon, left-handed
Glucan.
Element class is measured using ICP-MS, ICP-AES or XRF, and ionic species are measured using ion chromatography, and EC/OC class utilizes
The measurement of carbon component analysis instrument, organic species are measured with GC-MS.
The method of identification source class mark component: consulting literatures can directly fix that provide every class source normal according to existing research
Mark component；The higher component of content and each biggish component of source class difference, as Al, Si in Fig. 3, Ca, Fe, Ti, OC,
EC、Cl-、SO4 2-；In summary, determine that the mark group of soil dirt is divided into Al, Si, Fe, Ti, Ti be according to existing research, and
A source class difference is big.
The mark group of coal-fired dirt is divided into Al, Si, OC, EC, Cl-、SO4 2-, the mark group of building building cement dirt is divided into Si, Ca.Step
The particle diameter distribution that component is identified in rapid two is mark component in PM2.5Middle content and mark component are in PM10The ratio of middle content.
Step 3: acquiring and analyzes the component in ambient particle object pollution sources class；
Step 4: " multi-source spectrum matrix nesting parallel factor analysis " three-dimensional Factor Analysis Model is utilized to calculate ambient particle object
Concentration of component matrix, and the particle diameter distribution of the pollution sources class obtained in step 2 mark component is included in algorithm, obtain source spectrum
Matrix and source contribution matrix, identification pollution sources class and the contribution of quantitative contamination source class.
" multi-source spectrum matrix nesting parallel factor " analyzing three-dimensional Factor Analysis Model is to more partial size three dimensional receptors in step 4
Data carry out factorial analysis and obtain identical sources tribute to the particulate matter receptor of each partial size during Synchronization Analysis three-dimensional data
Trend matrix is offered, each partial size extracts respective derived components spectrum matrix respectively.
Multi-source spectrum matrix nesting parallel factor analysis " three-dimensional Factor Analysis Model:
In formula, xijkIt is the concentration of j-th of component in i-th of sample of kth different particle types object, unit μ g/m3；
aihIt is the contribution for simulating h-th of source class being calculated to i-th of sample, unit μ g/m3；
bjhkIt is the content of j-th of component in p-th of source constituents spectrum of kth different particle types object, unit g/g；
ejhkIt is the residual error of j-th of component in i-th of sample of kth different particle types object, unit μ g/m3。
Object concentration of component matrix is as follows up to formula in step 3: X1=c1AB1
X2=c2AB2
…
Xz=czABz
In formula, XzIt is the acceptor density data matrix of z-th of particles things；
A is the identical source contribution value matrix of each partial size；
BzIt is the derived components spectrum matrix of z-th of particles things；
czIt is scale factor.
SDABB source resolution system, it is characterised in that comprise the following modules:
Particulate count is according to acquisition and analysis module, for acquiring and analyzing acquisition and analyze Atmospheric Particulate Matter source class,
Obtain the module of the component data in the class of Atmospheric Particulate Matter source；
Component particle diameter distribution module is identified, for the component data according to the Atmospheric Particulate Matter source class of acquisition, identification
Mark component in the class of Atmospheric Particulate Matter source and calculate mark component particle diameter distribution module；
Ambient particle object data collection and analysis module, for acquiring and analyzing the component in ambient particle object pollution sources class
Module；The identification of source class and source class contribute quantitative module, utilize " multi-source spectrum matrix nesting parallel factor analysis " three-dimensional Factor minute
It analyses model and calculates ambient particle object concentration of component matrix, and the pollution sources class obtained in step 2 is identified to the particle diameter distribution of component
It is included in algorithm, obtains the module of source spectrum matrix and source contribution matrix, identification pollution sources class and the contribution of quantitative contamination source class.
SDABB source resolution terminal device, including memory, processor and storage are in memory and can be on a processor
The step of computer program of operation, processor realizes SDABB Source Apportionment when executing computer program.
The computer readable storage medium of SDABB line source analysis program is stored, computer-readable recording medium storage has meter
Calculation machine program realizes the step of source SDABB solves method when computer program is executed by processor.
Compared with prior art, the application will identify component particle diameter distribution and be included in three-dimensional Factor minute in airborne particulate material resource class
Source Apportionment is analysed, contribution of the source class to particulate matter is calculated, obtained result more meets the physics of particulate matter source particle diameter distribution
Meaning improves the accuracy of source resolution result.
One embodiment of the present invention has been described in detail above, but content is only presently preferred embodiments of the present invention,
It should not be considered as limiting the scope of the invention.Any changes and modifications in accordance with the scope of the present application,
It should still be within the scope of the patent of the present invention.
Claims (10)
1.SDABB Source Apportionment, it is characterised in that the following steps are included:
Step 1: acquiring and analyzes Atmospheric Particulate Matter source class, obtains the component data in the class of Atmospheric Particulate Matter source；
Step 2: it according to the component data of the Atmospheric Particulate Matter source class of acquisition, identifies in the class of Atmospheric Particulate Matter source
Mark component and the particle diameter distribution for calculating mark component；
Step 3: acquiring and analyzes the component in ambient particle object pollution sources class；
Step 4: " multi-source spectrum matrix nesting parallel factor analysis " three-dimensional Factor Analysis Model is utilized to calculate ambient particle object component
Concentration matrix, and the particle diameter distribution of the pollution sources class obtained in step 2 mark component is included in algorithm, obtain source spectrum matrix
With source contribution matrix, pollution sources class and the contribution of quantitative contamination source class are identified.
2. SDABB Source Apportionment according to claim 1, it is characterised in that Atmospheric particulates utilize in the step 1
The acquisition of Atmospheric particulates source sampling device obtains；
Preferably, Atmospheric particulates source sampling device is dilution tunnel sampler, ELPI sampler or settling flux sampler.
3. SDABB Source Apportionment according to claim 1 or 2, it is characterised in that pollution sources class packet in the step 1
Include city raised dust, coal-fired source, motor vehicle source, sulfate source, sources of nitrate, biomass combustion source, soil source, building building cement source.
4. SDABB Source Apportionment according to claim 1 or 2, it is characterised in that component includes member in the step 1
Plain class: Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Fe, Ni, Cu, Cd, Zn, Hg, Pb etc.；Ionic species: Cl-、NO3 -、SO4 2-、
NH4 +；EC/OC class: EC, OC；Organic matter species: organic matter hopance polycyclic aromatic hydrocarbon, levoglucosan.
Preferably, element class is measured using ICP-MS, ICP-AES or XRF, and ionic species are measured using ion chromatography, EC/OC class benefit
It is measured with carbon component analysis instrument, organic species are measured with GC-MS.
5. SDABB Source Apportionment according to claim 1 or 2, it is characterised in that identify component in the step 2
Particle diameter distribution is mark component in PM2.5Middle content and the mark component are in PM10The ratio of middle content.
6. SDABB Source Apportionment according to claim 1 or 2, it is characterised in that " multi-source spectrum matrix in the step 4
Nested parallel factor " analyzing three-dimensional Factor Analysis Model carries out factorial analysis to more partial size three dimensional receptor data, in Synchronization Analysis
During three-dimensional data, identical source contribution trend matrix is obtained to the particulate matter receptor of each partial size, each partial size extracts respectively respectively
From derived components spectrum matrix.
Preferably, the multi-source spectrum matrix nesting parallel factor analysis " three-dimensional Factor Analysis Model:
In formula, xijkIt is the concentration of j-th of component in i-th of sample of kth different particle types object, unit μ g/m3；
aihIt is the contribution for simulating h-th of source class being calculated to i-th of sample, unit μ g/m3；
bjhkIt is the content of j-th of component in p-th of source constituents spectrum of kth different particle types object, unit g/g；
ejhkIt is the residual error of j-th of component in i-th of sample of kth different particle types object, unit μ g/m3。
7. SDABB Source Apportionment according to claim 1 or 2, it is characterised in that object component is dense in the step 4
It is as follows up to formula to spend matrix: X1=c1AB1
X2=c2AB2
…
Xz=czABz
In formula, XzIt is the acceptor density data matrix of z-th of particles things；
A is the identical source contribution value matrix of each partial size；
BzIt is the derived components spectrum matrix of z-th of particles things；
czIt is scale factor.
8.SDABB source resolution system, it is characterised in that comprise the following modules:
Particulate count obtains Atmospheric particulates for acquiring and analyzing Atmospheric Particulate Matter source class according to acquisition and analysis module
The module of component data in pollution sources class；
Component particle diameter distribution module is identified, for the component data according to the Atmospheric Particulate Matter source class of acquisition, identifies atmosphere
Mark component in the class of Particulate Pollution source and calculate mark component particle diameter distribution module；
Ambient particle object data collection and analysis module, for acquiring and analyzing the mould of the component in ambient particle object pollution sources class
Block；
The identification of source class and source class contribute quantitative module, utilize " multi-source spectrum matrix nesting parallel factor analysis " three-dimensional factorial analysis mould
Type calculates ambient particle object concentration of component matrix, and the particle diameter distribution of the pollution sources class obtained in step 2 mark component is included in
In algorithm, the module of source spectrum matrix and source contribution matrix, identification pollution sources class and the contribution of quantitative contamination source class is obtained.
9.SDABB source resolution terminal device, including memory, processor and storage are in the memory and can be at the place
The computer program run on reason device, which is characterized in that realized when the processor executes the computer program as right is wanted
The step of seeking SDABB Source Apportionment described in 1-7.
10. storing the computer readable storage medium of SDABB line source analysis program, the computer-readable recording medium storage has
Computer program, which is characterized in that the SDABB as described in claim 1-7 is realized when the computer program is executed by processor
Source solves the step of method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810778540.0A CN108956881A (en) | 2018-07-16 | 2018-07-16 | SDABB Source Apportionment, system, terminal device and computer readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810778540.0A CN108956881A (en) | 2018-07-16 | 2018-07-16 | SDABB Source Apportionment, system, terminal device and computer readable storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN108956881A true CN108956881A (en) | 2018-12-07 |
Family
ID=64481340
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810778540.0A Pending CN108956881A (en) | 2018-07-16 | 2018-07-16 | SDABB Source Apportionment, system, terminal device and computer readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108956881A (en) |
Cited By (10)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110779843A (en) * | 2019-09-30 | 2020-02-11 | 成都市环境保护科学研究院 | Typing analytic system of atmosphere secondary particle pollution process |
| CN111739587A (en) * | 2020-06-19 | 2020-10-02 | 中科三清科技有限公司 | Processing method and device of particulate matter monitoring data, storage medium and terminal |
| CN111739588A (en) * | 2020-06-19 | 2020-10-02 | 中科三清科技有限公司 | Method and device for analyzing atmospheric pollutant source, storage medium and terminal |
| CN111797918A (en) * | 2020-06-29 | 2020-10-20 | 平安国际智慧城市科技股份有限公司 | Atmospheric pollution source identification method and device, computer equipment and storage medium |
| CN111798932A (en) * | 2020-07-03 | 2020-10-20 | 中科三清科技有限公司 | Pollutant monitoring data processing method and device, storage medium and terminal |
| CN112304814A (en) * | 2020-10-16 | 2021-02-02 | 南开大学 | High-particle-size-resolution atmospheric particulate risk source analysis method |
| CN112562794A (en) * | 2020-12-09 | 2021-03-26 | 河北先进环保产业创新中心有限公司 | Pollution source analysis method and device and terminal equipment |
| CN112562796A (en) * | 2020-12-07 | 2021-03-26 | 中科三清科技有限公司 | Method and device for analyzing atmospheric pollution source, storage medium and terminal |
| CN112613675A (en) * | 2020-12-29 | 2021-04-06 | 南开大学 | Analyzing pollution source and meteorological factor to PM of different degrees2.5Machine learning model of pollution impact contributions and effects |
| CN112858116A (en) * | 2021-01-15 | 2021-05-28 | 南开大学 | Refined particulate matter source analysis method |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103698255A (en) * | 2013-12-26 | 2014-04-02 | 中国环境科学研究院 | Real-time analytic method for atmospheric particle source |
| CN107239613A (en) * | 2017-06-05 | 2017-10-10 | 南开大学 | A kind of intelligent source class recognition methods based on online data and Factor Analysis Model |
| CN107402284A (en) * | 2017-07-13 | 2017-11-28 | 南开大学 | The polynary source resolution algorithm of atmosphere pollution based on Gas identification |
| CN108052486A (en) * | 2017-12-01 | 2018-05-18 | 南开大学 | A kind of Source Apportionment that becomes more meticulous based on particulate matter inorganic component and organic marker |
-
2018
- 2018-07-16 CN CN201810778540.0A patent/CN108956881A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103698255A (en) * | 2013-12-26 | 2014-04-02 | 中国环境科学研究院 | Real-time analytic method for atmospheric particle source |
| CN107239613A (en) * | 2017-06-05 | 2017-10-10 | 南开大学 | A kind of intelligent source class recognition methods based on online data and Factor Analysis Model |
| CN107402284A (en) * | 2017-07-13 | 2017-11-28 | 南开大学 | The polynary source resolution algorithm of atmosphere pollution based on Gas identification |
| CN108052486A (en) * | 2017-12-01 | 2018-05-18 | 南开大学 | A kind of Source Apportionment that becomes more meticulous based on particulate matter inorganic component and organic marker |
Non-Patent Citations (3)
| Title |
|---|
| 刘立: "东莞/武汉城市大气颗粒物的理化特性与来源解析", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
| 吴虹等: "青岛环境空气PM10和PM2.5污染特征与来源比较", 《环境科学研究》 * |
| 田瑛泽等: "基于大气颗粒物物理化学特征的三维源解析技术研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 * |
Cited By (17)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110779843A (en) * | 2019-09-30 | 2020-02-11 | 成都市环境保护科学研究院 | Typing analytic system of atmosphere secondary particle pollution process |
| CN110779843B (en) * | 2019-09-30 | 2022-06-07 | 成都市环境保护科学研究院 | Typing analytic system of atmosphere secondary particle pollution process |
| CN111739587A (en) * | 2020-06-19 | 2020-10-02 | 中科三清科技有限公司 | Processing method and device of particulate matter monitoring data, storage medium and terminal |
| CN111739588A (en) * | 2020-06-19 | 2020-10-02 | 中科三清科技有限公司 | Method and device for analyzing atmospheric pollutant source, storage medium and terminal |
| CN111739588B (en) * | 2020-06-19 | 2021-01-01 | 中科三清科技有限公司 | Method and device for analyzing atmospheric pollutant source, storage medium and terminal |
| CN111739587B (en) * | 2020-06-19 | 2021-02-02 | 中科三清科技有限公司 | Processing method and device of particulate matter monitoring data, storage medium and terminal |
| CN111797918A (en) * | 2020-06-29 | 2020-10-20 | 平安国际智慧城市科技股份有限公司 | Atmospheric pollution source identification method and device, computer equipment and storage medium |
| CN111797918B (en) * | 2020-06-29 | 2024-04-16 | 平安国际智慧城市科技股份有限公司 | Atmospheric pollution source identification method, apparatus, computer device and storage medium |
| CN111798932B (en) * | 2020-07-03 | 2021-06-15 | 中科三清科技有限公司 | Pollutant monitoring data processing method and device, storage medium and terminal |
| CN111798932A (en) * | 2020-07-03 | 2020-10-20 | 中科三清科技有限公司 | Pollutant monitoring data processing method and device, storage medium and terminal |
| CN112304814A (en) * | 2020-10-16 | 2021-02-02 | 南开大学 | High-particle-size-resolution atmospheric particulate risk source analysis method |
| CN112562796A (en) * | 2020-12-07 | 2021-03-26 | 中科三清科技有限公司 | Method and device for analyzing atmospheric pollution source, storage medium and terminal |
| CN112562794A (en) * | 2020-12-09 | 2021-03-26 | 河北先进环保产业创新中心有限公司 | Pollution source analysis method and device and terminal equipment |
| CN112562794B (en) * | 2020-12-09 | 2022-11-08 | 河北先进环保产业创新中心有限公司 | Pollution source analysis method and device and terminal equipment |
| CN112613675A (en) * | 2020-12-29 | 2021-04-06 | 南开大学 | Analyzing pollution source and meteorological factor to PM of different degrees2.5Machine learning model of pollution impact contributions and effects |
| CN112858116A (en) * | 2021-01-15 | 2021-05-28 | 南开大学 | Refined particulate matter source analysis method |
| CN112858116B (en) * | 2021-01-15 | 2023-10-17 | 南开大学 | Fine particle source analysis method |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108956881A (en) | SDABB Source Apportionment, system, terminal device and computer readable storage medium | |
| CN114424058B (en) | Tracing method for VOCs pollution | |
| Harrison et al. | A pragmatic mass closure model for airborne particulate matter at urban background and roadside sites | |
| Hopke | Review of receptor modeling methods for source apportionment | |
| Hudson et al. | Phosphate concentrations in lakes | |
| Lucarelli et al. | Is PIXE still a useful technique for the analysis of atmospheric aerosols? The LABEC experience | |
| CN105631203A (en) | Method for recognizing heavy metal pollution source in soil | |
| Wang et al. | Heavy metals in urban ambient PM 10 and soil background in eight cities around China | |
| CN108052486B (en) | Fine source analysis method based on inorganic components and organic markers of particulate matters | |
| Park et al. | Field evaluation of a near–real time elemental monitor and identification of element sources observed at an air monitoring supersite in Korea | |
| CN108446531B (en) | Receptor Si and Al reconstruction method based on actual measurement source component spectrum and source analysis technology | |
| Ching et al. | Aerosol mixing state revealed by transmission electron microscopy pertaining to cloud formation and human airway deposition | |
| CN102393401A (en) | Detection method of heavy metal element content in air pollutant particles | |
| CN109187287A (en) | Source of atmospheric particulate matter analytic method based on stable element particle diameter distribution information | |
| Samek et al. | Preliminary PM2. 5 and PM10 fractions source apportionment complemented by statistical accuracy determination | |
| Kajino et al. | Modeling transition metals in East Asia and Japan and its emission sources | |
| Orlović-Leko et al. | Physico-chemical characterization of an urban rainwater (Zagreb, Croatia) | |
| Varela et al. | On the use of epigaeic mosses to biomonitor atmospheric deposition of nitrogen | |
| Su et al. | Development and application of a mass closure PM 2.5 composition online monitoring system | |
| Wen et al. | Source apportionment of PM 2.5 in Tangshan, China—Hybrid approaches for primary and secondary species apportionment | |
| Cohen et al. | Ion beam techniques for source fingerprinting fine particle air pollution in major Asian-Pacific cities | |
| CN108805456A (en) | A kind of assessment method of the air major pollutants based on principal component analysis | |
| Gnauk et al. | Aerosol characterisation at the FEBUKO upwind station Goldlauter (I): Particle mass, main ionic components, OCEC, and mass closure | |
| Johnson et al. | Chemical and physical analyses of Houston aerosol for interlaboratory comparison of source apportionment procedures | |
| Tahri et al. | TXRF characterization and source identification by positive matrix factorization of airborne particulate matter sampled in Kenitra City (Morocco) |
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: 20181207 |