李　源,楊　寧，畢溫凱 ，沈　勁，高璟赟

(1. 天津市環(huán)境監(jiān)測(cè)中心，天津 300191；2.廣東省環(huán)境監(jiān)測(cè)中心, 國(guó)家環(huán)境保護(hù)區(qū)域空氣質(zhì)量監(jiān)測(cè)重點(diǎn)實(shí)驗(yàn)室，廣東 廣州 510045)

三維空氣質(zhì)量模型污染來源追因技術(shù)概述

李源1,楊寧1，畢溫凱1，沈勁2，高璟赟1

(1. 天津市環(huán)境監(jiān)測(cè)中心，天津 300191；2.廣東省環(huán)境監(jiān)測(cè)中心, 國(guó)家環(huán)境保護(hù)區(qū)域空氣質(zhì)量監(jiān)測(cè)重點(diǎn)實(shí)驗(yàn)室，廣東 廣州 510045)

摘要：為應(yīng)對(duì)區(qū)域大氣復(fù)合污染的嚴(yán)峻態(tài)勢(shì)，亟需建立健全空氣污染來源追因系統(tǒng)。三維空氣質(zhì)量模型可以克服外場(chǎng)觀測(cè)與煙霧箱模擬高成本與低時(shí)效的缺點(diǎn)。主要介紹了使用空氣質(zhì)量模型解決大氣污染物來源分析的方法，其中：過程分析(PA)可用于估算各物理過程與凈化學(xué)過程對(duì)污染物濃度的影響；臭氧來源解析技術(shù)(OSAT)和顆粒物來源解析技術(shù)(PSAT)通過標(biāo)記不同地區(qū)各類排放源一次污染物的排放，得出不同受體點(diǎn)位主要一次與二次污染物的來源或主要排放源區(qū)對(duì)區(qū)域各項(xiàng)污染物濃度的貢獻(xiàn)。

關(guān)鍵詞：大氣復(fù)合污染；三維空氣質(zhì)量模型；污染來源追因；技術(shù)概述

0引言

隨著經(jīng)濟(jì)的持續(xù)快速發(fā)展，我國(guó)不少區(qū)域已出現(xiàn)較為嚴(yán)重的以PM2.5和臭氧為主要污染物的大氣復(fù)合污染[1，2]，污染來源追因?qū)τ谙嚓P(guān)部門有針對(duì)性地控制污染有重要意義。

外場(chǎng)觀測(cè)可以直觀準(zhǔn)確地了解大氣污染情況[3]，煙霧箱實(shí)驗(yàn)對(duì)大氣化學(xué)的機(jī)理研究也有重要意義[4]，如在洛杉磯的觀測(cè)與煙霧箱研究中發(fā)現(xiàn)了光化學(xué)煙霧污染，臭氧是主要污染物，NOx和VOCs是臭氧重要的前體物[5]。然而，外場(chǎng)觀測(cè)與煙霧箱等研究手段都依賴于大量的監(jiān)測(cè)設(shè)備與分析實(shí)驗(yàn)，不但耗費(fèi)大量人力物力，而且比較耗時(shí)。外場(chǎng)觀測(cè)中，采樣點(diǎn)的代表性會(huì)影響對(duì)區(qū)域污染的理解，采樣的空間與時(shí)間尺度的連續(xù)性也較差。使用數(shù)值模擬有助于全面了解區(qū)域污染概況，對(duì)污染的產(chǎn)生、發(fā)展過程及影響因素都可以進(jìn)行深入研究，節(jié)約人力、物力及時(shí)間，是對(duì)目前有限的觀測(cè)與煙霧箱模擬研究的重要補(bǔ)充[6，7]。

利用數(shù)學(xué)的方法，綜合考慮各種過程和影響因素，定量描述污染物在大氣中的遷移、轉(zhuǎn)化規(guī)律的模式通常稱為空氣質(zhì)量模式(或模型)[8，9]。早期的模式以局地?zé)熈鲾U(kuò)散模式、盒子模式和拉格朗日軌跡模式為主[10]，基于擴(kuò)散的模式一般對(duì)化學(xué)反應(yīng)考慮得相對(duì)簡(jiǎn)單，不適用于二次污染物的研究。盒子模式雖然考慮了較復(fù)雜的化學(xué)反應(yīng)，但較少考慮污染物的空間差異及傳輸擴(kuò)散因素，僅適用于局地研究?？諝赓|(zhì)量數(shù)值模擬研究始于20世紀(jì)60年代，70年代末隨著各種大氣物理化學(xué)過程研究的深入，空氣質(zhì)量模型開始向精細(xì)化方向發(fā)展，逐步發(fā)展了以歐拉網(wǎng)格模型為主的空氣質(zhì)量模型[11]。90年代起，美國(guó)環(huán)保局開始致力開發(fā)綜合的第三代空氣質(zhì)量模擬系統(tǒng)，并提出了“一個(gè)大氣”的概念[12]?，F(xiàn)在，把大氣污染與氣候變化偶合的模型也有了很大的發(fā)展與較多的應(yīng)用，空氣質(zhì)量模型除了可用于模擬污染物的濃度外，不少用于污染來源追因的工具模塊也有所發(fā)展，如綜合空氣質(zhì)量模型CAMx的擴(kuò)展模塊[13]：過程分析(PA)、臭氧來源解析技術(shù)(OSAT)和顆粒物來源解析技術(shù)(PSAT)等。

1過程分析簡(jiǎn)介

污染物濃度的變化受化學(xué)過程、水平與垂直輸送、沉降等一系列復(fù)雜大氣過程的綜合影響，過程分析(PA)或綜合過程分析(IPR)可用于估算各物理過程與凈化學(xué)過程對(duì)污染物模擬結(jié)果的影響。過程分析通過分別計(jì)算歐拉連續(xù)性方程等號(hào)右邊的各項(xiàng)，便可以得出氣相化學(xué)過程、水平平流與擴(kuò)散、垂直對(duì)流與擴(kuò)散、沉降和其它過程對(duì)指定網(wǎng)格的污染物濃度的小時(shí)貢獻(xiàn)量[14]。

(1)

以臭氧為例，等式右邊第一項(xiàng)表示水平平流，東南西北邊界的臭氧水平通量可以區(qū)別對(duì)待，即可以定量某網(wǎng)格或區(qū)域四個(gè)邊界對(duì)該網(wǎng)格或區(qū)域的臭氧貢獻(xiàn)；第二項(xiàng)算出垂直對(duì)流對(duì)臭氧的貢獻(xiàn)；第三項(xiàng)是擴(kuò)散對(duì)臭氧的貢獻(xiàn)；第四項(xiàng)是排放對(duì)臭氧的貢獻(xiàn)，由于臭氧的直接排放很少，源清單并沒有考慮臭氧的一次排放，因此這一項(xiàng)為0；第五項(xiàng)是化學(xué)過程對(duì)臭氧的貢獻(xiàn)，其可以分為氣相化學(xué)與非均相化學(xué)兩部分，但由于臭氧的生成主要受氣相化學(xué)的影響，因此本研究并沒有考慮非均相化學(xué)對(duì)臭氧的貢獻(xiàn)；最后一項(xiàng)是沉降過程對(duì)臭氧的貢獻(xiàn)，可分為干沉降與濕沉降兩部分。

基于過程的分析對(duì)深入研究污染有重要意義，有助于定量不同物理化學(xué)過程的影響，以明確各邊界的臭氧輸通量、本地化學(xué)生成的強(qiáng)烈程度等，能深化區(qū)域污染的認(rèn)識(shí)。目前過程分析技術(shù)已在我國(guó)有所應(yīng)用。研究結(jié)果表明，水平輸送過程在1、4、10月是影響臭氧的主要因子；而在7月，氣相化學(xué)與垂直輸送是影響臭氧濃度的主要過程。對(duì)于近地面顆粒物，一次排放與氣溶膠過程是顆粒物濃度增長(zhǎng)的主要因素，而水平輸送則是顆粒物去除的主要過程[15]。在珠三角地區(qū)使用CMAQ模型和過程分析技術(shù)，結(jié)果表明白天光化學(xué)生成過程對(duì)大部分地區(qū)邊界層臭氧濃度升高貢獻(xiàn)最大[16]。

2臭氧來源解析及生成追蹤技術(shù)

追蹤臭氧的來源既有利于了解污染過程，也有利于控制方案的制定。Morris 等發(fā)展了臭氧來源解析技術(shù)OSAT[17，18]，并應(yīng)用于美國(guó)東北部的臭氧研究。目前OSAT已有了新發(fā)展，這一技術(shù)主要對(duì)不同地區(qū)各類源排放的NOx與VOCs分別進(jìn)行標(biāo)記，如第i類地區(qū)第j類源排放的NOx與VOCs分別表示為Ni,j與Vi,j。以H2O2與HNO3的生成速率比值PH2O2/PHNO3=0.35為劃分臭氧生成的NOx與VOC控制區(qū)的標(biāo)準(zhǔn)[19]，在珠三角的研究表明0.35作為閾值同樣適用于我國(guó)的南方地區(qū)[20]。OSAT對(duì)PH2O2/PHNO3=0.35這一閾值并不敏感，其在0.25～0.5波動(dòng)時(shí)對(duì)臭氧源解析結(jié)果的影響只在4%～6%[21]。臭氧的生成與消耗分別為PO3與DO3(PO3=ΔO3-DO3)，在NOx控制區(qū)下生成的臭氧表示為O3Ni,j，在VOC控制區(qū)下生成的臭氧表示為O3Vi,j，則

(2)

(3)

其中最大增量反應(yīng)性MIR可使用模型內(nèi)設(shè)的默認(rèn)值。同時(shí)，OSAT考慮了若干臭氧的消耗過程，分別是O3+VOCs，HOx+O3的反應(yīng)以及兩個(gè)與氧原子去除相關(guān)的反應(yīng)O(3P)+VOCs和O(1D)+O3。在不同的控制區(qū)下，有

(4)

其中X為N或V。

GOAT與OSAT類似，但GOAT不直接示蹤臭氧的前體物，而是對(duì)生成的臭氧進(jìn)行示蹤，如在第i類地區(qū)生成的臭氧標(biāo)記為O3i，以確定影響受體點(diǎn)的臭氧生成地。

GOAT解釋出的邊界(BC)對(duì)臭氧的貢獻(xiàn)為模擬區(qū)域外的臭氧直接傳輸?shù)绞荏w點(diǎn)的臭氧貢獻(xiàn)，約等于背景臭氧濃度；而OSAT解釋出的BC對(duì)臭氧貢獻(xiàn)為模擬區(qū)域外的臭氧及因前體物輸送而生成的臭氧對(duì)受體點(diǎn)的臭氧貢獻(xiàn)，兩者相減可得邊界前體物輸送對(duì)受體點(diǎn)的臭氧貢獻(xiàn)。

在珠三角地區(qū)，已有較多臭氧來源解析的相關(guān)研究[22-24]，按行政區(qū)把珠三角劃分為55個(gè)區(qū)域，排放源分為8類，分別是天然源、生物質(zhì)燃燒、道路交通尾氣、溶劑使用與存儲(chǔ)、生活與廢物處理、大點(diǎn)源、非溶劑類工業(yè)和其它排放源。按排放類型分，道路交通尾氣對(duì)珠三角各受體點(diǎn)的臭氧貢獻(xiàn)最大；在珠三角，排放源一般對(duì)下風(fēng)向40km范圍內(nèi)的地區(qū)臭氧貢獻(xiàn)最大。

3顆粒物來源識(shí)別方法

顆粒物來源識(shí)別技術(shù)(PSAT)是CAMX的一個(gè)擴(kuò)展功能，用于開展目標(biāo)區(qū)域顆粒物及其化學(xué)組分濃度的排放源類別和源地區(qū)來源研究[25]。主要原理是對(duì)一次顆粒物、二次顆粒物及其氣態(tài)前體物加入反應(yīng)性示蹤物進(jìn)行標(biāo)識(shí)，追蹤其在物理過程和化學(xué)過程中的生成、消除和轉(zhuǎn)化，以此量化各區(qū)域、各排放源對(duì)受體顆粒物濃度的貢獻(xiàn)大小，并識(shí)別其重要的污染源。

PSAT沿用CAMx氣溶膠模塊對(duì)顆粒物物種的劃分，分別對(duì)各個(gè)模型物種進(jìn)行來源分析。目前，PSAT支持的PM物種共有六類：硫酸鹽(PSO4)、硝酸鹽顆粒物(PNO3)、銨鹽(PNH4)、含汞顆粒物(Hg(p))、二次有機(jī)氣溶膠(SOA)和一次排放顆粒物(PPM2.5)。其中，一次排放顆粒物包括元素碳(EC)、一次有機(jī)氣溶膠(POA)、地殼細(xì)顆粒物(FCRS)、其他細(xì)顆粒物(FPRM)、粗地殼粗顆粒物(CCRS)和其他粗顆粒物(CPRM)。PSAT技術(shù)的一個(gè)基本假設(shè)是對(duì)每類顆粒物追蹤其主要前體物，對(duì)不同物種，PSAT采用不同的標(biāo)識(shí)物種(tracers)。

PSAT需要把模擬區(qū)域和排放源清單數(shù)據(jù)分成多個(gè)清單輸入組(emissiongroups)， 并將其全部應(yīng)用于所要追蹤的標(biāo)識(shí)物種，并根據(jù)物種的排放、沉降、遷移和化學(xué)過程分別計(jì)算每個(gè)物種在每個(gè)模擬網(wǎng)格中的變化量[26]。對(duì)于簡(jiǎn)單的一次排放源的PM標(biāo)識(shí)物種，PSAT直接采用公式(5)計(jì)算各清單輸入組中該物種的貢獻(xiàn)率。對(duì)于涉及化學(xué)反應(yīng)的模型標(biāo)識(shí)物種B，(例如A→B)，PSAT采用公式(6)計(jì)算各清單輸入組中該物種的貢獻(xiàn)率。此外，對(duì)于部分化學(xué)過程，需要考慮標(biāo)識(shí)物種在輸入組中的權(quán)重大小，例如在以芳香烴為標(biāo)志物時(shí)，由于不同輸入組的甲苯和二甲苯量不同，芳香烴的化學(xué)衰變也存在差異(此時(shí)采用公式(7)和公式(8))。以上公式均可用于計(jì)算標(biāo)識(shí)物種在排放、沉降和遷移過程的變化，對(duì)于排放過程，標(biāo)識(shí)物種ai為污染源排放輸入量；對(duì)于遷移過程，ai為上風(fēng)向網(wǎng)格中標(biāo)識(shí)物種的濃度。

ai(t+Δt)=ai(t)+ΔA(ai/Σai)

(5)

bi(t+Δt)=bi(t)+ΔB(ai/Σai)

(6)

ai(t+Δt)=ai(t)+ΔA(wiai/Σwiai)

(7)

bi(t+Δt)=bi(t)+ΔB(wiai/Σwiai)

(8)

其中：A為模型標(biāo)識(shí)物種，ai為各排放源清單輸入組中的A模型標(biāo)識(shí)物種，wi為ai的權(quán)重因子。

PSAT的主要功能是基于示蹤物方法來追蹤目標(biāo)區(qū)域網(wǎng)格顆粒物濃度的源區(qū)域和源類別貢獻(xiàn)，為建立污染源與環(huán)境受體點(diǎn)污染濃度的響應(yīng)關(guān)系，需要設(shè)置污染源類別和源追蹤區(qū)域，以把模擬區(qū)域和排放源清單數(shù)據(jù)分成多個(gè)清單輸入組，并選取受體點(diǎn)。在上海市曾應(yīng)用過這一技術(shù)[27]，結(jié)果表明：人為排放源中，工業(yè)鍋爐和窯爐、移動(dòng)源和電站鍋爐是對(duì)細(xì)顆粒物中硝酸鹽貢獻(xiàn)最大的3類排放源；工業(yè)源和移動(dòng)源是對(duì)硫酸鹽貢獻(xiàn)最大的兩類排放源。在灰霾、濕霾和過境這3次污染過程中，上海本地排放對(duì)PM2.5的濃度貢獻(xiàn)分別是35.3%、44.8%和22.7%。嚴(yán)重污染過程并非單一城市所致，區(qū)域聯(lián)防聯(lián)控，特別是重度污染期間的聯(lián)合減排對(duì)于緩解細(xì)顆粒物重度污染極為重要。

4小結(jié)

有效應(yīng)對(duì)和妥善處置大氣重污染過程，亟需盡快建立和完善空氣污染來源追因技術(shù)體系。本文旨在介紹三維空氣質(zhì)量模型中用于解析污染來源的技術(shù)原理及其應(yīng)用情況。其中，過程分析(PA)可用于估算各物理過程與凈化學(xué)過程對(duì)污染物濃度的影響，臭氧來源解析技術(shù)(OSAT)和顆粒物來源解析技術(shù)(PSAT)通過標(biāo)記不同地區(qū)各類排放源一次污染物的排放，得出不同受體點(diǎn)位主要一次與二次污染物的來源或主要排放源區(qū)對(duì)區(qū)域各項(xiàng)污染物濃度的貢獻(xiàn)。使用數(shù)值模型與來源追因模塊相結(jié)合的方式是高效解決我國(guó)空氣污染聯(lián)防聯(lián)控與提高環(huán)境管理能力的主要出路。

參考文獻(xiàn)：

[1]ZhangYH,HuM,ZhongLJ,etal.RegionalIntegratedExperimentsonAirQualityoverPearlRiverDelta2004 (PRIDE-PRD2004):Overview[J].AtmosphericEnvironment, 2008(42): 6157-6173.

[2]ChanCK,YaoX.AirpollutioninmegacitiesinChina[J].AtmosphericEnvironment, 2008, 42(1):1-42.

[3]ZhangYH,SuH,ZhongLJ,etal.Regionalozonepollutionandobservation-basedapproachforanalyzingozone-precursorrelationshipduringthePRIDE-PRD2004campaign[J].AtmosphericEnvironment, 2008(42): 6203-6218.

[4]王煒罡,李坤,周力, 等. 用于研究大氣氧化過程和機(jī)制的雙反應(yīng)器煙霧箱的評(píng)估和應(yīng)用 [J]. 物理化學(xué)學(xué)報(bào),2015(7):1251-1259.

[5]Haagen-SmitAJ.ChemistryandPhysiologyofLos-AngelesSmog[J].IndEngChem, 1952( 44): 1342-1346.

[6]王自發(fā), 謝付瑩, 王喜全, 等. 嵌套網(wǎng)格空氣質(zhì)量預(yù)報(bào)模式系統(tǒng)的發(fā)展與應(yīng)用[J]. 大氣科學(xué), 2006, 30(5): 778-790.

[7]ZhangY,PunB,WuSY,etal.ApplicationandevaluationoftwoairqualitymodelsforparticulatematterforasoutheasternU.S.episode[J].JournaloftheAir&WasteManagementAssociation, 2004, 54(12): 1478-1493.

[8]唐孝炎, 張遠(yuǎn)航,邵敏. 大氣環(huán)境化學(xué):第二版[M].北京：高等教育出版社,2006.

[9]劉峰, 張遠(yuǎn)航, 蘇杭, 等. 大氣化學(xué)傳輸模式CAMx的伴隨模式:構(gòu)建及應(yīng)用[J]. 北京大學(xué)學(xué)報(bào):自然科學(xué)版, 2007(43): 764-770.

[10]Holmes N S,Morawska L. A review of dispersion modelling and its application to the dispersion of particles: An overview of different dispersion models available [J]. Atmospheric Environment, 2006(40): 5902-5928.[11]Seinfeld J H. Ozone air quality models. A critical review [J]. Journal of the Air Pollution Control Association, 1988(38): 616-645.[12]Byun D W, Ching J K S. Science algorithms of the EPA Models-3 community multiscale air quality (CMAQ) modeling system[R]. Office of Research and Development, United States Environmental Protection Agency, 1999, EPA/600/R-99/030: ES1-ES8.

[13]Nopmongcol U, Koo B, Tai E, et al. Modeling Europe with CAMx for the Air Quality Model Evaluation International Initiative (AQMEII)[J]. Atmospheric Environment, 2012, 53(2):177-185.

[14]Shen J, Wang X S, Li J F, et al. Evaluation and intercomparison of ozone simulations by Models-3/CMAQ and CAMx over the Pearl River Delta[J]. Science China Chemistry, 2011, 54(11):1750-1762.

[15]Liu X H, Zhang Y, Xing J, et al. Understanding of regional air pollution over China using CMAQ, part II. Process analysis and sensitivity of ozone and particulate matter to precursor emissions [J]. Atmospheric Environment, 2010, 44(30):3719-3727.

[16]Wang X, Zhang Y, Hu Y, et al. Process analysis and sensitivity study of regional ozone formation over the Pearl River Delta, China, during the PRIDE-PRD2004 campaign using the CMAQ model [J]. Atmospheric Chemistry and Physics, 2010(10): 4423-4437.

[17]Morris R, Tai E, Wilson G, et al. Use of advanced ozone source apportionment techniques to estimate area of influence (AOI) of emissions contributions to elevated ozone concentrations[R]. Proc a Waste Manage Assoc Annu Meet Exhib, 1998: 31-53.

[18]Morris R E, Tai E, Hersey C, et al. A Methodology for Quantifying Ozone Transport and Assessing the Benefits of Alternative Control Strategies[R]. The A&WMA Annual Meeting & Exhibition, Orlando, FL, 2001.

[19]Sillman, S. The use of NO y, H2O2, and HNO3as indicators for ozone‐NO x ‐hydrocarbon sensitivity in urban locations [J]. Journal of Geophysical Research, 1995, 100(D7): 14175-14188.

[20]沈勁,鐘流舉,陳皓,等. H2O2與HNO3生成速率比值判別臭氧生成敏感性[J]. 中國(guó)科技論文,2014(6):725-728.

[21]Dunker A M, Yarwood G, Ortmann J P, et al. Comparison of Source Apportionment and Source Sensitivity of Ozone in a Three-Dimensional Air Quality Model [J]. Environ. Sci. Technol, 2002, 36(13): 2953-2964.

[22]沈勁,陳皓,鐘流舉. 珠三角秋季臭氧污染來源解析[J]. 環(huán)境污染與防治,2015(1):25-30.

[23]陳皓,王雪松,沈勁,陸克定,張遠(yuǎn)航. 珠江三角洲秋季典型光化學(xué)污染過程中的臭氧來源分析[J]. 北京大學(xué)學(xué)報(bào):自然科學(xué)版,2015(4):620-630.

[24]Li Y, Lau A K H, Fung J C H, et al. Ozone source apportionment (OSAT) to differentiate local regional and super‐regional source contributions in the Pearl River Delta region, China[J]. Journal of Geophysical Research Atmospheres, 2012, 117(D15), doi:10.1029/2011JD017340.

[25]Wagstrom K M, Pandis S N, Yarwood G, et al. Development and application of a computationally efficient particulate matter apportionment algorithm in a three-dimensional chemical transport model[J]. Atmospheric Environment, 2008, 42(22):5650-5659.

[26]潘月云.城市細(xì)顆粒物本地排放源識(shí)別與區(qū)域輸送影響研究[D].廣州:華南理工大學(xué),2015.

[27]李莉,安靜宇,嚴(yán)茹莎. 基于細(xì)顆粒物來源追蹤技術(shù)的2013年12月上海市嚴(yán)重污染過程中PM2.5的源貢獻(xiàn)分析[J]. 環(huán)境科學(xué),2015(10):3543-3553.

Overview of Air Pollutants Apportionment Technology by 3-D Air Quality Model

LI Yuan1, YANG Ning1,BI Wen-Kai1,SHEN Jin2, GAO Jing-Yun2

(1. Tianjin Environmental Monitoring Center，Tianjin 300191, China)

Abstract：To cope with the severe situation of the regional air pollution, it is urgent to establish a sound source apportionment system of air pollution. 3-D air quality model could overcome the shortcomings of high cost and inefficiency in field observation or smog chamber. Methods of using air quality model to solve the problem of air pollutants source apportionment were mainly introduced, in order to provide reference about establishing an environment management system in line with the international advanced level. Process analysis (PA) could be used to estimate the influence of the physical and chemical processes on the pollutants concentration.Through marking emissions of pollutants of different types and in different areas, ozone source apportionment technology (OSAT) and particulate source apportionment technology (PSAT) could figure out the sources of mainly primary and secondary pollutants in different receptors or the regional contribution of various pollutants by main emission areas.

Key words：air complex pollution; 3-D air quality model; pollution source apportionment; summary of technology

收稿日期：2015-12-08

基金項(xiàng)目：天津市重大科技專項(xiàng)項(xiàng)目(14ZCDGSF00027);國(guó)家科技支撐計(jì)劃項(xiàng)目(2014BAC23B01)。

作者簡(jiǎn)介：李源(1988-)，碩士，主要研究空氣質(zhì)量監(jiān)測(cè)與模擬。

通信作者：沈勁，工程師，主要研究空氣污染與監(jiān)測(cè)。

中圖分類號(hào)：X701

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1673-9655(2016)04-0072-04