• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    基于觀測的模型(OBM)的發(fā)展歷程及其在我國大氣化學研究中的應用與展望

    2022-03-24 08:36:12張英南薛麗坤陳天舒申恒青王文興
    環(huán)境科學研究 2022年3期
    關鍵詞:前體大氣觀測

    張英南,薛麗坤*,陳天舒,申恒青,李 紅,王文興,

    1. 山東大學環(huán)境研究院,山東 青島 266237

    2. 中國環(huán)境科學研究院,環(huán)境基準與風險評估國家重點實驗室,北京 100012

    自2013年以來,國務院先后頒布了《大氣污染防治行動計劃》《打贏藍天保衛(wèi)戰(zhàn)三年行動計劃》等一系列行動計劃,我國大氣污染治理工作不斷深入推進. 經(jīng)過長期努力,大氣監(jiān)測常規(guī)六參數(shù)中如SO2、NOx、CO、PM10、PM2.5的環(huán)境濃度已經(jīng)顯著降低[1].但近年來,我國臭氧(O3)污染呈惡化態(tài)勢[2-6],同時許多城市PM2.5濃度年均值尚遠未達標[7],且PM2.5中二次組分質(zhì)量濃度的下降是未來PM2.5控制的重點. 以PM2.5和O3為首要污染物的區(qū)域大氣復合污染問題已成為制約我國空氣質(zhì)量持續(xù)改善的關鍵問題和重要挑戰(zhàn)[8].

    準確確定二次污染物與前體物的非線性關系是有效控制二次污染的前提. OBM(observation-based model,基于觀測的模型)內(nèi)置大氣化學機制,運用數(shù)學方法和計算機語言并結合觀測數(shù)據(jù),能夠基于實際大氣環(huán)境模擬解析大氣化學反應的速率,是診斷二次污染生成機制的重要手段和工具之一[9]. 與基于排放的空氣質(zhì)量模型相比,OBM避免了排放清單不確定性的影響,但也對觀測數(shù)據(jù)質(zhì)量提出了較高要求,一般需要輸入詳細VOCs (揮發(fā)性有機物)組分的觀測濃度及其他相關參數(shù)[10]. 當前階段,我國嚴峻的大氣復合污染防控形勢和日臻完善的空氣質(zhì)量監(jiān)測網(wǎng)絡(特別是大氣污染監(jiān)測超級站和光化學監(jiān)測網(wǎng)的建設)為OBM在國內(nèi)大氣化學研究中的進一步應用提供了強烈的現(xiàn)實需求和有力的數(shù)據(jù)支撐. 同時,大氣化學機制和外部運行平臺的持續(xù)更新也進一步推動了OBM在大氣化學領域中的應用,使其逐步發(fā)展成為一種功能強大的大氣污染成因分析與化學機制解析工具,被用于解析國內(nèi)O3和二次氣溶膠的生成機制以及其他活性組分化學機制. 但是,OBM的應用也存在一些問題:一方面,我國可利用的大氣顆粒物組分和詳細VOCs組分的觀測數(shù)據(jù)仍有限且主要集中于經(jīng)濟發(fā)達的東部地區(qū),直接限制了OBM在大氣化學研究中(特別是二次有機氣溶膠生成機制解析)的進一步應用;另一方面,我國缺乏具有自主知識產(chǎn)權的大氣化學機制和OBM運行平臺,仍需加強對關鍵核心技術的攻堅.

    鑒于此,筆者在文獻調(diào)研的基礎上梳理了OBM的發(fā)展歷程,主要涉及OBM的產(chǎn)生、結構特征和內(nèi)置大氣化學機制的發(fā)展歷程;總結了國內(nèi)研究中基于OBM解析大氣化學過程的相關研究成果,重點關注O3和二次氣溶膠生成機制及其他活性成分化學機制的研究;最后,總結了OBM應用現(xiàn)狀下存在的問題并展望OBM未來的發(fā)展方向,以期進一步推進OBM的完善與應用,為我國區(qū)域大氣復合污染防治工作做出力所能及的貢獻.

    1 OBM的發(fā)展歷程

    1.1 OBM的產(chǎn)生

    20世紀50年代,以高濃度O3為標志的光化學煙霧籠罩洛杉磯地區(qū),揭開了美國數(shù)十年與光化學污染抗爭之路的帷幕. 20世紀90年代,根據(jù)《清潔空氣法》修改條例中的規(guī)定,美國環(huán)境保護局(US EPA)在24個城市地區(qū)建立了以O3、NOx和VOCs為監(jiān)測目 標的PAMS (photochemical assessment monitoring stations,光化學評估監(jiān)測站網(wǎng)絡),積累了豐富的O3污染案例,為評估O3及其前體物污染演變態(tài)勢、深入分析O3生成機制提供了有力的數(shù)據(jù)支撐.

    基于排放的空氣質(zhì)量模型綜合考慮大氣物理過程(氣象場)和大氣化學過程(大氣化學機制),運用數(shù)學方法結合排放清單模擬已知復雜的大氣污染過程,被廣泛用于建立區(qū)域O3與前體物減排的響應關系.然而,排放清單和邊界層動力學參數(shù)化具有不確定性,在很大程度上會影響模擬效果的準確度. 在此背景下,OBM應運而生,運用數(shù)學方法和計算機語言并結合實際大氣環(huán)境中O3和前體物的觀測濃度來診斷O3生成機制. OBM是美國佐治亞理工學院Cardelino等[11]于1995年依據(jù)相對增量反應活性概念RIR(relative incremental reactivity)開發(fā)的一個盒子模型.運行過程中,OBM將VOCs、NO、CO和O3的觀測濃度輸入盒子模型中,假設污染物充分混合,利用CBM-Ⅳ(第4代carbon bond mechanism)模擬原位大氣光化學反應過程,然后通過降低特定前體物的濃度(表征污染源的削減)來計算該前體物類別的RIR值,并據(jù)此判斷O3的生成對NOx以及不同VOCs類別的敏感性. RIR是應用OBM診斷O3與前體物關系的核心概念,其定義為目標前體物濃度變化的百分比與所引起的O3變化(生成量或凈生成速率)的百分比的比值[10]. RIR為正值,表明削減該前體物類別能有效降低O3的生成,且數(shù)值越大,表明O3生成對該前體物類別越敏感,即為活性物種;相反,RIR為負值則表明,控制該前體物類別會促進O3的生成(如NO滴定效應).

    與基于排放的空氣質(zhì)量模型相比,OBM的最大優(yōu)勢在于避免了排放清單和邊界層動力學參數(shù)化的不確定性所導致的誤差,能夠更加客觀地反映實際大氣的污染狀況以及O3與其前體物的關系. 但是,OBM也存在自身的局限性:一方面,它是一種診斷性的工具,不具備預測污染發(fā)生的功能;另一方面,觀測數(shù)據(jù)較粗的空間分辨率使得OBM的空間代表性有限,不適合用于解析耦合不同空間尺度的大氣化學過程.

    1.2 OBM結構的發(fā)展歷程

    OBM的基本原理是運用數(shù)學方法結合實際觀測數(shù)據(jù)來模擬解析大氣化學反應速率. 該原理間接規(guī)定了應用OBM的兩個前提,一是運用數(shù)學方法和計算機語言(如Matlab、IDL、Fortran和Facsimile)構建OBM運行系統(tǒng)[12],二是提供高質(zhì)量的觀測數(shù)據(jù)來約束OBM的積分計算. 對于前者,OBM的本質(zhì)是化學盒子模型,其運行系統(tǒng)的構建(主要包括初始化、觀測數(shù)據(jù)實時約束、大氣物理、大氣化學和輸出等模塊)是對化學盒子模型結構的進一步改進和延伸;對于后者,觀測數(shù)據(jù)的精確性與完整性越高,大氣化學過程的模擬量化越準確.

    如圖1所示,早期的化學盒子模型主要用于煙霧箱模擬試驗的對比驗證,結構相對簡單,一般包括初始化模塊(設置反應物的初始濃度和環(huán)境參數(shù))、大氣化學模塊(設置大氣化學機制)和輸出模塊(設置輸出反應物和產(chǎn)物的濃度及特定過程反應速率).Cardelino等[11]構建的OBM在改進以上模塊功能的基礎上增加了觀測數(shù)據(jù)實時約束模塊和大氣物理模塊. 其中,觀測數(shù)據(jù)實時約束模塊用于讀入化學成分的觀測數(shù)據(jù)以實時約束模型的積分計算;大氣物理模塊用于參數(shù)化關鍵的物理過程,如太陽輻射變化以及溫度、相對濕度和邊界層高度的日變化. 初始化模塊在初始化反應物濃度和氣象參數(shù)的基礎上還需設置觀測站點地理位置和模擬時間,為大氣物理模塊參數(shù)化大氣物理過程提供必要的信息. 隨著OBM應用需求的變化,大氣物理模塊的功能逐漸豐富,對干沉降和擴散等物理過程也進行了參數(shù)化.

    大氣化學模塊和輸出模塊的改進是升級OBM結構的核心. 相比于早期的化學盒子模型,Cardelino等[11]構建的OBM輸出模塊中增加了敏感性/情景分析功能,用于計算RIR值以診斷O3與前體物的關系.當前階段,大氣化學機制的不斷完善(特別是液相大氣化學機制的出現(xiàn))使得RIR的定義不再局限于O3污染,而被更新為目標前體物濃度變化的百分比與所引起的二次污染物變化(生成量或凈生成速率)的百分比的比值,可以用于診斷O3或二次氣溶膠與前體物的關系. 通過對初始化、觀測數(shù)據(jù)實時約束、大氣物理、大氣化學和輸出模塊的設置,OBM綜合考慮了重要的大氣物理和化學過程,可以建立任一活性組分在任何時刻的光穩(wěn)態(tài)關系,因此可以用于化學反應機制研究、高活性物種濃度模擬和化學過程分析等.

    考慮到O3污染問題既是OBM應用的熱點,也是當前國內(nèi)空氣質(zhì)量管理的重點,筆者以應用OBM診斷O3與前體物的關系為例進一步闡明OBM各模塊的具體設置. 常規(guī)氣體(如O3、NO、NO2、CO和SO2)和VOCs組分的觀測濃度以及氣象參數(shù)(溫度、相對濕度和氣壓)是OBM模擬的必需輸入數(shù)據(jù),其他數(shù)據(jù)可選擇性輸入. 首先進行基準情景的模擬(即無需對實際觀測數(shù)據(jù)進行任何調(diào)整):初始化模塊中設置污染物的初始濃度及氣象參數(shù)、OBM運行的起始時間和積分步長以及觀測站點的信息(經(jīng)緯度、海拔)等. 通常設置OBM運行的起始時間早于研究時段3~5天(模型預運行),使未被約束的化學物種(如自由基)達到穩(wěn)態(tài),為后續(xù)模擬提供更真實的大氣化學環(huán)境;觀測數(shù)據(jù)實時約束模塊按照一定的時間分辨率將O3、NO、NO2、CO、SO2和VOCs組分的觀測濃度以及溫度、相對濕度和氣壓讀入模型中. 需要注意的是,輸入數(shù)據(jù)必須是連續(xù)的,否則后續(xù)積分運算難以連續(xù)進行;大氣物理模塊中關于干沉降、擴散和太陽輻射變化的參數(shù)化方式并沒有統(tǒng)一的標準,該節(jié)不再詳細闡述;大氣化學模塊可根據(jù)模擬需求選取特定氣相化學機制進行應用,其中歸納化學機制需要對輸入的VOCs組分提前進行歸納處理;輸出模塊設置輸出O3的凈生成速率,以便后續(xù)計算RIR值.

    同理,還需要進行一系列減排情景的模擬以量化不同前體物類別的RIR值. 與基準情景的模擬相比,減排情景只需將某特定前體物類別(或單個前體物物種)降低一定比例(如10%),而其他條件保持不變.如式(1)所示,通過基準情景和減排情景量化的O3凈生成速率的比較,可以計算該前體物類別的RIR值[13]. 需要注意的是,由于O3與其前體物呈非線性關系,RIR會隨著減排比例的不同而發(fā)生變化.

    式中:RIR(X)為相對增量反應活性;Pnet-O3(X)為基準情景下的O3凈生成速率,10-9/h;Pnet-O3(X-ΔX)為減排情景下的O3凈生成速率,10-9/h;X為某個特定前體物組分;ΔS(X)為減排情景下假設的前體物濃度的變化,10-9;S(X)為基準情景下前體物組分的濃度,10-9.

    1.3 大氣化學機制的發(fā)展歷程

    大氣化學機制是OBM的重要組成部分,在很大程度上決定了OBM的應用范圍和推廣程度. 石玉珍等[14-15]已經(jīng)總結了幾種相對成熟的大氣化學機制的發(fā)展歷程,因此,該節(jié)簡要闡述大氣化學機制的發(fā)展歷程,并重點介紹OBM中常用大氣化學機制的適用條件及模擬結果對比.

    大氣化學機制的發(fā)展貫穿于OBM整個發(fā)展歷程(見圖1). 早期的化學盒子模型通常根據(jù)需要將與研究物種相關的已知化學反應參數(shù)化作為大氣化學模塊,隨著研究的不斷深入,先后出現(xiàn)了幾種比較成熟的大氣化學機制,涉及氣相反應〔如CBM、RACM(Regional Atmospheric Chemical Mechanism)、SAPRC(Statewide Air Pollution Research Center Mechanism)

    圖 1 OBM基本結構的演變Fig.1 The evolution of basic structures of OBM

    和MCM(Master Chemical Mechanism)等〕與液相反應〔如CAPRAM(Chemical Aqueous Phase RAdical Mechanism)等〕. 1995年,Cardelino等[11]開發(fā)的OBM利用CBM-Ⅳ模擬原位大氣光化學反應過程. CBM由Systems Applications International于1976年開發(fā),隨后數(shù)十年經(jīng)歷多次更新,主要版本有CBM-Ⅱ、CBM-Ⅲ、CBM-EX、CBM-Ⅳ(1989年)、CB05(2005年)和CB6(2010年)[16];RACM是在區(qū)域酸沉降機理RADM基礎上發(fā)展起來的,主要版本有RADM(1986年 )、RADM2 (1990年 )、RACM (1997年 )和RACM2(2012年)[17];SAPRC由Carter于1990年開發(fā),主要版本有SAPRC90(1990年)、SAPRC93(1993年)、SAPRC99(2000年)、SAPRC07(2010年)和SAPRC18(2020年)[18];MCM由Jenkin等于1997年開發(fā),主要版本有MCMv1(1996年)、MCMv2(1999年)、MCMv3(2002年)、MCMv3.1(2003年)、MCMv3.2(2012年)和MCMv3.3.1(2015年)[19-20]. 上述氣相化學機制詳細的發(fā)展歷程可參見文獻[14-15]. 液相化學機制CAPRAM由Herrmann等于1999年開發(fā),主要版本有CAPRAM2.3(2000年)、CAPRAM2.4(2003年)和CAPRAM3.0(2005年)[21]. 截至目前,以上大氣化學機制仍在持續(xù)更新與完善,為OBM的進一步發(fā)展提供了可靠的技術支撐.

    如表1所示,上述常用的大氣化學機制可以分為歸納化學機制(lumped mechanism)和詳細化學機制(explicit mechanism). 其中,氣相歸納化學機制主要包括CBM、SAPRC和RACM等,液相歸納化學機制主要包括CAPRAM等. 歸納化學機制涉及的化學反應式較少,因此計算效率較高,處理起來靈活方便,得以在模型中廣泛應用. 但是,歸納化學機制通常需要對VOCs進行歸類處理,會引發(fā)不確定性,如CBM根據(jù)分子結構類型對VOCs進行歸類,而SAPRC和RACM主要根據(jù)VOCs與OH·的反應活性進行歸類.相比之下,詳細化學機制(如MCM等)直接描述了單個VOC物種的大氣化學反應過程,避免了VOCs歸類造成的不確定性,但是同時增大了計算需求量,一定程度上限制了其在模型中的應用. 此外,OBM可以將氣相和液相化學機制聯(lián)用(如RACM/MCMCAPRAM等),內(nèi)部建立一套氣-液傳質(zhì)過程連通氣相和液相之間的物種轉化,適用于研究二次氣溶膠的生成機制. 需要注意的是,上述大氣化學機制均由國外研究團隊開發(fā),目前我國仍缺乏基于本國大氣化學環(huán)境開發(fā)的化學機制.

    化學機制內(nèi)部百家爭鳴的現(xiàn)象也使得基于不同化學機制(特別是不同氣相化學機制)的模擬結果對比成為熱點問題. 整體上,不同化學機制模擬結果(如O3、HOx·和PAN等污染物的濃度和凈生成速率)的日變化比較一致,但是量級存在不同程度的差距. Coates等[22-24]的研究結果表明,模擬情景中NOx的可利用性對于化學機制的模擬結果有較大影響. 低VOCs/NOx情景(即污染情景)下,活性氮化學機制占據(jù)主導地位,此時基于不同化學機制得到的模擬結果比較接近;相比之下,高VOCs/NOx情景(即清潔情景)下,基于不同化學機制得到的模擬結果則相差較大,這主要是由于化學機制內(nèi)部對有機化學(特別是VOCs物種)歸納和動力學參數(shù)的處理方式不同,很難判定究竟基于哪種化學機制的模擬結果更準確. 鑒于幾種比較成熟的大氣化學機制仍在持續(xù)發(fā)展與更新,建議在同類大氣化學機制中優(yōu)先選擇最新發(fā)布的版本進行應用. 準確的化學機制是開展模擬分析工作的必要前提,因此仍需進一步完善大氣化學機制.

    表 1 OBM中常用的大氣化學機制概述Table 1 Summary of atmospheric chemical mechanisms built in OBM

    2 OBM在中國大氣化學研究中的應用

    日臻完善的空氣質(zhì)量監(jiān)測網(wǎng)絡為應用OBM解析大氣化學過程提供了數(shù)據(jù)支撐,針對我國區(qū)域大氣復合污染成因開展了大量卓有成效的研究工作,取得了一系列重要進展.

    2.1 臭氧生成機制研究

    近年來日益嚴峻的O3污染形勢是OBM得以在國內(nèi)獲得關注的一個重要推力,同時O3生成機制的研究也是OBM應用的熱點. 筆者根據(jù)“observationbased model”“box model”“ozone”“O3”以及“China”“Chinese”“Pearl River Delta”“Yangtze River Delta”“North China Plain”“Beijing-Tianjin-Hebei”“Sichuan”“Beijing”“Shanghai”“Guangzhou”“Hong Kong”等檢索關鍵詞的組合、依托Web of Science、Scopus和Google Scholar等數(shù)據(jù)庫統(tǒng)計了2000-2021年9月公開發(fā)表的利用OBM開展國內(nèi)O3生成機制研究的英文論文數(shù)量(見圖2). 總體上,應用OBM診斷O3生成機制的論文數(shù)量隨時間尺度呈增加趨勢,特別是2017-2021年〔平均值±標準偏差為(13.6±4.9)篇〕比2007-2016年〔(3.6±1.9)篇〕高出10篇/a. 這一結果體現(xiàn)出,一方面,近年來我國對O3污染問題的重視程度與日俱增;另一方面,VOCs組分觀測數(shù)據(jù)的增多為OBM的推廣和應用提供了可靠的數(shù)據(jù)支撐. 在此基礎上,進一步根據(jù)研究區(qū)域對論文進行分類. 從空間分布來看,大部分OBM研究集中在京津冀及周邊、長三角和珠三角等經(jīng)濟相對發(fā)達的東部地區(qū)(77%),而廣大的中、西部地區(qū)研究相對較少(23%),且起步較晚,OBM還有比較廣闊的應用前景.

    圖 2 2000?2021年9月公開發(fā)表的利用OBM開展國內(nèi)O3生成機制研究的英文論文數(shù)量Fig.2 The number of papers (in English) regarding OBM studies in O3 formation regimes over China

    2.1.1臭氧生成敏感性的時空分布

    應用OBM診斷O3生成敏感性的研究成果[9,25]總結見表2. 不同區(qū)域O3生成敏感性的空間分布特征既存在一致也存在差異. 京津冀及周邊、長三角、珠三角和其他區(qū)域的城市站點O3生成敏感性基本一致,即主要處在VOCs控制區(qū);郊區(qū)站點的O3生成敏感性存在差異,如長三角和珠三角以及成渝地區(qū)(如重慶市)的郊區(qū)站點主要處在VOCs控制區(qū),而京津冀及周邊、西北地區(qū)(如蘭州市)和長江中下游地區(qū)(如武漢市)的郊區(qū)站點則主要處在NOx-VOCs協(xié)同控制區(qū)或NOx控制區(qū);鄉(xiāng)村站點的O3生成敏感性也存在差異,如長三角和珠三角的鄉(xiāng)村站點主要處在VOCs控制區(qū)或NOx-VOCs協(xié)同控制區(qū),而京津冀及周邊和成渝地區(qū)(如重慶市)的鄉(xiāng)村站點則主要處在NOx控制區(qū)或NOx-VOCs協(xié)同控制區(qū). 值得注意的是,目前發(fā)表的利用OBM診斷O3生成敏感性的研究大部分集中在城市地區(qū),郊區(qū)和鄉(xiāng)村地區(qū)則相對較少,這主要是由觀測試驗在不同功能區(qū)分布不均勻所致.對于特定區(qū)域(如京津冀及周邊地區(qū)),從城市到郊區(qū)再到鄉(xiāng)村站點(NOx濃度水平通常呈遞減趨勢),O3生成對于NOx的敏感性通常呈遞增趨勢,該結論與基于三維空氣質(zhì)量模型獲得的O3生成敏感性在水平尺度的空間分布[26-27]基本一致.

    表 2 基于OBM開展的O3生成機制研究總結[9,25]Table 2 Summary of the observation-based model studies in O3 formation regime over China[9,25]

    利用OBM探究垂直尺度上O3生成敏感性變化的研究相對較少. Benish等[28]利用OBM結合飛機航測數(shù)據(jù)分析河北省O3生成敏感性的垂直分布特征,發(fā)現(xiàn)靠近城市中心的低海拔地區(qū)處于VOCs控制區(qū),而隨著海拔的升高,O3生成對NOx的敏感性增加,逐漸轉變?yōu)镹Ox-VOCs協(xié)同控制區(qū)甚至NOx控制區(qū).這主要是由于NO的化學壽命非常短(<10 s),其濃度隨著海拔的升高而迅速降低,相比之下,部分VOCs物種壽命較長,甚至可以傳輸?shù)降妥杂蓪α鲗?約3 500 m),使得高空VOCs/NOx值較高,NOx成為限制O3生成的重要前體物. 該結論與基于三維空氣質(zhì)量模型獲得的O3生成敏感性在垂直尺度的空間分布特征[29]基本一致.

    在時間尺度上,Lu等[30-31]利用OBM揭示了O3生成敏感性的日變化特征. 早晨,光化學反應處于起步階段,自由基含量較低,O3生成主要由VOCs控制,而下午光化學反應較為旺盛,NOx被快速消耗,O3生成對NOx的敏感性增加,逐漸轉變?yōu)橛蒒Ox-VOCs共同控制或者由NOx控制. 上述日變化特征在城市(北京市北京大學站點)、郊區(qū)(北京市榆垡站點)和鄉(xiāng)村(長三角地區(qū)某一站點)大氣環(huán)境下均存在,為O3削峰行動提供了新思路.

    與OBM結果一致,基于觀測的光化學指示劑法〔如臭氧生成效率(OPE)、H2O2/NOz(或H2O2/HNO3)和HCHO/NO2比值〕的結果表明,我國典型地區(qū)大氣環(huán)境中O3生成敏感性存在明顯的時空分布特征[9]. 因此,僅基于單一站點的觀測結果來制定區(qū)域O3污染防控措施的策略可能會存在較大的不確定性,未來應建立有代表性的區(qū)域監(jiān)測網(wǎng)絡(建議包括常規(guī)六參數(shù)、詳細的VOCs組分和氣象參數(shù)等)以獲取高時空分辨率的觀測數(shù)據(jù),結合OBM或光化學指示劑法綜合反映區(qū)域尺度的O3生成敏感性.

    2.1.2控制臭氧生成的關鍵活性VOCs物種

    利用OBM計算主要VOCs類別(或單個VOC物種)的相對增量反應活性(RIR)或增量反應活性(incremental reactivity,IR )可以進一步識別關鍵活性VOCs物種.

    基于OBM量化RIR指標可以直接判斷優(yōu)控VOCs物種. 例如:Geng等[32-34]證實上海市的優(yōu)控物種主要是芳香烴;An等[35]證實南京市的優(yōu)控物種主要是烯烴;而珠三角地區(qū)控制O3生成的活性VOCs物種主要是芳香烴[9]. 個別站點(如禹城鄉(xiāng)村站點)植被覆蓋率較高,夏季高溫條件下會排放大量的生物源VOCs,從而在特定時間段對本地O3的生成產(chǎn)生重要影響[36].

    利用OBM計算IR指標可以用于量化臭氧生成潛勢( ozone formation potential,OFP ),進 而 估 算VOCs物種對O3生成的貢獻. Qiu等[37]和Zhang等[38]分別利用OBM耦合不同的大氣化學機制(前者使用RACM,后者使用MCM)計算了中國城市(如北京、上海、成都、深圳和廣州)本土化的MIR值(適用于高NOx情景),發(fā)現(xiàn)烯烴和芳香烴的MIR值較大,活性相對較高. 將上述結果與Carter[18]利用SAPRC針對美國城市污染情景計算出的MIR值進行對比,發(fā)現(xiàn)基于中、美大氣環(huán)境計算出的MIR值存在差異,一方面是由中、美大氣環(huán)境的差異(如VOCs化學組成方面,中國環(huán)境空氣中烯烴和OVOCs占比較大而烷烴占比較小)所致,另一方面是由不同的大氣化學機制(RACM、MCM、SAPRC)所致. 對于特定區(qū)域,直接應用其他地區(qū)的MIR值量化OFP可能會引發(fā)較大的不確定性.

    需要注意的是,基于OFP與RIR指標識別出的關鍵活性VOCs物種可能會存在差異,這是由于RIR指標根據(jù)實際NOx情景判斷O3生成對VOCs物種的敏感性,而OFP指標是在特定NOx情景(如MIR、MOR和EBIR)下估算VOCs物種經(jīng)過氧化生成O3的潛力. 基于RIR的結果更符合實際情景,但相比于RIR指標,OFP計算方式較為簡便. 目前我國大多數(shù)城市處于VOCs控制區(qū)且缺乏可用的本土化MIR值,考慮到大氣環(huán)境對量化MIR指標的影響以及當前國內(nèi)較為豐富的VOCs組分觀測數(shù)據(jù),應積極推進利用OBM計算城市MIR值的工作,為“十四五”期間VOCs的精細化管控提供科學支撐.

    2.2 二次氣溶膠生成機制研究

    OBM可以用于診斷二次氣溶膠的生成機制,但是相比于O3研究起步較晚. Xue等[39]利用OBM揭示了溶解性NO2氧化S(Ⅳ)和SO2在氣溶膠表面的非均相反應是我國特大城市霧霾條件下SO42-(硫酸鹽)的主要生成途徑,控制NOx和SO2排放可以有效抑制SO42-生成. 少數(shù)研究基于綜合觀測數(shù)據(jù)和OBM迭代模型計算了北京城市和鄉(xiāng)村以及中國香港鄉(xiāng)村站點大氣中的N2O5非均相反應攝取系數(shù),并將修正后的攝取系數(shù)直接用于OBM模擬,量化了N2O5非均相反應對夜間NO3-(硝酸鹽)生成的重要貢獻[40-42].Pathak等[43]利用OBM耦合RACM-CAPRAM得到了一致的結論,其中北京鄉(xiāng)村和上海郊區(qū)大氣中N2O5非均相反應對夜間NO3-的貢獻率為50%~100%.Wang等[44]在北京城市站點開展了垂直尺度上的觀測試驗,并結合OBM量化了N2O5非均相反應在不同高度對NO3-的貢獻,發(fā)現(xiàn)近地面N2O5非均相反應對NO3-的貢獻較小(高NO濃度環(huán)境使得N2O5濃度較低),而高層(如>150 m)大氣的貢獻較大,并會在次日清晨通過垂直混合加重近地面NO3-污染. Wen等[45]利用OBM耦合RACM-CAPRAM揭示了華北平原不同大氣環(huán)境(城市、鄉(xiāng)村和背景)中NO3-生成機制的共性,即白天主要以HNO3的氣相生成以及非均相傳質(zhì)為主(貢獻率為94%~96%),其中NO2是控制NO3-生成的關鍵前體物;夜間則主要以N2O5非均相反應為主(貢獻率為91%~98%),受NO2和O3共同控制. 以上研究揭示了我國SO42-和NO3-等主要二次無機離子組分的生成機制,對于制定科學的二次氣溶膠防控對策具有重要的參考價值.

    此外,Zhu等[46]采用國際上最詳細的MCMCAPRAM多相化學盒子模型并結合排放數(shù)據(jù)定量評估了液相反應對泰山大氣中二元羧酸等二次氣溶膠生成的重要貢獻(占比超過70%),揭示了控制液相二元羧酸生成的關鍵前體物以及乙二醛的氣-粒分配過程對二元羧酸生成的重要影響. 其中,異戊二烯、1,3-丁二烯和甲苯是乙二酸和乙醛酸的關鍵前體物,而異戊二烯是丙酮酸的關鍵前體物. 需要注意的是,當前階段應用OBM解析二次有機氣溶膠生成機制的研究仍相對較少,這主要是受到觀測數(shù)據(jù)可利用性的制約. 未來需提高對二次有機氣溶膠組分的觀測能力,為應用OBM深入解析二次有機氣溶膠生成機制提供可靠的數(shù)據(jù)支撐.

    2.3 其他活性組分化學機制研究

    厘清大氣自由基化學是理解區(qū)域復合污染成因的基礎和關鍵. 利用OBM量化自由基(OH·、HO2·和RO2·)初級來源的研究結果[47]如表3所示. OVOCs光解是HO2·和RO2·重要的初級來源. 不同大氣環(huán)境中OH·主要的初級來源不盡相同. Yang等[48-49]揭示了O3光解是城市站點(如北京和蘭州)OH·主要的初級來源;Yang等[50-51]揭示了HONO光解是郊區(qū)站點(如成都和中國香港) OH·主要的初級來源;不同區(qū)域鄉(xiāng)村站點OH·的初級來源存在差異,如O3光解是東營和中國香港鄉(xiāng)村站點OH·主要的初級來源[52-53],而HONO光解是泰山OH·主要的初級來源[54];O3光解是瓦里關背景站點OH·主要的初級來源[55]. 上述研究中多數(shù)缺乏自由基觀測數(shù)據(jù),僅基于OBM量化的自由基循環(huán)與收支. Lu等[56-57]將OH·觀測數(shù)據(jù)與OBM模擬結果進行對比,發(fā)現(xiàn)基于現(xiàn)有大氣化學機制模擬的OH·濃度被明顯低估,并進一步證實了過氧自由基參與的非傳統(tǒng)再生機制是低NOx情景下OH·主要的再生來源,其中異戊二烯有機過氧自由基的氫轉移反應貢獻顯著. 該結論在高VOCs和低NOx環(huán)境下具有普遍性,為大氣化學機制的進一步完善提供了新思路.

    表 3 基于OBM量化自由基初級來源和OH·反應活性的研究總結[47]Table 3 Summary of the observation-based model studies of radical sources over China[47]

    HONO光解是OH·重要的初級來源,在大氣化學中扮演著重要角色. 目前應用OBM解析HONO化學機制的研究主要關注其來源分析和HONO對自由基及大氣氧化能力的影響. 對于前者,Liu等[58-60]參數(shù)化潛在的HONO生成途徑(如土壤排放、顆粒態(tài)硝酸鹽光解和NO2水解的光增強反應等),并將其添加至現(xiàn)有大氣化學機制,同時結合OBM來量化HONO的來源,顯著改善了模型模擬HONO的能力,但是由于反應參數(shù)未定,敏感性試驗表明獲得的結果具有較大的不確定性. 對于后者,Jiang等[54]和Li等[61]利用OBM分別模擬了泰山和北京大氣環(huán)境中有、無HONO觀測數(shù)據(jù)約束兩種情景下的自由基收支與大氣氧化能力,定量評估了HONO光解對大氣OH·的初級來源貢獻,揭示了未知HONO來源對于定量OH·初級生成速率、OH·和HO2·濃度及大氣氧化能力的重要影響.

    活性鹵素自由基會氧化VOCs,推動大氣自由基循環(huán),從而促進O3和二次氣溶膠的生成. 目前應用OBM解析鹵素化學機制的研究主要關注鹵素大氣化學機制的完善和鹵素化學對空氣質(zhì)量的影響.早期研究對鹵素活化過程的認識不夠充分,所以當前的大氣化學反應機制很少或基本不考慮鹵素自由基與VOCs的化學反應,從而可能低估了大氣氧化能力. Xue等[62]和Peng等[63]先后更新了大氣氯化學模塊和溴化學模塊,并將其添加至現(xiàn)有大氣化學機制MCM,結合OBM模型揭示了Cl·和Br·對ROx· (ROx·=OH·+HO2·+RO2·)循環(huán)和O3生成的顯著促進作用. 以上我國自主開發(fā)的大氣氯化學模塊和溴化學模塊已經(jīng)在國際上得到證實與認可,并被國內(nèi)外多個課題組用于活性鹵素化學研究.

    3 問題與展望

    筆者回顧了OBM的發(fā)展歷程,總結了利用OBM解析國內(nèi)O3和二次氣溶膠生成機制及其他活性組分化學機制的相關研究成果. 當前,我國PM2.5污染防控工作進入瓶頸期,O3污染仍呈上升和蔓延態(tài)勢,區(qū)域大氣復合污染防控已成為空氣質(zhì)量管控的重大挑戰(zhàn). 日臻完善的空氣質(zhì)量監(jiān)測網(wǎng)絡較為全面地表征了我國大氣污染水平,同時也為利用OBM基于實際大氣環(huán)境解析區(qū)域大氣復合污染成因打下了堅實的基礎. 針對OBM的應用現(xiàn)狀,筆者指出如下問題并提出相應的建議,以供OBM進一步發(fā)展與完善.

    3.1 存在的問題

    a) OBM的結構決定了其在大氣化學研究中應用的局限性. OBM內(nèi)置的大氣化學機制較為詳細,但是基本不考慮或很少考慮大氣物理因素. 因此,OBM更多的是作為一種解析大氣化學過程的診斷性工具,而不具備預測大氣污染發(fā)生的功能. OBM不是對基于排放的空氣質(zhì)量模型的完全取代,而是在日臻完善的空氣質(zhì)量監(jiān)測網(wǎng)絡基礎上產(chǎn)生的對空氣質(zhì)量模型的補充.

    b) OBM的模擬效果直接取決于觀測數(shù)據(jù)的質(zhì)量. 目前我國可利用的空氣質(zhì)量觀測數(shù)據(jù)仍有限,主要表現(xiàn)在觀測數(shù)據(jù)空間分布不均勻(集中在經(jīng)濟發(fā)達的東部地區(qū))以及可公開獲取的VOCs和顆粒物組分數(shù)據(jù)較少,一定程度上限制了OBM的應用和推廣.特別地,現(xiàn)有研究中二次有機氣溶膠的觀測數(shù)據(jù)時間分辨率一般較低,不適合與OBM結合進行詳細的化學機制解析.

    c) 我國尚未開發(fā)出具有自主知識產(chǎn)權的OBM.一方面,大氣化學機制是OBM發(fā)展的核心,但是目前數(shù)值模型中常用的比較成熟的大氣化學機制(如CBM、RACM、SAPRC、MCM和CAPRAM)均由國外研究團隊開發(fā),我國仍缺乏基于本國大氣化學環(huán)境開發(fā)的大氣化學機制;另一方面,我國缺乏自主開發(fā)的開源OBM運行平臺,大多需要付出一定的經(jīng)濟成本從國外購買運行平臺(如Matlab、Facsimile和IDL)的許可證.

    3.2 未來展望

    a) OBM是一種強大的大氣污染成因研究與化學機理解析工具,其最大優(yōu)勢在于,基于觀測數(shù)據(jù)約束計算實際大氣環(huán)境的化學反應速率,避免了排放清單帶來的不確定性. 在實際應用中應注意揚長避短,充分利用OBM的優(yōu)勢,根據(jù)大氣化學過程解析需求選擇合適的模型.

    b) 突破制約OBM發(fā)展的關鍵技術. 一方面,要開發(fā)出具有中國自主知識產(chǎn)權的在線OBM運行系統(tǒng),對于充分挖掘觀測數(shù)據(jù)價值、快速診斷O3和二次氣溶膠生成機制具有重要意義;另一方面,要進一步完善、細化詳細大氣化學機制,發(fā)展本土化的歸納化學機制,開發(fā)具有中國自主知識產(chǎn)權的大氣化學機制.

    c) 不斷完善我國的大氣監(jiān)測網(wǎng)絡體系,使其覆蓋空間更廣泛(由東部地區(qū)向中西部地區(qū)延伸)、包含物種更全面(組建全國大氣顆粒物組分監(jiān)測網(wǎng)和大氣光化學監(jiān)測網(wǎng))、運行更高效(實現(xiàn)自動程序化監(jiān)測),為應用OBM解析區(qū)域大氣復合污染成因提供可靠的數(shù)據(jù)支撐.

    猜你喜歡
    前體大氣觀測
    觀測到恒星死亡瞬間
    軍事文摘(2023年18期)2023-11-03 09:45:42
    大氣的呵護
    軍事文摘(2023年10期)2023-06-09 09:15:06
    N-末端腦鈉肽前體與糖尿病及糖尿病相關并發(fā)癥呈負相關
    天測與測地VLBI 測地站周圍地形觀測遮掩的討論
    大氣古樸揮灑自如
    N-端腦鈉肽前體測定在高血壓疾病中的應用研究
    大氣、水之后,土十條來了
    可觀測宇宙
    太空探索(2016年7期)2016-07-10 12:10:15
    高分辨率對地觀測系統(tǒng)
    太空探索(2015年8期)2015-07-18 11:04:44
    茶葉香氣前體物研究進展
    茶葉通訊(2014年2期)2014-02-27 07:55:40
    成人无遮挡网站| 赤兔流量卡办理| 九草在线视频观看| 美女视频免费永久观看网站| 人妻一区二区av| 国产成人精品久久久久久| 一区在线观看完整版| 黑人高潮一二区| 18禁在线无遮挡免费观看视频| 免费看日本二区| 日韩成人伦理影院| 亚洲精品乱码久久久久久按摩| 黑人高潮一二区| 日韩欧美精品免费久久| 精品99又大又爽又粗少妇毛片| 亚洲无线观看免费| 日日啪夜夜撸| 色哟哟·www| 大香蕉97超碰在线| 国产免费又黄又爽又色| 中文天堂在线官网| a级毛片免费高清观看在线播放| 老司机亚洲免费影院| 久久久久久久大尺度免费视频| 欧美老熟妇乱子伦牲交| 国产精品福利在线免费观看| 99热国产这里只有精品6| 人人澡人人妻人| 乱系列少妇在线播放| 丁香六月天网| 欧美xxⅹ黑人| 色吧在线观看| 老司机亚洲免费影院| 老女人水多毛片| h日本视频在线播放| 如日韩欧美国产精品一区二区三区 | 日本黄色日本黄色录像| 黄片无遮挡物在线观看| 亚洲情色 制服丝袜| 在线免费观看不下载黄p国产| 大又大粗又爽又黄少妇毛片口| 全区人妻精品视频| 国产精品秋霞免费鲁丝片| 在线观看国产h片| 51国产日韩欧美| 亚洲三级黄色毛片| 女性生殖器流出的白浆| 国产伦精品一区二区三区视频9| 一区二区三区四区激情视频| 国产91av在线免费观看| 99精国产麻豆久久婷婷| 97超视频在线观看视频| 婷婷色综合www| videos熟女内射| 亚洲av二区三区四区| 少妇被粗大猛烈的视频| 热re99久久国产66热| 久久久久久久国产电影| 国产精品.久久久| 日韩中字成人| 成人无遮挡网站| 丁香六月天网| 精品人妻偷拍中文字幕| 亚洲情色 制服丝袜| 热re99久久精品国产66热6| 国产又色又爽无遮挡免| 欧美日韩视频精品一区| 免费观看性生交大片5| 国产精品久久久久久精品古装| a级片在线免费高清观看视频| 亚洲精品视频女| 亚洲欧洲日产国产| 国产精品伦人一区二区| 女性生殖器流出的白浆| 亚洲精品乱久久久久久| 美女内射精品一级片tv| 日韩 亚洲 欧美在线| 在线亚洲精品国产二区图片欧美 | 亚洲av福利一区| √禁漫天堂资源中文www| 精品久久久久久电影网| 一本—道久久a久久精品蜜桃钙片| 国产成人精品婷婷| av天堂久久9| 久久久久久久精品精品| 久久综合国产亚洲精品| 亚洲色图综合在线观看| av国产久精品久网站免费入址| 99九九在线精品视频 | 高清午夜精品一区二区三区| 午夜激情久久久久久久| 99久国产av精品国产电影| 精品亚洲成a人片在线观看| 亚洲在久久综合| 如日韩欧美国产精品一区二区三区 | 亚洲av.av天堂| 性色avwww在线观看| 亚洲av成人精品一区久久| 亚洲欧美一区二区三区黑人 | 国产精品免费大片| 日韩熟女老妇一区二区性免费视频| 中国美白少妇内射xxxbb| 女性生殖器流出的白浆| 免费人妻精品一区二区三区视频| 午夜福利网站1000一区二区三区| 少妇高潮的动态图| 日本免费在线观看一区| 亚洲国产精品一区二区三区在线| 亚洲第一区二区三区不卡| 婷婷色麻豆天堂久久| 国产精品福利在线免费观看| 大片电影免费在线观看免费| 精品一区在线观看国产| 午夜福利视频精品| 老司机影院成人| 一级黄片播放器| 国产日韩欧美亚洲二区| 国产精品国产三级国产av玫瑰| 成人18禁高潮啪啪吃奶动态图 | 国产 一区精品| 亚洲内射少妇av| 51国产日韩欧美| 久久久久久久久久成人| 黑人高潮一二区| 人妻少妇偷人精品九色| 国产一区二区三区综合在线观看 | 亚洲精品,欧美精品| 免费观看性生交大片5| 亚洲av在线观看美女高潮| 女性生殖器流出的白浆| 欧美日韩精品成人综合77777| 精品人妻熟女毛片av久久网站| 日韩一区二区三区影片| 国产伦在线观看视频一区| 亚洲激情五月婷婷啪啪| 免费久久久久久久精品成人欧美视频 | 美女福利国产在线| 亚洲国产最新在线播放| 亚洲精品成人av观看孕妇| 久久人人爽人人片av| 丰满饥渴人妻一区二区三| 欧美 日韩 精品 国产| 伊人久久国产一区二区| 80岁老熟妇乱子伦牲交| 日本黄色日本黄色录像| 大陆偷拍与自拍| 免费高清在线观看视频在线观看| 亚洲美女视频黄频| a级毛片免费高清观看在线播放| 亚洲成人手机| a级毛色黄片| 国产成人精品一,二区| 欧美激情国产日韩精品一区| 22中文网久久字幕| 少妇被粗大猛烈的视频| 人人妻人人添人人爽欧美一区卜| 下体分泌物呈黄色| 国产精品人妻久久久久久| 日日啪夜夜撸| 亚洲国产精品专区欧美| 成年人午夜在线观看视频| 另类亚洲欧美激情| 桃花免费在线播放| 欧美少妇被猛烈插入视频| 九草在线视频观看| 国产成人精品久久久久久| 久久久久久伊人网av| 国产男人的电影天堂91| 中国国产av一级| 日韩 亚洲 欧美在线| 一级毛片aaaaaa免费看小| 成人毛片60女人毛片免费| 在线观看三级黄色| 婷婷色综合www| 在线观看人妻少妇| 国产男女超爽视频在线观看| 成人国产麻豆网| 亚洲色图综合在线观看| 欧美成人精品欧美一级黄| 蜜桃在线观看..| 制服丝袜香蕉在线| 各种免费的搞黄视频| 日日爽夜夜爽网站| 一级黄片播放器| 亚洲图色成人| 18+在线观看网站| 国产男女超爽视频在线观看| av播播在线观看一区| 国产熟女午夜一区二区三区 | 视频区图区小说| 97在线视频观看| 日日摸夜夜添夜夜爱| 天堂俺去俺来也www色官网| 国产欧美日韩精品一区二区| 国产成人精品久久久久久| 99热国产这里只有精品6| 日韩亚洲欧美综合| 欧美最新免费一区二区三区| 国产精品一区二区在线不卡| 免费少妇av软件| 街头女战士在线观看网站| 好男人视频免费观看在线| 国产在线男女| av又黄又爽大尺度在线免费看| 一级黄片播放器| 夫妻性生交免费视频一级片| 黑人猛操日本美女一级片| 高清毛片免费看| 男的添女的下面高潮视频| 久久久久久久精品精品| 国产精品国产av在线观看| 亚洲欧美成人综合另类久久久| 国产片特级美女逼逼视频| 久久精品国产亚洲av涩爱| 97超视频在线观看视频| 内地一区二区视频在线| 国产午夜精品久久久久久一区二区三区| 一区二区三区免费毛片| 不卡视频在线观看欧美| 各种免费的搞黄视频| 婷婷色av中文字幕| 人人妻人人爽人人添夜夜欢视频 | 国内揄拍国产精品人妻在线| 久久久久久久久大av| 日本vs欧美在线观看视频 | 大又大粗又爽又黄少妇毛片口| 久久av网站| 在线观看av片永久免费下载| 人体艺术视频欧美日本| 国产精品一区二区在线观看99| 一区二区三区精品91| 欧美日韩在线观看h| 夫妻午夜视频| 精品99又大又爽又粗少妇毛片| 久久久久久久久久久久大奶| 免费观看无遮挡的男女| 国产精品99久久久久久久久| 午夜久久久在线观看| 亚洲美女搞黄在线观看| 亚洲欧美精品专区久久| 人妻少妇偷人精品九色| 亚洲欧美一区二区三区黑人 | 亚洲第一区二区三区不卡| 久久久久久久精品精品| 国产精品不卡视频一区二区| 国产亚洲欧美精品永久| 观看免费一级毛片| 免费播放大片免费观看视频在线观看| 国产精品久久久久久精品古装| 国产成人精品无人区| 久久6这里有精品| 午夜视频国产福利| 日韩大片免费观看网站| 亚洲精品一区蜜桃| av免费观看日本| 亚洲人与动物交配视频| 午夜福利网站1000一区二区三区| 最近手机中文字幕大全| 看十八女毛片水多多多| 亚洲,一卡二卡三卡| 久久毛片免费看一区二区三区| 一级毛片 在线播放| 啦啦啦啦在线视频资源| 2022亚洲国产成人精品| 最近2019中文字幕mv第一页| 日韩欧美 国产精品| av卡一久久| 18+在线观看网站| 亚洲国产欧美在线一区| 丝袜喷水一区| 亚洲国产av新网站| 99热这里只有精品一区| 最新的欧美精品一区二区| 亚洲精品国产av成人精品| 欧美日本中文国产一区发布| 青春草亚洲视频在线观看| 亚洲第一av免费看| 亚洲精品乱久久久久久| 91成人精品电影| 3wmmmm亚洲av在线观看| 91精品国产九色| 国产黄色免费在线视频| 亚洲欧美中文字幕日韩二区| 久久人妻熟女aⅴ| 免费观看无遮挡的男女| 看免费成人av毛片| 最近中文字幕高清免费大全6| 久久毛片免费看一区二区三区| 国产欧美另类精品又又久久亚洲欧美| 久久精品国产鲁丝片午夜精品| 一区二区三区精品91| 亚洲婷婷狠狠爱综合网| 国产中年淑女户外野战色| 我要看黄色一级片免费的| a级一级毛片免费在线观看| 免费人妻精品一区二区三区视频| 国内精品宾馆在线| 国产精品国产三级国产专区5o| 日韩在线高清观看一区二区三区| 欧美精品一区二区免费开放| 王馨瑶露胸无遮挡在线观看| 极品少妇高潮喷水抽搐| 全区人妻精品视频| 久久综合国产亚洲精品| 亚洲三级黄色毛片| 最黄视频免费看| 久久人人爽av亚洲精品天堂| 18禁动态无遮挡网站| 成人亚洲精品一区在线观看| 日本vs欧美在线观看视频 | 亚洲,一卡二卡三卡| 亚洲一级一片aⅴ在线观看| 日韩免费高清中文字幕av| 亚洲欧美中文字幕日韩二区| 三级国产精品片| 国产精品蜜桃在线观看| 久久婷婷青草| 一个人看视频在线观看www免费| 亚洲第一区二区三区不卡| 国产成人精品婷婷| 日本欧美国产在线视频| 亚洲三级黄色毛片| 成年美女黄网站色视频大全免费 | 亚洲精品中文字幕在线视频 | 日产精品乱码卡一卡2卡三| 亚洲国产精品一区三区| 欧美日韩亚洲高清精品| 国产熟女欧美一区二区| 久久久久人妻精品一区果冻| 亚洲图色成人| 国产免费一级a男人的天堂| 国产精品欧美亚洲77777| 狠狠精品人妻久久久久久综合| 99精国产麻豆久久婷婷| 噜噜噜噜噜久久久久久91| 午夜精品国产一区二区电影| 亚洲国产日韩一区二区| 欧美精品国产亚洲| 欧美 日韩 精品 国产| 三级国产精品欧美在线观看| 久久青草综合色| 亚洲国产毛片av蜜桃av| 成人毛片60女人毛片免费| 777米奇影视久久| 乱码一卡2卡4卡精品| 国产一区二区在线观看av| 久久久久久人妻| 亚洲天堂av无毛| 麻豆精品久久久久久蜜桃| 亚洲av国产av综合av卡| 最近的中文字幕免费完整| 永久免费av网站大全| 老司机亚洲免费影院| 日本-黄色视频高清免费观看| 国产免费一区二区三区四区乱码| 久久女婷五月综合色啪小说| 少妇猛男粗大的猛烈进出视频| 国产欧美亚洲国产| av卡一久久| 色94色欧美一区二区| 毛片一级片免费看久久久久| 黄色日韩在线| 精品国产国语对白av| 国产白丝娇喘喷水9色精品| 啦啦啦中文免费视频观看日本| 伦精品一区二区三区| 亚洲av电影在线观看一区二区三区| 欧美变态另类bdsm刘玥| a级毛片免费高清观看在线播放| 午夜视频国产福利| 丰满迷人的少妇在线观看| 中国美白少妇内射xxxbb| 成人18禁高潮啪啪吃奶动态图 | 人妻制服诱惑在线中文字幕| 亚洲精品久久午夜乱码| 成年人免费黄色播放视频 | a级片在线免费高清观看视频| 中国国产av一级| 中国美白少妇内射xxxbb| 9色porny在线观看| 内射极品少妇av片p| 国产黄片视频在线免费观看| 成人亚洲精品一区在线观看| 亚洲va在线va天堂va国产| 夜夜看夜夜爽夜夜摸| 欧美日韩在线观看h| 男人狂女人下面高潮的视频| 日韩不卡一区二区三区视频在线| 精品亚洲成a人片在线观看| 99久久中文字幕三级久久日本| 亚洲国产精品国产精品| 亚洲一级一片aⅴ在线观看| 在线观看一区二区三区激情| 国产精品人妻久久久影院| 少妇精品久久久久久久| 又粗又硬又长又爽又黄的视频| 国产成人精品一,二区| 午夜激情久久久久久久| 欧美高清成人免费视频www| 一级二级三级毛片免费看| 日日摸夜夜添夜夜添av毛片| 天美传媒精品一区二区| 大香蕉久久网| 成人无遮挡网站| av.在线天堂| 国产精品久久久久成人av| 高清黄色对白视频在线免费看 | 国产一区二区在线观看日韩| 国产乱人偷精品视频| 97超视频在线观看视频| 久热这里只有精品99| 黑人猛操日本美女一级片| 欧美老熟妇乱子伦牲交| 日日爽夜夜爽网站| 九九在线视频观看精品| 老司机影院成人| 男女啪啪激烈高潮av片| 日本黄色片子视频| 亚洲欧美成人综合另类久久久| 中文字幕人妻熟人妻熟丝袜美| 丰满人妻一区二区三区视频av| 中文在线观看免费www的网站| av又黄又爽大尺度在线免费看| h视频一区二区三区| 国产高清国产精品国产三级| 国产探花极品一区二区| 日日摸夜夜添夜夜爱| 国产成人一区二区在线| 少妇丰满av| 亚洲天堂av无毛| 夜夜骑夜夜射夜夜干| 人人妻人人爽人人添夜夜欢视频 | 中文精品一卡2卡3卡4更新| 久久人人97超碰香蕉20202| 最新的欧美精品一区二区| 亚洲自偷自拍图片 自拍| 亚洲av日韩精品久久久久久密| 久久热在线av| 午夜福利一区二区在线看| 久久久国产一区二区| 国产在线一区二区三区精| 亚洲精品自拍成人| 精品高清国产在线一区| 淫妇啪啪啪对白视频 | 色94色欧美一区二区| 制服人妻中文乱码| 亚洲精品国产精品久久久不卡| 免费少妇av软件| 精品一区二区三区四区五区乱码| 91字幕亚洲| 日本一区二区免费在线视频| 国产无遮挡羞羞视频在线观看| 汤姆久久久久久久影院中文字幕| 亚洲欧美清纯卡通| 男女无遮挡免费网站观看| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美精品亚洲一区二区| 美女脱内裤让男人舔精品视频| 日本av手机在线免费观看| a级片在线免费高清观看视频| 夫妻午夜视频| 美女中出高潮动态图| 99热网站在线观看| 韩国精品一区二区三区| 欧美日韩一级在线毛片| 亚洲精品国产区一区二| 免费看十八禁软件| 午夜激情av网站| 满18在线观看网站| 我要看黄色一级片免费的| 午夜福利在线免费观看网站| 久久 成人 亚洲| 亚洲精品粉嫩美女一区| 日本av手机在线免费观看| 亚洲男人天堂网一区| 国产人伦9x9x在线观看| 久久天堂一区二区三区四区| 高清在线国产一区| 80岁老熟妇乱子伦牲交| 高清av免费在线| 精品少妇内射三级| 欧美日韩亚洲高清精品| 大香蕉久久成人网| 色婷婷久久久亚洲欧美| 亚洲国产欧美一区二区综合| 久久精品亚洲熟妇少妇任你| 狠狠婷婷综合久久久久久88av| 在线观看免费高清a一片| 十八禁网站免费在线| 日日摸夜夜添夜夜添小说| 亚洲中文av在线| 天堂8中文在线网| 国产欧美日韩精品亚洲av| 91老司机精品| 久久久水蜜桃国产精品网| 久久人人爽av亚洲精品天堂| 欧美变态另类bdsm刘玥| 国产高清videossex| 久久亚洲精品不卡| 99热网站在线观看| 每晚都被弄得嗷嗷叫到高潮| 超色免费av| 九色亚洲精品在线播放| 日韩视频一区二区在线观看| 99国产精品99久久久久| 美女高潮到喷水免费观看| 一级毛片女人18水好多| 亚洲精品美女久久久久99蜜臀| 人人妻人人添人人爽欧美一区卜| 亚洲久久久国产精品| 色视频在线一区二区三区| 男女床上黄色一级片免费看| 50天的宝宝边吃奶边哭怎么回事| 国产一卡二卡三卡精品| 亚洲国产精品一区三区| 亚洲,欧美精品.| 人人澡人人妻人| 国产在线视频一区二区| 黄色毛片三级朝国网站| 日韩熟女老妇一区二区性免费视频| 91精品国产国语对白视频| 亚洲av日韩精品久久久久久密| 91精品伊人久久大香线蕉| 久久国产精品大桥未久av| 国产在视频线精品| 亚洲七黄色美女视频| 国产精品久久久久久精品电影小说| 老鸭窝网址在线观看| 侵犯人妻中文字幕一二三四区| 老鸭窝网址在线观看| 在线观看www视频免费| 在线av久久热| 欧美激情 高清一区二区三区| 搡老熟女国产l中国老女人| 汤姆久久久久久久影院中文字幕| av在线app专区| 成人18禁高潮啪啪吃奶动态图| 最近中文字幕2019免费版| 国产精品一区二区精品视频观看| 亚洲欧美成人综合另类久久久| 久热这里只有精品99| 91精品伊人久久大香线蕉| 最近最新免费中文字幕在线| 久久久久久久国产电影| svipshipincom国产片| 国产三级黄色录像| 美女国产高潮福利片在线看| 国产国语露脸激情在线看| 亚洲精品久久成人aⅴ小说| 在线观看免费高清a一片| 午夜福利在线观看吧| 十八禁高潮呻吟视频| 亚洲av电影在线观看一区二区三区| 亚洲九九香蕉| 国产精品一二三区在线看| 日本一区二区免费在线视频| 韩国精品一区二区三区| 人人妻人人爽人人添夜夜欢视频| 美女主播在线视频| 久久久久久免费高清国产稀缺| 欧美成狂野欧美在线观看| 狠狠精品人妻久久久久久综合| 免费看十八禁软件| 亚洲欧美日韩另类电影网站| 亚洲精品久久久久久婷婷小说| 丰满饥渴人妻一区二区三| 如日韩欧美国产精品一区二区三区| 婷婷成人精品国产| 日本一区二区免费在线视频| 99国产极品粉嫩在线观看| 大码成人一级视频| 一区二区日韩欧美中文字幕| 国产一区二区 视频在线| 黄色视频不卡| 女人爽到高潮嗷嗷叫在线视频| 最近最新中文字幕大全免费视频| 少妇 在线观看| 成人黄色视频免费在线看| 黑人巨大精品欧美一区二区蜜桃| 亚洲成国产人片在线观看| 久久久精品免费免费高清| 人人澡人人妻人| 国产一区二区三区在线臀色熟女 | 久久久久国产精品人妻一区二区| 大香蕉久久成人网| 不卡一级毛片| 在线观看免费午夜福利视频| 无遮挡黄片免费观看| 欧美少妇被猛烈插入视频| 精品少妇黑人巨大在线播放| 亚洲中文日韩欧美视频| 成年美女黄网站色视频大全免费| 男男h啪啪无遮挡| 性色av一级| 亚洲精品久久久久久婷婷小说| 欧美亚洲日本最大视频资源| 欧美人与性动交α欧美精品济南到| 波多野结衣av一区二区av| 91大片在线观看| 久久精品人人爽人人爽视色| 在线观看舔阴道视频| 中文字幕人妻丝袜制服| 亚洲专区国产一区二区| 青草久久国产| 久久久久久久大尺度免费视频| 日本撒尿小便嘘嘘汇集6| 午夜影院在线不卡| 欧美精品av麻豆av| 捣出白浆h1v1| 日韩免费高清中文字幕av| 一级,二级,三级黄色视频| 亚洲国产欧美日韩在线播放| 久久久久久亚洲精品国产蜜桃av| 80岁老熟妇乱子伦牲交| 一区二区三区四区激情视频|