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

    天津地區(qū)霾天氣特征研究

    2019-02-13 06:09:08蔡子穎韓素芹邱曉濱吳彬貴王雪蓮
    中國(guó)環(huán)境科學(xué) 2019年1期
    關(guān)鍵詞:天津地區(qū)能見度中度

    蔡子穎,姚 青,韓素芹,邱曉濱,張 敏,吳彬貴,王雪蓮

    ?

    天津地區(qū)霾天氣特征研究

    蔡子穎1,2,3*,姚 青1,3,韓素芹3,邱曉濱3,張 敏1,吳彬貴3,王雪蓮3

    (1.天津市環(huán)境氣象中心,天津 300074;2.中國(guó)氣象局大氣化學(xué)重點(diǎn)開放實(shí)驗(yàn)室,北京 100081;3.天津市氣象科學(xué)研究所,天津 300074)

    基于2014~2017年天津地區(qū)PM2.5質(zhì)量濃度,能見度和相對(duì)濕度監(jiān)測(cè)數(shù)據(jù)開展霾天氣特征研究.結(jié)果表明:天津中度以上霾過(guò)程分為五類:高壓后部型,北部弱高壓型,低壓槽型,均壓場(chǎng)型和鋒前低壓型.在現(xiàn)行標(biāo)準(zhǔn)下,中度霾一般對(duì)應(yīng)重度污染天氣,重度霾對(duì)應(yīng)重度到嚴(yán)重污染天氣;五級(jí)重度污染天氣一般有中-重度霾發(fā)生,六級(jí)嚴(yán)重污染天氣有重度霾天氣發(fā)生.2013年“大氣污染防治行動(dòng)計(jì)劃”開展以后,天津PM2.5質(zhì)量濃度和霾日均顯著減少,2017年相比2013年霾日減少了55%,中度及其以上霾日由2013年的41d下降到2017年的20d,下降幅度超過(guò)50%.基于實(shí)況監(jiān)測(cè)的PM2.5質(zhì)量濃度,能見度和相對(duì)濕度,可以較好的構(gòu)建區(qū)域能見度計(jì)算方程.統(tǒng)計(jì)數(shù)據(jù)顯示,其估算的能見度和實(shí)況值相關(guān)系數(shù)為0.94,相對(duì)誤差為18.6%,非霾日辨識(shí)準(zhǔn)確率為85%,霾日辨識(shí)準(zhǔn)確率為95.6%,輕微霾辨識(shí)準(zhǔn)確率為83%,輕度霾辨識(shí)準(zhǔn)確率為78%,中度霾辨識(shí)準(zhǔn)確率為93%,重度及以上霾辨識(shí)準(zhǔn)確率為94%,對(duì)于判斷霾等級(jí),有較強(qiáng)的適用性.將該方程與空氣質(zhì)量模式結(jié)合開展霾等級(jí)預(yù)報(bào),2015~2017年24h預(yù)報(bào)產(chǎn)品檢驗(yàn)顯示:能見度預(yù)報(bào)值與實(shí)況值相關(guān)系數(shù)為0.75,預(yù)報(bào)均值13.9km,實(shí)況均值14.1km,相對(duì)誤差為29.6%, FAC2(預(yù)報(bào)值在實(shí)況值兩倍范圍內(nèi)百分比)為98.1%,霾日預(yù)報(bào)準(zhǔn)確率81.4%,霾日漏報(bào)率18.6%,霾日空?qǐng)?bào)率20.6%,如果容錯(cuò)1級(jí),輕微霾日預(yù)報(bào)準(zhǔn)確率為96%,輕度霾日預(yù)報(bào)準(zhǔn)確率為85%,中度及以上霾日預(yù)報(bào)準(zhǔn)確率為69%,可有效支撐天津霾等級(jí)預(yù)報(bào)的開展.

    霾等級(jí)預(yù)報(bào);重污染天氣;天津

    近年伴隨人民日益增長(zhǎng)的美好生活需要,霾受政府重視和公眾關(guān)注,成為研究的熱點(diǎn)問(wèn)題.21世紀(jì)初吳兌[1]基于珠江三角洲開展霾與霧區(qū)別和灰霾天氣預(yù)警的討論,并研究珠三角洲大氣灰霾導(dǎo)致能見度下降問(wèn)題[2],陳歡歡[3]和吳蒙[4]延續(xù)吳兌工作,認(rèn)為廣東省典型灰霾過(guò)程分布在每年的10月~次年4月,造成灰霾的天氣形勢(shì)主要是高壓變性出海.在長(zhǎng)江三角洲地區(qū)霾的影響明顯,其垂直結(jié)構(gòu)[5-6],光學(xué)特性[7],機(jī)理[8]被細(xì)致的研究.京津冀地區(qū)是我國(guó)霾最高發(fā)的地區(qū),受地形影響[9],導(dǎo)致山前暖區(qū)空氣流動(dòng)性較小形成氣流停滯區(qū),污染物和水汽容易聚集從而有利于霾形成.綜合近年研究,霾天氣發(fā)生,排放是內(nèi)因[10],氣象是外因[11].靜穩(wěn)天氣[12]、小風(fēng)速[13]、高相對(duì)濕度[14]、大氣層結(jié)穩(wěn)定[15]、混合層厚度低[16-18]和弱氣壓場(chǎng)[19]都是霾天氣發(fā)生的重要?dú)庀笳T因.能見度的高低是霾天氣判斷的重要標(biāo)識(shí),影響能見度的天氣現(xiàn)象有霧、霾、沙塵和降水等,霧天氣能見度的大小與霧滴液態(tài)含水量和霧滴尺度密切相關(guān)[20];霾天氣能見度與氣溶膠和相對(duì)濕度密切相關(guān)[21-22],其貢獻(xiàn)達(dá)到大氣消光的85%以上[23-24];降水天氣出現(xiàn)低能見度,主要來(lái)自雨滴的消光作用[25].在霾判斷時(shí),需要扣除霧、沙塵和降水等視程障礙天氣的影響[26].關(guān)于霾天氣的預(yù)報(bào),一般可以分為天氣學(xué)方法[27],統(tǒng)計(jì)學(xué)方法[28]和數(shù)值預(yù)報(bào)[29-31]方法.隨著近年數(shù)值模式和高性能計(jì)算機(jī)技術(shù)的不斷發(fā)展,越來(lái)越多的數(shù)值模式被用于霾和能見度預(yù)報(bào).基于數(shù)值模式開展霾天氣和能見度預(yù)報(bào),需要解決大氣消光系數(shù)計(jì)算的問(wèn)題,常見的方法可以分為3類,一類是基于質(zhì)量濃度或者體積濃度觀測(cè)數(shù)據(jù),構(gòu)建區(qū)域統(tǒng)計(jì)方程和數(shù)值模式相互銜接計(jì)算大氣消光系數(shù)[32-34],其優(yōu)勢(shì)是方程構(gòu)建相對(duì)簡(jiǎn)單,且往往能取得不錯(cuò)的效果,缺點(diǎn)是需要長(zhǎng)序列觀測(cè)數(shù)據(jù)建立方程,方程構(gòu)建區(qū)域特征明顯,移植其它地區(qū)需要重建方程,物理意義略差;二是基于化學(xué)組分計(jì)算大氣消光系數(shù),典型的如IMPROVE方程[35-37],其優(yōu)點(diǎn)是物理意義明確性強(qiáng),且對(duì)大氣消光的成因可以有效分解,缺點(diǎn)是成分?jǐn)?shù)據(jù)采樣分析成本高,空氣質(zhì)量模型化學(xué)組分模擬相比質(zhì)量濃度模擬偏差大;三是基于MIE散射和粒徑譜的消光系數(shù)計(jì)算[38-39],其優(yōu)點(diǎn)是物理意義明確,且無(wú)需觀測(cè)數(shù)據(jù)新建方程;缺點(diǎn)是模式粒徑譜與實(shí)況存在偏差,導(dǎo)致大氣消光系數(shù)計(jì)算存在誤差.

    相比傳統(tǒng)天氣業(yè)務(wù).環(huán)境氣象業(yè)務(wù)發(fā)展滯后.中國(guó)氣象局《應(yīng)對(duì)氣候變化保障生態(tài)文明建設(shè)規(guī)劃》明確要求,各省氣象局應(yīng)提高大氣污染防治氣象保障服務(wù)水平,到2020年建立集約化,0~10d無(wú)縫隙的環(huán)境氣象預(yù)報(bào)業(yè)務(wù),重污染天氣和霾預(yù)警時(shí)效提前至48h,霾、空氣質(zhì)量24h時(shí)效預(yù)報(bào)TS評(píng)分達(dá)到0.35和0.55,發(fā)布72h效逐3h霾預(yù)報(bào)產(chǎn)品及霾、能見度預(yù)報(bào)格點(diǎn)化產(chǎn)品.為有效支撐天津地區(qū)霾等級(jí)預(yù)報(bào)開展,本文在收集2014~2017年天津地區(qū)PM2.5質(zhì)量濃度、能見度和相對(duì)濕度監(jiān)測(cè)數(shù)據(jù)的基礎(chǔ)上,從監(jiān)測(cè)數(shù)據(jù)出發(fā),分析現(xiàn)行標(biāo)準(zhǔn)下天津霾天氣與重污染天氣(以AQI為標(biāo)準(zhǔn))的區(qū)別和聯(lián)系,以及兩者的時(shí)空演變規(guī)律特征.結(jié)合天氣學(xué)分析,構(gòu)建天津中度及以上霾天氣概念模型.在此背景下研究天津PM2.5質(zhì)量濃度、能見度和相對(duì)濕度三者的相互關(guān)系.構(gòu)建能見度計(jì)算方程并與天津氣象部門空氣質(zhì)量數(shù)值模式相互銜接,試驗(yàn)性提供霾等級(jí)預(yù)報(bào)客觀支撐產(chǎn)品,并基于2015~2017產(chǎn)品對(duì)方法進(jìn)行檢驗(yàn),最終以期提升天津氣象部門霾天氣能見度和霾等級(jí)預(yù)報(bào)能力,服務(wù)天津生態(tài)文明建設(shè).

    1 研究方法

    1.1 監(jiān)測(cè)數(shù)據(jù)

    圖1 2014~2016年中國(guó)氣象局天津大氣邊界層觀測(cè)站與天津市生態(tài)環(huán)境監(jiān)測(cè)中心發(fā)布PM2.5質(zhì)量濃度對(duì)比

    本文研究氣象數(shù)據(jù)來(lái)自天津市13個(gè)區(qū)級(jí)氣象站;空氣質(zhì)量數(shù)據(jù)來(lái)自生態(tài)環(huán)境保護(hù)部門和氣象部門監(jiān)測(cè),其中2014年1月1日~2017年12月31日為天津市生態(tài)環(huán)境監(jiān)測(cè)中心發(fā)布的逐小時(shí)PM2.5質(zhì)量濃度監(jiān)測(cè)數(shù)據(jù),監(jiān)測(cè)站點(diǎn)27個(gè),監(jiān)測(cè)方法為震蕩天平法,增加揮發(fā)有機(jī)物補(bǔ)償測(cè)量;2009年1月1日~2013年12月31日為中國(guó)氣象局天津大氣邊界層觀測(cè)站(北緯39°06¢,東經(jīng)117°10¢,海拔高度2.2m,臺(tái)站編號(hào):54517)逐小時(shí)監(jiān)測(cè)數(shù)據(jù),監(jiān)測(cè)方法為震蕩天平法,兩者數(shù)據(jù)通過(guò)2014~2016年3a同步觀測(cè)(圖1),進(jìn)行歸一化修正.本文使用數(shù)據(jù)經(jīng)過(guò)嚴(yán)格控制比對(duì).

    1.2 模式設(shè)置

    WRF/chem模式[40-42]是NCAR(national center for atmospheric research)和NOAA(national oceanic and atmospheric administration )聯(lián)合一些大學(xué)和研發(fā)機(jī)構(gòu)開發(fā)的中尺度在線大氣化學(xué)模式,模式考慮大氣污染的化學(xué)過(guò)程,平流輸送,湍流擴(kuò)散,干濕沉降過(guò)程,在全球空氣質(zhì)量預(yù)報(bào)和模擬中有廣泛的運(yùn)用.本文采用WRF/ chem3.8.1版本,氣相化學(xué)過(guò)程采用CBMZ機(jī)制,氣溶膠過(guò)程采用MOSAIC模型,主要物理過(guò)程設(shè)置如下:積云對(duì)流方案采用Grell-3D,微物理過(guò)程采用WSM5,長(zhǎng)波輻射方案和短波輻射方案均采用RRTMG,邊界層方案使用YSU方案,模式水平分辨率15km,水平網(wǎng)格121×121,中心經(jīng)緯度為39N°,117E°,垂直方向分為41層.模式的人為排放源清單使用清華大學(xué)MEIC (Multi-resolution Emission Inventory for China)2012,分辨率0.25°×0.25°,在天津地區(qū)使用27個(gè)空氣質(zhì)量監(jiān)測(cè)站實(shí)況數(shù)據(jù)和相關(guān)排放源統(tǒng)計(jì)信息進(jìn)行時(shí)空的細(xì)化,氣象初始場(chǎng)和背景場(chǎng)均使用NECP的FNL全球1°×1°數(shù)據(jù),模擬時(shí)間為2015年1月1日~2017年12月31日,模擬采用24h滾動(dòng)計(jì)算,每24h重新使用一次FNL氣象初始場(chǎng),而污染初始場(chǎng)則為上一次的模擬值.

    2 結(jié)果與討論

    2.1 霾天氣與重污染天氣的聯(lián)系和區(qū)別

    2010年氣象行業(yè)標(biāo)準(zhǔn)《霾的觀測(cè)和預(yù)報(bào)等級(jí)》[43]首次引入了大氣成分指標(biāo)作為霾等級(jí)的判識(shí)標(biāo)準(zhǔn),認(rèn)為當(dāng)能見度低于10km,相對(duì)濕度低于80%時(shí)是霾,相對(duì)濕度為80%~95%時(shí)PM2.5質(zhì)量濃度大于75μg/m3為霾.且依據(jù)能見度()對(duì)霾等級(jí)進(jìn)行劃分:輕微霾(5km£<10km),輕度霾(3km£<5km),中度霾(2km£<3km),重度霾(<2km).2013年中國(guó)氣象局在充分考慮能見度、相對(duì)濕度和PM2.5質(zhì)量濃度的基礎(chǔ)上制定霾預(yù)警(暫行)標(biāo)準(zhǔn),明確霾等級(jí)劃分方法(表1),并以此為依托發(fā)布霾預(yù)警.

    依據(jù)重污染天氣應(yīng)急預(yù)案定義,重污染天氣為空氣質(zhì)量指數(shù)(AQI)級(jí)別達(dá)到五級(jí)(重度污染)及以上污染程度的天氣.2013~2017年監(jiān)測(cè)數(shù)據(jù)顯示,天津地區(qū)AQI大于200的重污染天氣共計(jì)169d,其中首要污染物為PM2.5的163d,為PM10的6d(5d受到顯著的沙塵天氣影響),未出現(xiàn)O3和其它污染物為首要污染物的重污染天氣,即在天津地區(qū)重污染天氣一般為PM2.5為首要污染物的五級(jí)及以上重污染天氣,即PM2.5質(zhì)量濃度大于150mg/m3.

    表1 中國(guó)氣象局霾預(yù)警標(biāo)準(zhǔn)(暫行)

    綜上所述,基于現(xiàn)行業(yè)務(wù)規(guī)范和天津大氣環(huán)境現(xiàn)狀,天津霾等級(jí)的判斷依賴于能見度、相對(duì)濕度和PM2.5質(zhì)量濃度.重污染天氣的判斷主要為PM2.5質(zhì)量濃度,當(dāng)PM2.5質(zhì)量濃度大于150mg/m3即認(rèn)為出現(xiàn)重污染天氣.研究顯示,在非霧天氣,大氣能見度的下降主要為細(xì)顆粒物的吸收和散射,其中PM2.5在其中發(fā)揮最為主要的作用,在濕度不變的情況,PM2.5質(zhì)量濃度與能見度基本呈現(xiàn)冪指數(shù)關(guān)系[33],確定一個(gè)地區(qū)PM2.5體積濃度(如果體積譜確定,體積濃度和質(zhì)量濃度可進(jìn)行換算)和相對(duì)濕度,可以計(jì)算出該區(qū)域的大氣消光系數(shù)和能見度.所以確定一個(gè)區(qū)域PM2.5質(zhì)量濃度和相對(duì)濕度,可對(duì)霾的等級(jí)進(jìn)行粗略辨識(shí).基于上述討論,依托2014~ 2017年天津地區(qū)實(shí)測(cè)的資料,分析天津霾天氣和重污染天氣的區(qū)別和聯(lián)系.

    由表2~3統(tǒng)計(jì)顯示,基于上述標(biāo)準(zhǔn)判斷,天津地區(qū)輕度霾天氣,一般對(duì)應(yīng)輕度-中度污染,其平均能見度4.15km,PM2.5質(zhì)量濃度112mg/m3;中度霾天氣平均能見度3.40km,PM2.5質(zhì)量濃度177mg/m3,90%為重度污染天氣,10%為中度污染天氣;重度及以上霾天氣,一般對(duì)應(yīng)重度到嚴(yán)重污染天氣,其概率為重污染天氣51%,嚴(yán)重污染天氣49%.當(dāng)空氣質(zhì)量為優(yōu),能見度均在10km以上,沒有霾天氣發(fā)生;空氣質(zhì)量為良,有11%的天氣能見度介于5~10km,出現(xiàn)輕微霾天氣;當(dāng)空氣質(zhì)量為輕度污染,有35%天氣由于相對(duì)濕度較低,氣溶膠消光能力弱,能見度在10km以上,65%天氣出現(xiàn)輕微到輕度霾天氣;當(dāng)空氣質(zhì)量達(dá)到中度污染等級(jí),平均能見度降至6.42km,90%的情況有霾發(fā)生,其中輕度霾占26%,有7%的天氣由于相對(duì)濕度較高,氣溶膠吸濕增長(zhǎng)明顯,能見度低于3km,出現(xiàn)中度霾天氣;當(dāng)空氣質(zhì)量達(dá)到重度污染水平,66%的天氣將出現(xiàn)中度霾,20%的天氣將出現(xiàn)重度霾,但也有部分天氣,由于相對(duì)濕度非常低(空氣干燥),能見度大于5km,甚至大于10km,這樣的天氣給公眾的感官較好,但大氣中高負(fù)載的氣溶膠對(duì)人體也將產(chǎn)生嚴(yán)重傷害,其占比為13%.當(dāng)空氣質(zhì)量達(dá)到嚴(yán)重污染程度,平均能見度將降至1.50km,在所有研究樣本中僅有1d,能見度在5km以上,大部分個(gè)例能見度介于1~2km之間,一般對(duì)應(yīng)重度霾天氣.綜合而言:天津地區(qū)霾天氣和重污染天氣雖然定義不同,霾天氣側(cè)重于能見度和PM2.5質(zhì)量濃度的綜合判斷,重污染天氣主要依賴于PM2.5質(zhì)量濃度判斷,但由于霾天氣能見度的高低主要為細(xì)顆粒物的散射和吸收作用決定,所以霾等級(jí)與空氣質(zhì)量等級(jí)存在較為固定的關(guān)系.一般而言,中度霾對(duì)應(yīng)重度污染天氣(90%),重度霾對(duì)應(yīng)重度到嚴(yán)重污染天氣(95%);反之,五級(jí)重度污染天氣一般有中-重度霾發(fā)生(87%),六級(jí)嚴(yán)重污染天氣有重度霾發(fā)生(95%).重污染天氣(包含五級(jí)、六級(jí))88%有中度-重度霾發(fā)生,僅有12%由于大氣過(guò)于干燥,能見度大于5km.

    表2 不同等級(jí)霾對(duì)應(yīng)的氣象要素和空氣質(zhì)量等級(jí)

    表3 不同等級(jí)空氣質(zhì)量對(duì)應(yīng)的氣象要素和霾等級(jí)

    2.2 霾天氣和重污染天氣的統(tǒng)計(jì)特征

    統(tǒng)計(jì)2014~2017年天津地區(qū)霾天氣和重污染天氣特征.圖2a顯示由于氣象條件以及采暖季和非采暖季排放源差異,PM2.5質(zhì)量濃度呈現(xiàn)單峰結(jié)構(gòu),每年12月~次年1月為峰值,5月~9月為谷值.對(duì)應(yīng)PM2.5質(zhì)量濃度單峰結(jié)構(gòu),重污染天氣也呈現(xiàn)單峰結(jié)構(gòu),每年10月~次年3月多有重污染天氣發(fā)生,出現(xiàn)概率15.5%,其余時(shí)間重污染天氣發(fā)生概率0.7%.霾天氣與重污染天氣略有區(qū)別,天津能見度峰值出現(xiàn)在5~6月,7月雨季開始,相對(duì)濕度明顯增加,雖然PM2.5質(zhì)量濃度較低,但能見度小于10km天數(shù)明顯增加,霾出現(xiàn)概率34%,一般以輕微霾為主,未出現(xiàn)中度及以上霾天氣.對(duì)應(yīng)PM2.5質(zhì)量濃度分布,10月~次年3月,也是霾高發(fā)季節(jié),霾發(fā)生概率41%,尤其中度及以上霾,96%發(fā)生在每年10月~次年3月(圖2b).綜合而言,天津地區(qū)霾天氣全年均有發(fā)生,但中度及以上霾呈現(xiàn)明顯的季節(jié)差異,96%發(fā)生在每年10月~次年3月,與重污染天氣月季分布規(guī)律基本一致.

    圖2 2014~2017年天津地區(qū)霾天氣統(tǒng)計(jì)特征

    a 霾月變化,b中度以上霾發(fā)生概率,c PM2.5質(zhì)量濃度年際變化,d 霾日數(shù)年際變化,e 日變化,f 發(fā)生概率日變化

    年際變化分析,天津地區(qū)PM2.5質(zhì)量濃度2009~ 2012年呈現(xiàn)平穩(wěn)變化,維持在86mg/m3左右,受極端不利氣象條件影響,2013年P(guān)M2.5質(zhì)量濃度躍遷至96mg/m3,《大氣污染防治行動(dòng)計(jì)劃》頒布以后PM2.5質(zhì)量濃度呈現(xiàn)逐年下降的趨勢(shì)(圖2c),相比峰值濃度,2017年P(guān)M2.5質(zhì)量濃度下降了35.4%,相比2009~ 2012年的穩(wěn)定濃度,PM2.5質(zhì)量濃度下降28%.隨著大氣中PM2.5質(zhì)量濃度的減少,霾日呈現(xiàn)顯著下降趨勢(shì),2017年相比2013年霾日下降了55%,中度及其以上霾日由2013年的41d下降到2017年的20d,下降幅度超過(guò)50%,重污染天氣與中度及以上霾日變化趨勢(shì)一致,在2009~2013年保持45~50d范圍變化,經(jīng)過(guò)5a的治理,2017年僅為23d,重污染日顯著減少(圖2d).

    基于小時(shí)數(shù)據(jù)分析霾時(shí)和重污染時(shí)日變化.天津地區(qū)PM2.5質(zhì)量濃度呈現(xiàn)單峰型分布(圖2e),每日15h~17h為谷值,8h~10h為峰值,傍晚18h雖然PM2.5質(zhì)量濃度有所增加,但晚高峰并不明顯,18h~22h有持續(xù)增加,然后保持平穩(wěn),直到早晨再次出現(xiàn)峰值.能見度日變化規(guī)律與PM2.5質(zhì)量濃度相反,但波動(dòng)幅度明顯更大,其差異主要因?yàn)槟芤姸仁芟鄬?duì)濕度日變化影響.下半夜到清晨太陽(yáng)升起前相對(duì)濕度出現(xiàn)每日高值,此時(shí)能見度也為全天最低,日出后雖然PM2.5質(zhì)量濃度有小幅增加,但伴隨氣溫的增加,導(dǎo)致濕度降低,能見度呈現(xiàn)增加趨勢(shì).基于上述規(guī)律,天津地區(qū)霾最易發(fā)生在下半夜到上午,峰值出現(xiàn)在6h~9h,有接近50%的概率出現(xiàn)霾,且出現(xiàn)中度及以上霾概率達(dá)到13%(圖2f),午后由于大氣垂直擴(kuò)散條件改善以及相對(duì)濕度下降,霾發(fā)生概率大幅度減小,其發(fā)生概率約在25%~30%.以小時(shí)計(jì)算,2017年天津地區(qū)霾約為2164h,其中中度及以上霾發(fā)生590h,相比2014年分別減少1944h和657h,但霾時(shí)仍占所有時(shí)間的25%,霾天氣依舊是影響天津人民群眾生活的重要問(wèn)題.與霾天氣對(duì)應(yīng)的是重污染天氣,重污染天氣日變化分布相對(duì)平緩,9h~11h發(fā)生的概率略高,為11.6%,每日20h~次日12h發(fā)生概率為10%~11%,午后受垂直擴(kuò)散條件改善影響,在8%~10%之間.對(duì)比霾日間變化,重污染受日間變化影響更小,持續(xù)性更強(qiáng).

    2.3 霾天氣和重污染天氣的天氣概念模型

    2014~2017年發(fā)生中度及以上霾為111d,重污染天氣為118d,重合天數(shù)為104d.出現(xiàn)中度及以上霾天氣,未出現(xiàn)重污染天氣7d(出現(xiàn)中度污染).地面以高壓后(東風(fēng)增濕)或者鋒前低壓(弱北風(fēng)降溫增濕)為主,平均相對(duì)濕度為82%,由于相對(duì)濕度較高,雖然未出現(xiàn)重度污染,但能見度較低(平均2.5km),達(dá)到中度霾等級(jí);出現(xiàn)重污染天氣,未出現(xiàn)中度霾天氣為14d,其中850hPa呈現(xiàn)偏西或者西北氣流為10d,平均相對(duì)濕度為50%,低相對(duì)濕度導(dǎo)致能見度保持在5~10km,且其中50%的過(guò)程地面已經(jīng)轉(zhuǎn)為偏西或者西北風(fēng),西北風(fēng)、西風(fēng)使得天津地區(qū)相對(duì)濕度快速下降,但北部、西北部輸送的污染物仍然使得天津有較高的細(xì)顆粒物質(zhì)量濃度,空氣質(zhì)量改善時(shí)間明顯滯后于能見度改善時(shí)間.

    500hPa環(huán)流分析,當(dāng)高空環(huán)流為平直西風(fēng)氣流、弱西北氣流和槽前西南氣流時(shí)易于霾天氣發(fā)生.具體而言,當(dāng)高空環(huán)流為平直西風(fēng)氣流時(shí),西風(fēng)帶無(wú)明顯波動(dòng),無(wú)明顯天氣系統(tǒng)影響華北地區(qū),氣壓場(chǎng)較弱,易于霾天氣出現(xiàn),占所有過(guò)程的34%;當(dāng)高空處于弱西北氣流控制時(shí),華北地區(qū)受脊前或者槽后弱西北氣流影響,地面以高壓后、弱高壓為主,也有部分前傾槽,地面位于低壓后部或者鋒前低壓區(qū),易于霾天氣出現(xiàn),占所有過(guò)程34%;當(dāng)高空為槽前西南氣流時(shí),華北地區(qū)位于高空槽前,地面多為鋒前低壓區(qū)或均壓場(chǎng)控制,輻合風(fēng)場(chǎng)有利于污染物的累積,占所有過(guò)程32%.

    850hPa流場(chǎng)分析,西南氣流是霾天氣最重要的影響氣流,占61%.當(dāng)850hPa為西南氣流時(shí),西南氣流將攜帶河北中南部的污染物輸送到天津,且如果華北地區(qū)受明顯的暖舌控制,暖舌的影響利于逆溫形成,或者西南氣流會(huì)帶來(lái)暖平流,暖平流利于逆溫形成,引起污染物累積.除西南氣流以外,850hPa呈現(xiàn)偏西氣流占13.5%,其多對(duì)應(yīng)地形槽或者華北小低壓天氣,南北向太行山與西風(fēng)下沉,易形成輻合低壓污染.850hPa呈西北流占23%,地面多對(duì)應(yīng)鋒前低壓區(qū),雖然已經(jīng)轉(zhuǎn)為偏北風(fēng)影響,但風(fēng)場(chǎng)較弱,北部輸送和鋒面逆溫導(dǎo)致霾天氣加劇.比較特殊的是前傾槽過(guò)程,高空850hPa為西北氣流,呈現(xiàn)下沉趨勢(shì),地面仍然位于低壓槽前,呈現(xiàn)上升氣流,下沉氣流抑制污染物的垂直擴(kuò)散,而地面弱上升,意味著仍處于低壓輻合區(qū),水平和垂直擴(kuò)散條件的雙重不利,導(dǎo)致污染過(guò)程加劇.

    結(jié)合500hPa環(huán)流分析和850hPa流場(chǎng)分析,基于地面天氣形勢(shì),可將天津中度及以上霾天氣劃分為5種類型.分別為均壓場(chǎng)型、弱高壓型、鋒前低壓型、高壓后型和低壓槽型.高壓類型2類,分別為高壓后和弱高壓,其中高壓后部型是華北地區(qū)較為常見的一種污染天氣形勢(shì)(占比27%),由于河北中南部相比京津冀中北部大氣污染物排放量和氣象條件均有利污染物積累,當(dāng)高壓后-低壓前地面形勢(shì)出現(xiàn)時(shí),天津地區(qū)一般為暖平流控制,易于逆溫出現(xiàn),西南輸送和本地垂直擴(kuò)散條件不利,雙重疊加有利于霾天氣發(fā)生.弱高壓(占比9.9%),主要為北部弱高壓型,海平面氣壓場(chǎng)呈現(xiàn)北高南低格局,前期污染積累以及冷空氣的強(qiáng)度對(duì)于弱高壓型污染是否出現(xiàn)中度及以上霾天氣至關(guān)重要.當(dāng)前期冷空氣較弱時(shí),對(duì)區(qū)域性污染清除不徹底,冷空氣減弱后,易于大氣污染回流或者上游輸送滯留,導(dǎo)致霾天氣出現(xiàn),此類天氣需要注意弱高壓天氣下沉氣流對(duì)大氣污染物垂直擴(kuò)散的影響.低壓類型有3類,占比最多的為鋒前低壓(占比28.8%),該天氣一般為污染過(guò)程的最后階段,也是污染過(guò)程峰值階段,未來(lái)有冷鋒過(guò)境,目前處于低壓系統(tǒng)或弱氣壓場(chǎng)內(nèi).在此階段大氣污染受前期積累和上游輸送的共同的影響,呈現(xiàn)較高的峰值濃度.其次為均壓場(chǎng)(占比22%),弱風(fēng)場(chǎng)是此類污染天氣的最大特征,當(dāng)濕度較高時(shí),此類天氣也易出現(xiàn)霧霾交替產(chǎn)生的現(xiàn)象.低壓槽也是華北地區(qū)典型的低壓污染類型(占比11.7%),其屬于輻合風(fēng)場(chǎng)污染類型.由于我國(guó)處于西風(fēng)帶,太行山又呈現(xiàn)南北向,氣流過(guò)山后,氣柱伸長(zhǎng),空氣發(fā)生輻合,氣旋性渦度增加,下沉氣流絕熱增溫,在對(duì)流層低層產(chǎn)生暖溫度脊,使低層減壓,在華北平原形成低槽(華北地形槽).該地形槽常為地面氣壓場(chǎng)的暖性低槽,有時(shí)在地形槽內(nèi)出現(xiàn)地形低壓,強(qiáng)度弱、不發(fā)展.受其影響在華北平原常有低壓區(qū)輻合區(qū)的存在,當(dāng)?shù)蛪狠椇蠀^(qū)閉合時(shí),我們稱為華北小低壓,當(dāng)?shù)蛪狠椇蠀^(qū)不閉合時(shí),稱為地形槽.地形槽與低壓系統(tǒng)相互融合,使得華北平原地區(qū)風(fēng)場(chǎng)輻射,污染滯留,有利于霾天氣的出現(xiàn).

    圖3 中度及以上霾天氣概念模型

    根據(jù)上述分析,構(gòu)建天津地區(qū)霾天氣概念模型,I為高壓后型;II為北部弱高壓型;III低壓槽型;IV均壓場(chǎng)型;V鋒前低壓型.具體如下(圖3):

    2.4 PM2.5、相對(duì)濕度和能見度的關(guān)系

    大氣低能見度的形成原因復(fù)雜,且存在明顯的地域性和季節(jié)性,其對(duì)應(yīng)的天氣現(xiàn)象有霧、霾、沙塵和降水等.霾天氣能見度與氣溶膠、相對(duì)濕度密切相關(guān),其貢獻(xiàn)達(dá)到大氣消光的85%以上,氣溶膠質(zhì)量濃度、粒徑譜分布、化學(xué)組分、吸濕特性和黑碳及其混合狀態(tài)等性質(zhì)均對(duì)其產(chǎn)生影響,在PM2.5中硫酸鹽、硝酸鹽和含碳?xì)馊苣z是引起能加度降低的主要成分.如果固定相對(duì)濕度,清晰顯示PM2.5質(zhì)量濃度與能見度呈現(xiàn)冪指數(shù)關(guān)系(表4),參考馬楠等[33]和蔡子穎等[32]在天津武清和城區(qū)的工作,基于2014~ 2016年P(guān)M2.5質(zhì)量、相對(duì)濕度和能見度建立擬合式,共計(jì)25828h,以2017年數(shù)據(jù)作為檢驗(yàn).擬合式1,其中為大氣消光系數(shù),可通過(guò)3912/計(jì)算大氣能見度(最高為30km,超過(guò)30km取30km),值為7.63,值為0.9,值為0.479,相比使用2012年城區(qū)觀測(cè)擬合的數(shù)據(jù),本研究擬合數(shù)據(jù)采用天津全市樣本,且2014年后環(huán)保部門使用RP1405F測(cè)量PM2.5質(zhì)量,相比原有的氣象部門RP1400,其揮發(fā)性有機(jī)物的補(bǔ)償使得監(jiān)測(cè)PM2.5質(zhì)量濃度略高于以前觀測(cè).基于此、和略有調(diào)整(原擬合式[21]中為6.72,為0.89,為0.838).基于2017年全年樣本進(jìn)行檢驗(yàn),模擬能見度和實(shí)況能見度相關(guān)系數(shù)0.94(圖4),相對(duì)誤差18.6%,比使用原有公式模擬2017年能見度略有改進(jìn)(相關(guān)系數(shù)0.93,相對(duì)誤差22%).在不考慮PM2.5質(zhì)量和相對(duì)濕度預(yù)報(bào)誤差的情況下,使用式1計(jì)算霾等級(jí)(2014~2017年),霾日辨識(shí)準(zhǔn)確率為85%,霾日辨識(shí)準(zhǔn)確為95.6%,樣本500d;輕微霾辨識(shí)準(zhǔn)確率為83%,樣本338;輕度霾辨識(shí)準(zhǔn)確率為78%,樣本55d;中度霾辨識(shí)準(zhǔn)確率93%,樣本72d;重度及以上霾辨識(shí)準(zhǔn)確率為94%,樣本34d.假定誤差一個(gè)等級(jí)辨識(shí)記為正確,則各等級(jí)平均辨識(shí)準(zhǔn)確率為99.5%.由此表明該式對(duì)霾等級(jí)辨識(shí)具有較高的準(zhǔn)確性.同時(shí)也說(shuō)明在確定PM2.5質(zhì)量濃度和相對(duì)濕度的基礎(chǔ)上,依據(jù)式(1)計(jì)算能見度日值雖然有接近20%的誤差,但對(duì)于判斷霾等級(jí),有較強(qiáng)的適用性.尤其在中度及以上霾天氣,其等級(jí)辨識(shí)準(zhǔn)確率可達(dá)到93%以上.

    =′PM2.5(1-RH)-c′RH(1)

    表4 不同相對(duì)濕度段天津能見度(y)與PM2.5質(zhì)量濃度(x)的定量關(guān)系式

    圖4 大氣水平能見度模擬值和實(shí)測(cè)值對(duì)比

    2.5 霾天氣數(shù)值預(yù)報(bào)

    2.4節(jié)從理論確定了使用式(1)可較好的進(jìn)行霾天氣數(shù)值預(yù)報(bào).基于WRF/chem.模型構(gòu)建霾天氣數(shù)值模式,收集2015~2017年的24h預(yù)報(bào)產(chǎn)品進(jìn)行檢驗(yàn).從統(tǒng)計(jì)數(shù)據(jù)反映,基于式(1)和WRF/chem.模型可以較好的開展天津霾天氣的能見度模擬和預(yù)報(bào),其能見度預(yù)報(bào)值與實(shí)況值相關(guān)系數(shù)0.75,預(yù)報(bào)均值13.9km,實(shí)況均值14.1km,相對(duì)誤差29.6%,FAC2(預(yù)報(bào)值在實(shí)況值兩倍范圍內(nèi)百分比)為98.1%.預(yù)報(bào)值與實(shí)況出現(xiàn)誤差最主要的影響因素為PM2.5質(zhì)量濃度預(yù)報(bào)和實(shí)況的偏差,2015~2017年P(guān)M2.5質(zhì)量濃度預(yù)報(bào)值與實(shí)況值相關(guān)系數(shù)0.73,預(yù)報(bào)均值68μg/m3,實(shí)況均值69μg/m3,相對(duì)誤差33.5%,FAC2值為86%,相對(duì)濕度相比PM2.5質(zhì)量濃度預(yù)報(bào),無(wú)論相對(duì)誤差,還是相關(guān)系數(shù)均表現(xiàn)的更為優(yōu)異,其預(yù)報(bào)值與實(shí)況值相關(guān)系數(shù)為0.89,預(yù)報(bào)均值為53,實(shí)況均值為57,相對(duì)誤差為12.7%,FAC2值為98.1%.基于上述方法,開展霾天氣預(yù)報(bào)和模擬,霾日預(yù)報(bào)準(zhǔn)確率為81.4%,霾日漏報(bào)率為18.6%,霾日空?qǐng)?bào)率20.6%,可有效支撐霾預(yù)報(bào)的開展.對(duì)于精細(xì)化的霾等級(jí)預(yù)報(bào)開展,表5顯示大部分情況下霾等級(jí)預(yù)報(bào)的準(zhǔn)確率在50%以上,如果容錯(cuò)1級(jí),輕微霾日預(yù)報(bào)準(zhǔn)確率為96%,輕度霾日預(yù)報(bào)準(zhǔn)確率為85%,中度及以上霾日預(yù)報(bào)準(zhǔn)確率為69%,相比基于PM2.5的空氣質(zhì)量等級(jí)預(yù)報(bào),其等級(jí)預(yù)報(bào)準(zhǔn)確率35%,容錯(cuò)1級(jí)等級(jí)預(yù)報(bào)準(zhǔn)確率86%,基于式1和數(shù)值模式的天津霾等級(jí)預(yù)報(bào)準(zhǔn)確率與基于數(shù)值模式的空氣質(zhì)量預(yù)報(bào)準(zhǔn)確率大體相當(dāng),后期需要通過(guò)改進(jìn)數(shù)值模式PM2.5質(zhì)量濃度預(yù)報(bào)能力,提升天津地區(qū)霾等級(jí)預(yù)報(bào)水平.

    表5 基于數(shù)值模式天津地區(qū)霾等級(jí)預(yù)報(bào)準(zhǔn)確率(2015~2017年)

    3 結(jié)論

    3.1 基于2014~2017年天津地區(qū)實(shí)測(cè)的資料,分析霾天氣和重污染天氣的聯(lián)系.一般而言,中度霾對(duì)應(yīng)重度污染天氣,重度霾對(duì)應(yīng)重度到嚴(yán)重污染天氣;反之,五級(jí)重度污染天氣一般有中-重度霾發(fā)生,六級(jí)嚴(yán)重污染天氣有重度霾發(fā)生.

    3.2 天津地區(qū)霾天氣全年均有發(fā)生,中度及以上霾呈現(xiàn)明顯季節(jié)差異,96%發(fā)生在每年10月~次年3月,與重污染天氣月季分布規(guī)律基本一致.隨著大氣中PM2.5質(zhì)量濃度的減少,霾日呈現(xiàn)顯著下降趨勢(shì),2017年相比2013年霾日下降了55%,中度及其以上霾日由2013年的41d下降到2017年的20d,下降幅度超過(guò)50%.分析2013~2017年所有中度及以上霾過(guò)程,天津地區(qū)霾天氣可以劃分為5類,分別為高壓后部型、北部弱高壓型、低壓槽型、均壓場(chǎng)型和鋒前低壓型.

    3.3 PM2.5質(zhì)量濃度與能見度呈現(xiàn)冪指數(shù)關(guān)系,基于實(shí)況監(jiān)測(cè)可清晰描繪天津PM2.5質(zhì)量濃度、相對(duì)濕度和能見度三者關(guān)系,并建立能見度計(jì)算方程.實(shí)況監(jiān)測(cè)數(shù)據(jù)反算,該方程計(jì)算的能見度和實(shí)況能見度相關(guān)系數(shù)0.94,相對(duì)誤差18.6%,非霾日辨識(shí)準(zhǔn)確率85%,霾日辨識(shí)準(zhǔn)確95.6%,輕微霾辨識(shí)準(zhǔn)確率為83%,輕度霾辨識(shí)準(zhǔn)確率為78%,中度霾辨識(shí)準(zhǔn)確率93%,重度及以上霾辨識(shí)準(zhǔn)確率94%.假定誤差一個(gè)等級(jí)辨識(shí)記為正確,則各等級(jí)平均辨識(shí)準(zhǔn)確率為99.5%,對(duì)于判斷霾等級(jí),有較強(qiáng)的適用性.

    3.4 基于本地化能見度計(jì)算方程,結(jié)合WRF/chem數(shù)值模式開展天津地區(qū)能見度和霾等級(jí)預(yù)報(bào).2015~ 2017年24h預(yù)報(bào)產(chǎn)品檢驗(yàn)顯示:能見度預(yù)報(bào)值與實(shí)況值相關(guān)系數(shù)為0.75,預(yù)報(bào)均值為13.9km,實(shí)況均值為14.1km,相對(duì)誤差為29.6%,FAC2為98.1%;霾日預(yù)報(bào)準(zhǔn)確率為81.4%,霾日漏報(bào)率為18.6%,霾日空?qǐng)?bào)率為20.6%,如果容錯(cuò)1級(jí),輕微霾日預(yù)報(bào)準(zhǔn)確率為96%,輕度霾日預(yù)報(bào)準(zhǔn)確率為85%,中度及以上霾日預(yù)報(bào)準(zhǔn)確率為69%,其準(zhǔn)確率與基于數(shù)值模式的空氣質(zhì)量預(yù)報(bào)準(zhǔn)確率大體相當(dāng),可有效支撐霾等級(jí)預(yù)報(bào)的開展,在后期研究中,也需要通過(guò)改進(jìn)數(shù)值模式PM2.5質(zhì)量濃度預(yù)報(bào)能力,提升天津霾等級(jí)預(yù)報(bào)水平.

    [1] 吳 兌.關(guān)于霾與霧的區(qū)別和灰霾天氣預(yù)警的討論[J]. 氣象, 2005, 31(4):3-7. Wu D. A Discussion on Difference between Haze and Fog and Warning of Ash Haze Weather [J]. Meteorological monthly. 2005, 31(4):3-7.

    [2] 吳 兌,畢雪巖,鄧雪嬌,等.珠江三角洲大氣灰霾導(dǎo)致能見度下降問(wèn)題研究[J]. 氣象學(xué)報(bào), 2006,64(4):510-517.Wu D, Bi X, Deng X J, et al. Effect of atmospheric haze on the deterioration of visibility over the pear river detla [J]. Acta Meteorological sinica, 2006,64(4):510-517.

    [3] 陳歡歡,吳 兌,譚浩波,等.珠江三角洲2001~2008年灰霾天氣過(guò)程特征分析[J]. 熱帶氣象學(xué)報(bào), 2010,26(2):147-155.Chen H H, Wu D, Tan H B, et al. Study on the character of haze weather process from the year 2001 to 2008 over the pearl river delta [J]. Journal of tropical meteorology, 2010,26(2):147-155.

    [4] 吳 蒙,范紹佳,吳 兌,等.廣州地區(qū)灰霾與清潔天氣變化特征及影響因素分析[J]. 中國(guó)環(huán)境科學(xué), 2012,32(8):1409-1415.Wu M, Fan S J, Wu D, et al. Characteristics and influence factors study of haze and cleaning weather over Guangzhou Area [J]. China Environmental Science. 2012,32(8):1409-1415.

    [5] 劉 瓊,耿福海,陳勇航,等.上海不同強(qiáng)度干霾期間氣溶膠垂直分布特征[J]. 中國(guó)環(huán)境科學(xué), 2012,32(2):207-213.Liu Q, Geng F H, Chen Y H, et al. Vertical distribution of aerosols during different intense dry haze periods around Shanghai [J]. China Environmental Science, 2012,32(2):207-213.

    [6] 王 苑,耿福海,陳勇航,等.基于微脈沖激光雷達(dá)的上海浦東地區(qū)不同強(qiáng)度霾研究[J]. 中國(guó)環(huán)境科學(xué), 2013,33(1):21-29.Wang Y, Geng F H, Chen Y H, et al. Study of haze with different intensity categories at Pudong District in Shanghai based on micro- pulse lidar [J]. China Environmental Science, 2013,33(1):21-29.

    [7] 王 靜,牛生杰,許 丹,等.南京一次典型霧霾天氣氣溶膠光學(xué)特性[J]. 中國(guó)環(huán)境科學(xué), 2013,33(2):201-208.Wang J, Niu S J, Xu D, et al. Aerosol optical properties during the typical haze/fog event in Nanjing [J]. China Environmental Science. 2013,33(2):201-208.

    [8] 朱佳雷,王體健,邢 莉,等.江蘇省一次重霾污染天氣的特征和機(jī)理分析[J]. 中國(guó)環(huán)境科學(xué), 2011,31(12):1943-1950.Zhu J L, Wang T J, Xing L, et al. Analysis on the characteristics and mechanism of a heavy haze episode in Jiangsu Province [J]. China Environmental Science. 2011,31(12):1943-1950.

    [9] 吳 兌,廖碧婷,吳 蒙,等.環(huán)首都圈霾和霧的長(zhǎng)期變化特征與典型個(gè)例的近地層輸送條件[J]. 環(huán)境科學(xué)學(xué)報(bào), 2014,34(1):1-11.Wu D, Liao B, Wu M, et al. The long-term trend of haze and fog days and the surface layer transport conditions under haze weather in North China [J]. Acta Scientiae Circumstantiae, 2014,34(1):1-11.

    [10] 薛亦峰,周 震,聶 滕,等.2015年12月北京市空氣重污染過(guò)程分析及污染源排放變化[J]. 環(huán)境科學(xué), 2016,37(5):1593-1601. Xue Y F, Zhou Z, Nie T, et al. Exploring the Severe Haze in Beijing During December, 2015: Pollution Process and Emissions Variation [J]. Environmental Science. 2016,37(5):1593-1601.

    [11] 穆 穆,張人禾.應(yīng)對(duì)灰霾天氣:氣象科學(xué)與技術(shù)大有可為[J]. 中國(guó)科學(xué):地球科學(xué), 2014,44(1):1-2. Mu M, Zhang R H.Dealing with Haze Weather: Great Prospects for Meteorological Science and Technology [J].Chinese Science: Geosciences, 2014,44(1):1-2.

    [12] 董春卿,鄭有飛,武永利,等.邊界層方案對(duì)山西冬季一次靜穩(wěn)天氣PM2.5濃度模擬的影響[J]. 中國(guó)環(huán)境科學(xué), 2016,36(6):1669-1680.Dong C Q, Zheng Y F, Wu Y L, et al. The effects of different planetary boundary layer schemes on PM2.5concentration simulations in winter stable weather of Shanxi [J]. China Environmental Science. 2016, 36(6):1669-1680.

    [13] 李令軍,王占山,張大偉,等.2013~2014年北京大氣重污染特征研究[J]. 中國(guó)環(huán)境科學(xué), 2016,36(1):27-35.Ling J, Wang Z S, Zhang D W, et al. Analysis of heavy air pollution episodes in Beijing during 2013~2014 [J]. China Environmental Science, 2016,36(1):27-35.

    [14] 劉興瑞,馬 嫣,崔芬萍,等.南京北郊一次重污染事件期間PM2.5理化特性及其對(duì)大氣消光的影響[J]. 環(huán)境化學(xué), 2016,35(6):1164- 1171. LiuX R, Ma Y, Cui F P, et al. Physicochemical characteristics of PM2.5and impacts on light extinction during the heavy pollution period at North Suburban Nanjing [J]. Environmental Chemistry, 2016,35(6): 1164-1171.

    [15] 蔡子穎,韓素芹,張 敏,等.天津地區(qū)污染天氣分析中垂直擴(kuò)散指標(biāo)構(gòu)建及運(yùn)用[J]. 環(huán)境科學(xué), 2018,39(6):2548-2556.Cai Z Y, Han S Q, Zhang M, et al. Construction and Application of Vertical Diffusion Index for Analyzing Weather During Pollution Events in Tianjin [J]. Environmental Science, 2018,39(6):2548-2556.

    [16] Tang G, Zhu X, Hu B, et al. Impact of emission controls on air quality in Beijing during APEC 2014: lidar ceilometer observations [J]. Atmospheric Chemistry & Physics, 2015,15(21):12667-12680.

    [17] Han S, Hao T, Zhang Y, et al. Vertical observation and analysis on rapid formation and evolutionary mechanisms of a prolonged haze episode over central-eastern China. Science of total environment [J]. Science of the total environment, 2018,2018(616):135-146.

    [18] 李 夢(mèng),唐貴謙,黃 俊,等.京津冀冬季大氣混合層高度與大氣污染的關(guān)系[J]. 環(huán)境科學(xué), 2015,36(6):1935-1943.Li M, Tang G Q, Huang J, et al. Characteristics of Winter Atmospheric Mixing Layer Height in Beijing-Tianjin-Hebei Region and Their Relationship with the Atmospheric Pollution [J]. Environmental Science, 2015,36(6):1935-1943.

    [19] 蔡子穎,韓素芹,汪 靖,等.基于天氣背景天津地區(qū)重污染天氣特征分析[J]. 環(huán)境科學(xué)學(xué)報(bào), 2017,37(10):3906-3917.Cai Z Y, Han S Q, Wang J, et al. Analysis of synoptic characteristics of heavy pollution in Tianjin based on weather background [J]. Acta Scientiae Circumstantiae, 2017,37(10):3906-3917.

    [20] 周 鑫,丁菊麗,田 偉,等.能見度參數(shù)化方案概述[J]. 氣象研究與應(yīng)用, 2009,30(3):22-25.Zhou X, Ding J, Tian W, et al. Introduction of Visibility Parameterization Scheme [J]. Journal of meteorological research and application. 2009,30(3):22-25.

    [21] 邊 海,張?jiān)7?韓素芹,等.天津市大氣能見度與顆粒物污染的關(guān)系[J]. 中國(guó)環(huán)境科學(xué), 2012,32(3):406-410. Bian H, Han S Q, Zhang Y F, et al. Relationship between atmospheric visibility and particulate matter pollution in Tianjin [J]. China Environmental Science. 2012,32(3):406-410.

    [22] 陳義珍,趙 丹,柴發(fā)合,等.廣州市與北京市大氣能見度與顆粒物質(zhì)量濃度的關(guān)系[J]. 中國(guó)環(huán)境科學(xué), 2010,30(7):967-971.Chen Y Z, Zhao D, Chai F H, et al. Correlation between the atmospheric visibility and aerosol fine particle concentrations in Guangzhou and Beijing [J]. China Environmental Science. 2010,30(7): 967-971.

    [23] 姚 青,蔡子穎,韓素芹,等.2009年秋冬季天津低能見度天氣下氣溶膠污染特征[J]. 氣象, 2012,38(9):1096-1102. Yao Q, Cai Z Y, Han S Q, et al. Study on Pollution Characteristics of aerosols during low visibility weather at Tianjin City in Autumn and Winter 2009 [J]. Meteorological monthly. 2012,38(9):1096-1102.

    [24] 姚 青,韓素芹,蔡子穎,等.天津城區(qū)春季大氣氣溶膠消光特性研究[J]. 中國(guó)環(huán)境科學(xué), 2012,32(5):795-802.Yao Q, Han S Q, Cai Z Y, et al. Study on characteristic of aerosol extinction at Tianjin City in the spring [J]. China Environmental Science, 2012,32(5):795-802.

    [25] 劉西川,高太長(zhǎng),劉 磊,等.降水現(xiàn)象對(duì)大氣消光系數(shù)和能見度的影響[J]. 應(yīng)用氣象學(xué)報(bào), 2010,21(4):433-441. Liu X C, Gao T C, Liu L, et al. Influences of Precipitation on atmospheric Extinction coefficient and Visibility [J]. Journal of applied metrological science, 2010,21(4):433-441.

    [26] 劉寧微,馬雁軍,王揚(yáng)鋒.遼寧中部城市群夏季大氣能見度的觀測(cè)研究[J]. 氣象學(xué)報(bào), 2012,70(4):814-820.Liu N W, Ma Y J, Wang Y F, et al. Observational study of atmospheric visibility in summer in the area of multicities in central Liaoning Province [J]. Acta Meteorological Sinica, 2012,70(4):814-820.

    [27] 饒曉琴,馬學(xué)款,黃 威.2014年12月大氣環(huán)流和天氣分析[J]. 氣象, 2015,41(3):380-387.Wang J K, Gui H L, Ma X K, et al. Analysis of the December 2017 Atmospheric circulation and weather [J]. Meteorological Monthly. 2015,41(3):380-387.

    [28] 熊亞軍,廖曉農(nóng),李梓銘,等.KNN數(shù)據(jù)挖掘算法在北京地區(qū)霾等級(jí)預(yù)報(bào)中的應(yīng)用[J]. 氣象, 2015,41(1):98-104.Xiong Y J, Liao X N, Li Z M, et al. Application of KNN data mining algorithm to haze grade forecasting in Beijing [J]. Meteorological Monthly. 2015,41(1):98-104.

    [29] 趙秀娟,徐 敬,張自銀,等.北京區(qū)域環(huán)境氣象數(shù)值預(yù)報(bào)系統(tǒng)及PM2.5預(yù)報(bào)檢驗(yàn)[J]. 應(yīng)用氣象學(xué)報(bào), 2016,27(2):160-172.Zhao X J, Xu J, Zhang Z Y, et al. Beijing Regional Environmental Metrology prediction system and its performance test of PM2.5concentration [J]. Journal of applied metrological science. 2016,27(2): 160-172.

    [30] 鄧 濤,鄧雪嬌,吳 兌,等.珠三角灰霾數(shù)值預(yù)報(bào)模式與業(yè)務(wù)運(yùn)行評(píng)估[J]. 氣象科技進(jìn)展, 2012,2(6):38-44.Deng T, Deng X J, Wu D, et al. Study on Numerical Forecast Model of Haze over Pearl River Delta Region and Routine Business Assessment [J]. Advances in Meteorological Science and Technology. 2012,2(6): 38-44.

    [31] 李 鋒,朱 彬,安俊嶺,等.2013年12月初長(zhǎng)江三角洲及周邊地區(qū)重霾污染的數(shù)值模擬[J]. 中國(guó)環(huán)境科學(xué), 2015,35(7):1965-1974.Li F, Zhu B, An J L, et al. Modeling study of a severe haze episode occurred over the Yangtze River Delta and its surrounding regions during early December, 2013 [J]. China Environmental Science, 2015, 35(7):1965-1974.

    [32] 蔡子穎,劉愛霞,韓素芹,等.天津低能見度特征初探[J]. 氣象, 2014,40(1):114-118.Cai Z Y, Liu A X, Han S Q, et al. Research of Characteristics on low visibility in Tianjin[J].Meteorological monthly. 2014,40(1):114-118.

    [33] 馬 楠,趙春生,陳 靜,等.基于實(shí)測(cè)PM2.5、能見度和相對(duì)濕度分辨霧霾的新方法[J]. 中國(guó)科學(xué):地球科學(xué), 2015,45(2):227-235. Ma N, Zhao C S, Chen J, et al. A novel method for distinguishing fog and haze based on PM2.5, visibility, and relative humidity [J]. Science China: Earth Sciences, 2014(57):2156–2164.

    [34] 白永清,祁海霞,劉 琳,等.武漢大氣能見度與PM2.5濃度及相對(duì)濕度關(guān)系的非線性分析及能見度預(yù)報(bào)[J]. 氣象學(xué)報(bào), 2016,74(2): 189-199.Bai Y Q, Qi H X, Liu L, et al. Study on the nonlinear relationship among the visibility, PM2.5concentration and relative humidity in Wuhan and visibility prediction [J]. Acta Meteorological sinica. 2016, 74(2):189-199.

    [35] 胡 俊,趙天良,張澤鋒,等.霾污染環(huán)境大氣能見度參數(shù)化方案的改進(jìn)[J]. 環(huán)境科學(xué)研究, 2017,30(11):1680-1688.Hu J, Zhao T L, Zhang Z F, et al. Upgradeding atmospheric visibility parameterization scheme for haze pollution environment [J]. Research of Environmental Sciences,2017,30(11):1680-1688.

    [36] 沈鐵迪,王體健,陳璞瓏,等.南京城區(qū)夏秋季能見度與PM2.5化學(xué)成分的關(guān)系[J]. 中國(guó)環(huán)境科學(xué), 2015,35(3):652-658.Shen T D, Wang T J, Chen P L,et al. Relationship between atmospheric visibility and chemical composition of PM2.5in the summer and autumn of Nanjing [J]. China Environmental Science. 2015,35(3):652-658.

    [37] 古金霞,董海燕,白志鵬,等.天津市顆粒物散射消光特征及化學(xué)組分貢獻(xiàn)研究[J]. 中國(guó)環(huán)境科學(xué), 2012,32(1):17-22.Gu J X, Dong H Y, Bai Z P,et al. Light scattering extinction characterization of particulates and chemical composition contribution in Tianjin [J]. China Environmental Science. 2012,32(1):17-22.

    [38] 尚 倩,李子華,楊 軍,等.南京冬季大氣氣溶膠粒子譜分布及其對(duì)能見度的影響[J]. 環(huán)境科學(xué), 2011,32(9):2750-2760.Shan Q, Li Z H, Yang J, et al. Size Distributions of Aerosol Particles and the Impact on Visibility in Winter of Nanjing [J]. Environmental Science. 2011,32(9):2750-2760.

    [39] 韓素芹,張?jiān)7?李英華,等.天津市春季氣溶膠消光特征和輻射效應(yīng)的數(shù)值模擬[J]. 中國(guó)環(huán)境科學(xué), 2011,31(1):8-12.Han S Q, Zhang Y F, Li Y H, et al. Simulation of extinction and radiant effect of aerosol in spring of Tianjin City [J]. China Environmental Science, 2011,31(1):8-12.

    [40] Mu Q, Lammel G, Gencarelli C N, et al. Regional modelling of polycyclic aromatic hydrocarbons: WRF/Chem-PAH model development and East Asia case studies [J]. Atmospheric Chemistry & Physics, 2017,2017(17):12253-12267.

    [41] Ryu Y H, Hodzic A, Barre J, et al. Quantifying errors in surface ozone predictions associated with clouds over CONUS: A WRF-Chem modeling study using satellite cloud retrievals [J]. Atmospheric Chemistry & Physics, 2018,(18):7509-7525.

    [42] Zhang L, Li Q, Wang T, et al. Combined impacts of nitrous acid and nitryl chloride on lower-tropospheric ozone: new module development in WRF-Chem and application to China [J]. Atmospheric Chemistry & Physics, 2017,(17):9733-9750.

    [43] 中國(guó)氣象局.QX/T 113中華人民共和國(guó)氣象行業(yè)標(biāo)準(zhǔn):霾的觀測(cè)和預(yù)報(bào)等級(jí) [S]. 2010.China Meteorological Administration. QX/T 113. Meteorological industry standard of the People's Republic of China: Haze observation and forecast level

    Research on characteristics of haze weather in Tianjin.

    CAI Zi-ying1,2,3*, YAO Qing1,3, HAN Su-qin3, QIU Xiao-bin3, ZHANG-Min1, WU Bin-gui3, WANG Xue-lian3

    (1.Tianjin Environmental Meteorological Center, Tianjin 300074, China;2.Key Laboratory of Atmospheric Chemistry, China Meteorological Administration, Beijing 100081, China;3.Tianjin Institute of Meteorology, Tianjin 300074, China)., 2019,39(1):50~60

    A research was carried out on characteristics of the haze weather in Tianjin based on the monitoring data of PM2.5mass concentration, visibility and relative humidity in the area from 2014 to 2017, The findings showed that there are five types of haze processes above moderate degree in Tianjin: that in the rear of high pressure, that caused by weak high pressure from north, that caused by low pressure trough, that caused by uniform pressure field and that caused by low pressure before front. As per the current standards, moderate haze generally corresponds to severely polluted weather; severe haze corresponds to severely to untra-severely polluted weather. Moderate to severe haze usually occurs in Level V severely polluted weather, and severe haze occurs in Level VI ultra-severely polluted weather. After "Air Pollution Prevention Action Plan" was implemented in 2013, PM2.5mass concentration and haze days were significantly reduced, by 2017, the number of haze days had decreased by 55% compared with 2013, moderate and above haze days decreased by 50% from 41d in 2013 to 20d in 2017. Based on the monitoring data of PM2.5mass concentration, visibility and relative humidity, the calculation equation of regional visibility can be well constructed. Statistics showed that the correlation coefficient between visibility and real value in the estimation is 0.94, with the relative error 18.6%, identification accuracy rate of the non-haze day is 85%, that of the haze days 95.6%, that of the ultra-slight haze days 83%, that of the slight haze days 78%, that of moderate haze days 93% and that of severe and above haze days 94%, which is greatly applicable for judgement of haze level. This equation was combined with the air quality model to predict the haze level. Inspection of the forecasting product showed: the correlation coefficient between the forecast values and real values of visibility is 0.75, with average forecast value 13.9km and the average real value 14.1km, and the relative error rate 29.6%, the FAC2 (the ratio of the forecast value to two times the real value) is 98.1%, the forecast accuracy rate, forecast missing rate and false rate of the haze days are 81.4%, 18.6% and 20.6%, respectively. Let the tolerance be of Level I, the forecast accuracy rates of ultra-slight haze days, the slight haze days and the moderate and above haze days are 96%, 85% and 69%, respectively, which may effectively support the haze level forecast in Tianjin.

    forecast of haze level;severely polluted weather;Tianjin

    X131.1

    A

    1000-6923(2019)01-0050-11

    蔡子穎(1984-),男,江蘇揚(yáng)州人,南京信息工程大學(xué)碩士研究生,高級(jí)工程師,主要從事大氣環(huán)境和氣候效應(yīng)研究.發(fā)表論文14篇.

    2018-05-24

    國(guó)家自然科學(xué)基金資助項(xiàng)目(41771242,41675018);國(guó)家重點(diǎn)研發(fā)計(jì)劃(2016YFC0203302);中國(guó)氣象科學(xué)研究院中國(guó)氣象局大氣化學(xué)重點(diǎn)實(shí)驗(yàn)室開放課題(2018B06);天津市自然科學(xué)基金資助項(xiàng)目(16JCYBJC21500,17JCYBJC23400)

    * 責(zé)任作者, 高級(jí)工程師, 120078030@163.com

    猜你喜歡
    天津地區(qū)能見度中度
    天津地區(qū)口感型蔬菜品種介紹
    天津地區(qū)黏蟲遷飛發(fā)生特征及監(jiān)測(cè)預(yù)報(bào)
    小米手表
    2005—2017年白云機(jī)場(chǎng)能見度變化特征及其與影響因子關(guān)系研究
    LEEP治療中度宮頸上皮內(nèi)瘤變對(duì)HPV感染的影響
    低能見度下高速公路主動(dòng)誘導(dǎo)技術(shù)的應(yīng)用
    前向散射能見度儀的常見異?,F(xiàn)象處理及日常維護(hù)
    前向散射能見度儀故障實(shí)例分析
    天津地區(qū)地?zé)崃黧w化學(xué)成分垂向演化特征研究
    輿論引導(dǎo)中度的把握
    聲屏世界(2014年6期)2014-02-28 15:18:07
    高清午夜精品一区二区三区| 精品午夜福利在线看| 亚洲国产精品一区三区| 午夜激情福利司机影院| 久久久久久久亚洲中文字幕| 在线看a的网站| 美女主播在线视频| 精品一区在线观看国产| 国产黄色视频一区二区在线观看| 中文字幕人妻熟人妻熟丝袜美| 看免费成人av毛片| 欧美日韩视频精品一区| 男女国产视频网站| 男人爽女人下面视频在线观看| 久久久久久久国产电影| h日本视频在线播放| 五月伊人婷婷丁香| 91精品伊人久久大香线蕉| 极品教师在线视频| 色综合色国产| 亚洲国产精品专区欧美| 亚洲欧美清纯卡通| 欧美一级a爱片免费观看看| 99视频精品全部免费 在线| 日韩精品有码人妻一区| 一区二区三区免费毛片| 丝瓜视频免费看黄片| 岛国毛片在线播放| 26uuu在线亚洲综合色| 综合色丁香网| 午夜免费观看性视频| 天堂中文最新版在线下载| 最黄视频免费看| 人人妻人人添人人爽欧美一区卜 | 久久久a久久爽久久v久久| 网址你懂的国产日韩在线| 久久国产乱子免费精品| 日本av手机在线免费观看| 久久精品国产亚洲网站| 欧美日本视频| 日本与韩国留学比较| 九色成人免费人妻av| 一区二区三区乱码不卡18| 亚洲精品一区蜜桃| 亚洲欧美精品自产自拍| 国产精品欧美亚洲77777| 国产人妻一区二区三区在| 精品一区二区三卡| 亚洲一区二区三区欧美精品| 天美传媒精品一区二区| 99久久人妻综合| 三级经典国产精品| 黑人高潮一二区| 一本久久精品| 日本欧美视频一区| 制服丝袜香蕉在线| 成人影院久久| 美女视频免费永久观看网站| 欧美日本视频| 亚洲欧美成人综合另类久久久| 久久6这里有精品| 欧美成人精品欧美一级黄| 欧美人与善性xxx| 久久99精品国语久久久| 女人十人毛片免费观看3o分钟| 亚洲精品456在线播放app| 内射极品少妇av片p| 亚洲美女黄色视频免费看| 青青草视频在线视频观看| 妹子高潮喷水视频| 三级经典国产精品| 大片电影免费在线观看免费| 简卡轻食公司| 久久av网站| 只有这里有精品99| 亚洲高清免费不卡视频| 涩涩av久久男人的天堂| 我的老师免费观看完整版| 欧美精品国产亚洲| 三级国产精品片| 国产毛片在线视频| 日本与韩国留学比较| 黄色怎么调成土黄色| 91精品国产九色| 成人高潮视频无遮挡免费网站| 婷婷色综合大香蕉| 高清欧美精品videossex| 国产亚洲av片在线观看秒播厂| 九九在线视频观看精品| 五月玫瑰六月丁香| 国产乱来视频区| 亚洲国产精品999| 少妇人妻一区二区三区视频| av在线播放精品| 欧美极品一区二区三区四区| 精品熟女少妇av免费看| 在线观看国产h片| 成年免费大片在线观看| 伊人久久精品亚洲午夜| 99热网站在线观看| 亚洲国产色片| 不卡视频在线观看欧美| 亚洲美女视频黄频| 国产淫片久久久久久久久| 国产伦精品一区二区三区四那| 午夜福利在线在线| 王馨瑶露胸无遮挡在线观看| 国产色婷婷99| 亚洲精品乱码久久久久久按摩| 国产精品女同一区二区软件| 久久精品国产鲁丝片午夜精品| 熟女电影av网| 欧美性感艳星| 中文在线观看免费www的网站| 1000部很黄的大片| 国产精品精品国产色婷婷| 嘟嘟电影网在线观看| 国产一区二区三区综合在线观看 | 老师上课跳d突然被开到最大视频| 制服丝袜香蕉在线| 观看美女的网站| 高清毛片免费看| 人体艺术视频欧美日本| 久久人人爽人人爽人人片va| 男人狂女人下面高潮的视频| 日本黄色日本黄色录像| 国产欧美另类精品又又久久亚洲欧美| 久久久久性生活片| 亚洲av二区三区四区| 成人国产av品久久久| 国产欧美另类精品又又久久亚洲欧美| 国产美女午夜福利| 国产 一区 欧美 日韩| 九九在线视频观看精品| 美女xxoo啪啪120秒动态图| 久久久久性生活片| 久久久久久久久久人人人人人人| 亚洲国产欧美人成| 嘟嘟电影网在线观看| 欧美日韩综合久久久久久| 亚洲一级一片aⅴ在线观看| 2021少妇久久久久久久久久久| 99久久人妻综合| 国产免费一级a男人的天堂| 搡老乐熟女国产| 亚洲精品中文字幕在线视频 | 欧美三级亚洲精品| 噜噜噜噜噜久久久久久91| 男人狂女人下面高潮的视频| 国产精品一区二区在线观看99| 一本—道久久a久久精品蜜桃钙片| 简卡轻食公司| 91狼人影院| 欧美亚洲 丝袜 人妻 在线| 天美传媒精品一区二区| 又黄又爽又刺激的免费视频.| 亚洲第一区二区三区不卡| 日本猛色少妇xxxxx猛交久久| 免费在线观看成人毛片| 你懂的网址亚洲精品在线观看| 精品久久国产蜜桃| 欧美日韩亚洲高清精品| 国产精品一区二区在线观看99| 蜜桃在线观看..| 最黄视频免费看| 大香蕉97超碰在线| 久久精品国产亚洲av涩爱| 日本wwww免费看| av国产精品久久久久影院| 老熟女久久久| 国产有黄有色有爽视频| 五月玫瑰六月丁香| 国产老妇伦熟女老妇高清| 亚洲人成网站在线播| 美女内射精品一级片tv| 91久久精品国产一区二区三区| 亚洲欧美日韩卡通动漫| 美女国产视频在线观看| 女人十人毛片免费观看3o分钟| 欧美丝袜亚洲另类| 黄色欧美视频在线观看| 午夜福利在线观看免费完整高清在| 国产有黄有色有爽视频| 最新中文字幕久久久久| 狂野欧美激情性xxxx在线观看| 国产精品久久久久成人av| 精品人妻熟女av久视频| 黄色视频在线播放观看不卡| 欧美成人精品欧美一级黄| 亚洲精品乱码久久久v下载方式| 欧美激情极品国产一区二区三区 | 深夜a级毛片| 一级毛片电影观看| 久久精品夜色国产| 日日摸夜夜添夜夜添av毛片| 国产成人精品久久久久久| 精品久久国产蜜桃| 在线免费十八禁| 久久精品国产亚洲av涩爱| 一级毛片 在线播放| 久久人人爽人人片av| 国产黄片美女视频| 综合色丁香网| 免费人成在线观看视频色| 亚洲欧美日韩无卡精品| a 毛片基地| 人妻制服诱惑在线中文字幕| 伦理电影免费视频| 国产 一区 欧美 日韩| 日日摸夜夜添夜夜添av毛片| 一级片'在线观看视频| 亚洲精品456在线播放app| 制服丝袜香蕉在线| 国产免费一区二区三区四区乱码| 免费看光身美女| av又黄又爽大尺度在线免费看| 成年免费大片在线观看| av网站免费在线观看视频| 啦啦啦视频在线资源免费观看| 国产免费又黄又爽又色| 国产伦精品一区二区三区视频9| 简卡轻食公司| 亚洲第一区二区三区不卡| 亚洲精品久久午夜乱码| 午夜免费观看性视频| 国产精品av视频在线免费观看| 久久久久精品性色| 全区人妻精品视频| 91久久精品国产一区二区成人| 欧美97在线视频| 在线观看一区二区三区| 久久99热6这里只有精品| 国产av精品麻豆| 成人漫画全彩无遮挡| 一级爰片在线观看| 菩萨蛮人人尽说江南好唐韦庄| 国产精品国产三级国产专区5o| 全区人妻精品视频| 好男人视频免费观看在线| 熟女av电影| 久久久久久伊人网av| 欧美丝袜亚洲另类| 亚洲精品亚洲一区二区| 男人爽女人下面视频在线观看| 国产精品一区二区三区四区免费观看| 亚洲国产精品一区三区| 精品一区二区三卡| 日韩三级伦理在线观看| 99久久人妻综合| 在线观看免费日韩欧美大片 | 成人影院久久| 观看免费一级毛片| 精品久久久久久电影网| 国产中年淑女户外野战色| 国产在线免费精品| 一级二级三级毛片免费看| 黄色欧美视频在线观看| 一本—道久久a久久精品蜜桃钙片| 国产精品一二三区在线看| 欧美日韩一区二区视频在线观看视频在线| 啦啦啦视频在线资源免费观看| 免费看不卡的av| 爱豆传媒免费全集在线观看| 97在线人人人人妻| 99热国产这里只有精品6| 最近2019中文字幕mv第一页| 亚洲经典国产精华液单| 国产黄片视频在线免费观看| 熟女av电影| 2018国产大陆天天弄谢| 中文字幕人妻熟人妻熟丝袜美| 国产 一区精品| 久久久久精品性色| 1000部很黄的大片| 国产精品嫩草影院av在线观看| 国产成人精品福利久久| 偷拍熟女少妇极品色| 一级毛片aaaaaa免费看小| 国产精品av视频在线免费观看| 日本黄色片子视频| 国内少妇人妻偷人精品xxx网站| 国产免费一区二区三区四区乱码| 女人久久www免费人成看片| 男女国产视频网站| 日韩视频在线欧美| 中文精品一卡2卡3卡4更新| 亚洲av在线观看美女高潮| 欧美精品一区二区免费开放| 日韩 亚洲 欧美在线| 夫妻性生交免费视频一级片| 啦啦啦中文免费视频观看日本| 免费看不卡的av| 国产亚洲欧美精品永久| 夜夜爽夜夜爽视频| 日韩中文字幕视频在线看片 | 亚洲国产欧美人成| 久久这里有精品视频免费| 国产成人免费观看mmmm| 老熟女久久久| 亚洲国产av新网站| 国产在线一区二区三区精| 一区在线观看完整版| 国产黄色视频一区二区在线观看| 午夜福利高清视频| 亚洲精品日本国产第一区| 久久久久久久久久人人人人人人| 好男人视频免费观看在线| 久久鲁丝午夜福利片| 男女边吃奶边做爰视频| 一级毛片aaaaaa免费看小| 97超视频在线观看视频| 观看免费一级毛片| 爱豆传媒免费全集在线观看| 一级爰片在线观看| 国产欧美日韩精品一区二区| 五月玫瑰六月丁香| 精品亚洲成国产av| 爱豆传媒免费全集在线观看| 一级爰片在线观看| 中文精品一卡2卡3卡4更新| 精品亚洲成a人片在线观看 | 精品亚洲成国产av| 18禁裸乳无遮挡动漫免费视频| 精品99又大又爽又粗少妇毛片| 亚洲欧美精品专区久久| 男人狂女人下面高潮的视频| 免费观看a级毛片全部| 免费在线观看成人毛片| 一区二区三区精品91| 国产免费视频播放在线视频| 欧美激情极品国产一区二区三区 | 肉色欧美久久久久久久蜜桃| 最近中文字幕高清免费大全6| 国产精品不卡视频一区二区| 国产精品精品国产色婷婷| 18+在线观看网站| 午夜激情久久久久久久| 国产色婷婷99| 一级a做视频免费观看| 最近最新中文字幕大全电影3| 亚洲欧美精品专区久久| 国产精品久久久久久精品古装| 日本与韩国留学比较| 99热全是精品| 亚洲欧美精品专区久久| 嘟嘟电影网在线观看| 街头女战士在线观看网站| 我要看黄色一级片免费的| 在线天堂最新版资源| 午夜老司机福利剧场| 精品少妇黑人巨大在线播放| 久久97久久精品| 日本黄色日本黄色录像| 国产成人免费无遮挡视频| 精品少妇黑人巨大在线播放| 女人十人毛片免费观看3o分钟| 亚洲最大成人中文| 久久久久人妻精品一区果冻| 少妇的逼水好多| 精品久久国产蜜桃| 国产精品人妻久久久久久| 国产精品三级大全| 97在线人人人人妻| 亚洲精品色激情综合| 国产日韩欧美亚洲二区| 日日撸夜夜添| 夜夜爽夜夜爽视频| 一级毛片电影观看| 人妻一区二区av| 噜噜噜噜噜久久久久久91| 97超碰精品成人国产| 中国国产av一级| 中国美白少妇内射xxxbb| 18禁在线播放成人免费| 国产欧美亚洲国产| 嫩草影院入口| 精品国产露脸久久av麻豆| 人人妻人人添人人爽欧美一区卜 | 成人影院久久| 国产探花极品一区二区| av在线播放精品| 国产成人精品福利久久| 97热精品久久久久久| 精品国产一区二区三区久久久樱花 | 欧美性感艳星| 国产欧美亚洲国产| 久久午夜福利片| 欧美激情极品国产一区二区三区 | 午夜福利高清视频| 人妻系列 视频| 中文欧美无线码| 嘟嘟电影网在线观看| 国产一区二区在线观看日韩| 日日啪夜夜爽| 草草在线视频免费看| 极品少妇高潮喷水抽搐| 99久久中文字幕三级久久日本| videos熟女内射| 美女国产视频在线观看| 欧美97在线视频| 亚洲国产欧美人成| 亚洲,欧美,日韩| 亚洲欧美日韩卡通动漫| 人体艺术视频欧美日本| 亚洲国产色片| 亚洲精品,欧美精品| 一级av片app| 亚洲精品日韩av片在线观看| 欧美高清性xxxxhd video| 少妇丰满av| 久久99热6这里只有精品| 精品国产一区二区三区久久久樱花 | 黄色一级大片看看| 国产精品女同一区二区软件| 另类亚洲欧美激情| 久久久国产一区二区| tube8黄色片| 97超视频在线观看视频| 少妇精品久久久久久久| 成人国产麻豆网| 亚洲精品成人av观看孕妇| 亚洲精品视频女| 国产黄色视频一区二区在线观看| 国产女主播在线喷水免费视频网站| 国内精品宾馆在线| 黑丝袜美女国产一区| 欧美成人一区二区免费高清观看| 免费av不卡在线播放| av国产精品久久久久影院| 日韩,欧美,国产一区二区三区| 久久人人爽av亚洲精品天堂 | 青春草亚洲视频在线观看| 久久国产精品男人的天堂亚洲 | 国产成人免费无遮挡视频| 欧美激情国产日韩精品一区| 欧美性感艳星| 久久婷婷青草| 免费久久久久久久精品成人欧美视频 | 99久国产av精品国产电影| 夫妻性生交免费视频一级片| 欧美zozozo另类| 国产成人午夜福利电影在线观看| a级毛片免费高清观看在线播放| 少妇的逼水好多| 黑人猛操日本美女一级片| 青春草视频在线免费观看| 2018国产大陆天天弄谢| av专区在线播放| 国产精品三级大全| 久热这里只有精品99| 激情五月婷婷亚洲| 中文资源天堂在线| 久久ye,这里只有精品| 黄色欧美视频在线观看| 777米奇影视久久| 晚上一个人看的免费电影| 国产永久视频网站| 日本欧美国产在线视频| 免费观看在线日韩| 精品久久久精品久久久| 亚洲国产高清在线一区二区三| 丰满少妇做爰视频| 成人免费观看视频高清| 免费不卡的大黄色大毛片视频在线观看| 乱码一卡2卡4卡精品| 国产免费福利视频在线观看| 免费观看的影片在线观看| 日韩av在线免费看完整版不卡| 熟女电影av网| 18禁在线无遮挡免费观看视频| 久久亚洲国产成人精品v| 亚洲内射少妇av| 一个人免费看片子| 最新中文字幕久久久久| 晚上一个人看的免费电影| 欧美精品一区二区免费开放| 国产精品蜜桃在线观看| 欧美日韩视频精品一区| 嘟嘟电影网在线观看| 欧美精品一区二区免费开放| 国产精品伦人一区二区| 大话2 男鬼变身卡| 国产淫片久久久久久久久| 日本黄大片高清| 国内精品宾馆在线| 永久免费av网站大全| 成年人午夜在线观看视频| 精品国产一区二区三区久久久樱花 | 精品视频人人做人人爽| 99热国产这里只有精品6| 视频中文字幕在线观看| 国产免费又黄又爽又色| 菩萨蛮人人尽说江南好唐韦庄| 香蕉精品网在线| 久久99蜜桃精品久久| 少妇熟女欧美另类| 男女啪啪激烈高潮av片| 国产亚洲精品久久久com| 欧美三级亚洲精品| 日本欧美国产在线视频| 亚洲欧美一区二区三区黑人 | 国产视频首页在线观看| 久久久精品免费免费高清| 大片电影免费在线观看免费| 新久久久久国产一级毛片| 久久久久视频综合| 男人添女人高潮全过程视频| 亚洲av福利一区| 亚洲美女搞黄在线观看| 大话2 男鬼变身卡| 26uuu在线亚洲综合色| 亚洲成人一二三区av| 亚洲真实伦在线观看| 亚洲欧美中文字幕日韩二区| 我的老师免费观看完整版| 亚洲va在线va天堂va国产| 一级毛片aaaaaa免费看小| 各种免费的搞黄视频| 直男gayav资源| 精品久久久久久久久亚洲| 成人二区视频| 亚洲欧美日韩无卡精品| 美女视频免费永久观看网站| av国产免费在线观看| 精品99又大又爽又粗少妇毛片| 最新中文字幕久久久久| 久久精品国产鲁丝片午夜精品| 中文字幕av成人在线电影| 亚洲精品第二区| 欧美精品人与动牲交sv欧美| 久久久久网色| 亚洲激情五月婷婷啪啪| 免费在线观看成人毛片| 亚洲精品久久久久久婷婷小说| 国产精品国产三级国产专区5o| 国产精品成人在线| 国产精品久久久久久精品古装| 观看av在线不卡| 丝瓜视频免费看黄片| av免费在线看不卡| 成人美女网站在线观看视频| 一边亲一边摸免费视频| 免费看不卡的av| 亚洲国产最新在线播放| 婷婷色麻豆天堂久久| av在线app专区| 美女xxoo啪啪120秒动态图| 亚洲电影在线观看av| 亚洲精品中文字幕在线视频 | 最近最新中文字幕免费大全7| 久久av网站| 一本一本综合久久| 哪个播放器可以免费观看大片| 美女中出高潮动态图| 少妇猛男粗大的猛烈进出视频| 久久精品人妻少妇| 中文字幕免费在线视频6| 亚洲内射少妇av| 国产欧美亚洲国产| 日韩 亚洲 欧美在线| 久久久成人免费电影| 成人综合一区亚洲| 一边亲一边摸免费视频| 中国美白少妇内射xxxbb| 色视频在线一区二区三区| 岛国毛片在线播放| 午夜老司机福利剧场| 高清毛片免费看| 色吧在线观看| 久久婷婷青草| 一级a做视频免费观看| 久久毛片免费看一区二区三区| 少妇裸体淫交视频免费看高清| 亚洲最大成人中文| 国产在视频线精品| 身体一侧抽搐| 51国产日韩欧美| 91aial.com中文字幕在线观看| 高清毛片免费看| 免费不卡的大黄色大毛片视频在线观看| 大又大粗又爽又黄少妇毛片口| 在线观看三级黄色| 日韩欧美一区视频在线观看 | 哪个播放器可以免费观看大片| 纯流量卡能插随身wifi吗| 日本vs欧美在线观看视频 | 中文字幕久久专区| 精品久久国产蜜桃| 亚洲性久久影院| 一二三四中文在线观看免费高清| 十八禁网站网址无遮挡 | av不卡在线播放| 午夜福利高清视频| 91精品国产九色| 18禁裸乳无遮挡动漫免费视频| 久久这里有精品视频免费| 成人漫画全彩无遮挡| 十八禁网站网址无遮挡 | 国产亚洲精品久久久com| 男的添女的下面高潮视频| 性色av一级| freevideosex欧美| 纯流量卡能插随身wifi吗| 亚洲一级一片aⅴ在线观看| 国产av一区二区精品久久 | 亚洲精品国产av蜜桃| 七月丁香在线播放| 国产在线男女| 91在线精品国自产拍蜜月| 国产有黄有色有爽视频| 51国产日韩欧美| 少妇丰满av| 夜夜爽夜夜爽视频| 舔av片在线| 波野结衣二区三区在线|