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

    石家莊地區(qū)能見(jiàn)度變化特征及其與相對(duì)濕度和顆粒物濃度的關(guān)系

    2016-09-22 01:08:10王曉敏韓軍彩鈐偉妙岳艷霞
    干旱氣象 2016年4期
    關(guān)鍵詞:能見(jiàn)度實(shí)測(cè)值石家莊

    王曉敏,韓軍彩,陳 靜,鈐偉妙,岳艷霞

    (河北省石家莊市氣象局,河北 石家莊 050081)

    ?

    石家莊地區(qū)能見(jiàn)度變化特征及其與相對(duì)濕度和顆粒物濃度的關(guān)系

    王曉敏,韓軍彩,陳靜,鈐偉妙,岳艷霞

    (河北省石家莊市氣象局,河北石家莊050081)

    利用1980—2013年石家莊地區(qū)12個(gè)氣象臺(tái)站能見(jiàn)度資料,結(jié)合相對(duì)濕度和PM2.5、PM10濃度數(shù)據(jù),分析了石家莊地區(qū)能見(jiàn)度的時(shí)空分布特征,通過(guò)研究能見(jiàn)度與相對(duì)濕度和PM2.5、PM10濃度的關(guān)系,建立大氣能見(jiàn)度的多元非線性預(yù)報(bào)模型。結(jié)果表明:(1)1980年以來(lái)石家莊地區(qū)年平均能見(jiàn)度以-1.0 km·(10 a)-1的速率呈下降趨勢(shì),夏季下降趨勢(shì)最明顯,春季下降趨勢(shì)最?。?2)1998年前后石家莊地區(qū)能見(jiàn)度變化較大,1999—2013年平均能見(jiàn)度較1980—1998年下降了15.3%,且空間變化也較明顯,1998年之前分別在中北部和中南部存在2個(gè)高值中心,在市區(qū)和趙縣存在2個(gè)低值中心,1998年之后則呈由東向西逐漸遞減的分布形勢(shì);(3)能見(jiàn)度與相對(duì)濕度存在顯著的指數(shù)函數(shù)關(guān)系,而與PM2.5和PM10濃度均呈冪函數(shù)關(guān)系。據(jù)此建立的能見(jiàn)度與相對(duì)濕度和PM2.5、PM10濃度的多元非線性擬合模型能較好地反映能見(jiàn)度的變化規(guī)律,并對(duì)能見(jiàn)度具有一定的預(yù)報(bào)能力。

    能見(jiàn)度;相對(duì)濕度;顆粒物濃度;時(shí)空特征;非線性擬合

    引 言

    能見(jiàn)度是大氣透明度的表征[1],可以反映城市的大氣環(huán)境質(zhì)量[2]。隨著經(jīng)濟(jì)的發(fā)展,城市大氣污染越來(lái)越嚴(yán)重,能見(jiàn)度也隨之發(fā)生顯著變化[3]。影響大氣能見(jiàn)度的因素包括自然因素和人為活動(dòng),自然因素主要指降水、霧、霾、揚(yáng)塵等天氣現(xiàn)象,而人為活動(dòng)包括煤炭、化工等工業(yè)生產(chǎn)的煙氣排放、汽車尾氣、焚燒等[4]。由于能見(jiàn)度觀測(cè)記錄開(kāi)始較早,資料時(shí)間序列較長(zhǎng),可以通過(guò)分析能見(jiàn)度的長(zhǎng)期變化趨勢(shì)來(lái)研究空氣質(zhì)量的變化[5]。早在1960年代,國(guó)外學(xué)者開(kāi)展了能見(jiàn)度變化趨勢(shì)的研究,提出了累積百分率、Ridit分析法等研究方法[6-7]。國(guó)內(nèi)大量學(xué)者也相繼對(duì)能見(jiàn)度及其影響因子展開(kāi)研究[8-20]:付桂琴等[17]分析了河北低能見(jiàn)度的區(qū)域分布和變化趨勢(shì)特征,認(rèn)為相對(duì)濕度與能見(jiàn)度呈顯著負(fù)相關(guān);劉寧微等[18]對(duì)遼寧中部城市群夏季大氣能見(jiàn)度變化特征進(jìn)行分析,并討論了細(xì)顆粒物、水汽和風(fēng)速等要素對(duì)能見(jiàn)度的影響;吳兌等[19]分析指出造成廣州低能見(jiàn)度年代際變化的主要原因是人類活動(dòng)產(chǎn)生的細(xì)粒子;蔡子穎等[20]分析表明,相對(duì)濕度、PM2.5和PM10質(zhì)量濃度與能見(jiàn)度關(guān)系密切,并嘗試?yán)么髿饣瘜W(xué)模式和擬合公式對(duì)能見(jiàn)度進(jìn)行等級(jí)模擬和預(yù)報(bào)。華北地區(qū)是中國(guó)四大灰霾污染嚴(yán)重區(qū)域之一[21],石家莊由于其特殊的地理位置和地形地貌,成為華北地區(qū)大氣污染嚴(yán)重、霧霾等低能見(jiàn)度事件多發(fā)的典型代表性城市[22]。前人研究主要集中在能見(jiàn)度變化特征及與氣象條件之間的關(guān)系,以及基于模式的等級(jí)預(yù)報(bào),由于擬合公式和相對(duì)濕度模擬的誤差,能見(jiàn)度等級(jí)預(yù)報(bào)準(zhǔn)確率仍有待提高。本文利用近30 a石家莊地區(qū)能見(jiàn)度資料,分析能見(jiàn)度的時(shí)空變化特征,通過(guò)探討相對(duì)濕度、PM2.5和PM10濃度與能見(jiàn)度的關(guān)系,利用較復(fù)雜的非線性方法進(jìn)行量化分析和擬合,從而進(jìn)行能見(jiàn)度預(yù)報(bào)預(yù)測(cè),可為大氣能見(jiàn)度的業(yè)務(wù)化定量預(yù)報(bào)提供方法和參考。

    1 資料與方法

    1980年以前能見(jiàn)度觀測(cè)采用0~9的等級(jí)記錄方法,自1980年開(kāi)始按照《地面氣象觀測(cè)規(guī)范(1979年版)》規(guī)定,以km為記錄單位。2014年之后,行唐、高邑、元氏、鹿泉4個(gè)氣象臺(tái)站的能見(jiàn)度改為自動(dòng)在線連續(xù)觀測(cè),觀測(cè)方法與之前不同。為便于分析,資料選取時(shí)段為1980—2013年,區(qū)域?yàn)槭仪f地區(qū),剔除該時(shí)段能見(jiàn)度資料缺測(cè)較多的5個(gè)站點(diǎn),共選取了12個(gè)站點(diǎn)(圖1)。能見(jiàn)度采用人工每天08:00、14:00、20:00(北京時(shí),下同)3個(gè)時(shí)次的觀測(cè)數(shù)據(jù),日能見(jiàn)度為3個(gè)時(shí)次的算術(shù)平均,區(qū)域能見(jiàn)度為12個(gè)站點(diǎn)的算術(shù)平均。另外,還使用了2014—2015年石家莊市氣象站能見(jiàn)度和相對(duì)濕度逐時(shí)數(shù)據(jù)(2014年石家莊市國(guó)家基本氣象站改為自動(dòng)觀測(cè))。

    PM2.5和PM10濃度數(shù)據(jù)來(lái)源于2014年1月—2015年8月石家莊市環(huán)境監(jiān)測(cè)站逐時(shí)監(jiān)測(cè)數(shù)據(jù),為市區(qū)國(guó)控點(diǎn)數(shù)據(jù)經(jīng)過(guò)質(zhì)量控制后的平均值。利用1980—2013年河北省能源消耗總量和煤炭占能源消耗總量百分比的逐年數(shù)據(jù)(由河北省統(tǒng)計(jì)局提供)和2000—2010年每5 a統(tǒng)計(jì)的工業(yè)源廢棄物總排放量和污染物總排放效率數(shù)據(jù)(來(lái)源于河北省統(tǒng)計(jì)年鑒),以及石家莊市統(tǒng)計(jì)局提供的1980—2010年每5 a統(tǒng)計(jì)的3大產(chǎn)業(yè)占比數(shù)據(jù)。冷空氣發(fā)生次數(shù)采用1980—2013年石家莊地區(qū)12個(gè)氣象觀測(cè)站(與能見(jiàn)度站點(diǎn)一致)逐日氣溫資料,區(qū)域冷空氣發(fā)生次數(shù)為12個(gè)站點(diǎn)的算術(shù)平均。

    圖1 石家莊地區(qū)12個(gè)氣象觀測(cè)站分布(陰影為地形高度,單位:km)

    2 能見(jiàn)度時(shí)空變化特征

    2.1時(shí)間分布

    圖2給出1980—2013年石家莊地區(qū)平均能見(jiàn)度的年變化,可以看出,石家莊地區(qū)年均能見(jiàn)度氣候值(1981—2010年)為11.9 km,1980年以來(lái)平均能見(jiàn)度整體呈顯著下降趨勢(shì)(通過(guò)α=0.001的顯著性檢驗(yàn)),變化率為-1.02 km·(10 a)-1,從3 a滑動(dòng)平均來(lái)看,能見(jiàn)度的變化具有一定的階段性,1980—1990年代能見(jiàn)度持續(xù)下降明顯,2000—2007年波動(dòng)變化,無(wú)明顯趨勢(shì),而后呈上升趨勢(shì)。四季平均能見(jiàn)度也均呈下降趨勢(shì)(圖略),其中夏季下降趨勢(shì)最明顯,趨勢(shì)變化率為-1.3 km·(10 a)-1,秋季次之,春季下降趨勢(shì)最小,為-0.5 km·(10 a)-1。

    圖2 1980—2013年石家莊地區(qū)平均能見(jiàn)度的年變化

    1—5月隨著氣溫回暖,大氣湍流活動(dòng)增強(qiáng),石家莊地區(qū)能見(jiàn)度逐漸增大(圖3),5月達(dá)到全年最大值(平均13.6 km);6—8月隨著雨季的來(lái)臨,空氣濕度增大,能見(jiàn)度下降,8月平均值降至11.2 km;

    圖3 1980—2013年石家莊地區(qū)各月平均能見(jiàn)度的年變化(單位:km)

    9月降水明顯減弱,空氣濕度下降,能見(jiàn)度稍有上升(平均12.2 km);10—12月氣溫下降,大氣層結(jié)較穩(wěn)定,能見(jiàn)度下降,12月降至全年最低值(平均10.4 km)。另外,1998年前后能見(jiàn)度變化較大,1998年之前各月平均能見(jiàn)度均≥11 km,1980—1998年平均能見(jiàn)度達(dá)12.7 km;1998年之后能見(jiàn)度明顯下降,1999—2013年平均能見(jiàn)度為10.8 km,較前期下降了15.3%,其中7月—翌年2月能見(jiàn)度<10 km。

    2.2空間分布

    從1980—2013年平均能見(jiàn)度空間分布來(lái)看(圖4a),石家莊地區(qū)平均能見(jiàn)度基本在12 km左右,其中在正定縣和新樂(lè)縣交界、贊皇縣分別存在2個(gè)高值中心,平均能見(jiàn)度分別為14 km和13.3 km;而在鹿泉市和市區(qū)交界、趙縣分別有2個(gè)低值中心,平均能見(jiàn)度僅9.8 km。四季平均能見(jiàn)度的空間分布與年平均能見(jiàn)度基本相似(圖略),均在正定縣和新樂(lè)縣交界處、贊皇縣存在2個(gè)高值中心,而在市區(qū)與鹿泉市交界處、趙縣存在2個(gè)低值中心。不同的是,從春季到冬季平均能見(jiàn)度逐漸降低,高值中心范圍逐漸縮小,冬季贊皇縣的高值中心基本消失;低值中心的范圍逐漸擴(kuò)大,秋、冬季趙縣的低值中心已擴(kuò)大至趙縣、藁城市、無(wú)極縣及以東地區(qū),而冬季鹿泉市與市區(qū)交界的低值中心也擴(kuò)大至平山、靈壽、欒城、井陘等縣(市)區(qū)。

    圖4 1980—2013年石家莊地區(qū)不同時(shí)期平均能見(jiàn)度空間分布(單位:km)

    由前面分析可知,1998年前后石家莊地區(qū)能見(jiàn)度發(fā)生明顯變化,因此以1998年為界分析1980—1998年(圖4b)和1999—2013年(圖4c)平均能見(jiàn)度空間分布。發(fā)現(xiàn),1980—1998年與1980—2013年整個(gè)時(shí)段平均能見(jiàn)度的空間分布較相似,也在正定縣和新樂(lè)縣交界、贊皇縣分別存在2個(gè)高值中心,平均能見(jiàn)度分別達(dá)16.8 km和16.1 km;而在鹿泉市與市區(qū)交界、趙縣分別為2個(gè)低值中心,平均能見(jiàn)度分別為9 km和10.3 km。不同之處是,1980—1998年在井陘縣和平山縣能見(jiàn)度也較高,平均能見(jiàn)度為13 km(圖4b)。1998年之后,除東部地區(qū)無(wú)變化和市區(qū)與鹿泉市交界處略增加外,其余地區(qū)平均能見(jiàn)度均較前期有所減少,且能見(jiàn)度由東向西逐漸遞減,西部靈壽、平山、井陘、贊皇等縣(市)區(qū)平均能見(jiàn)度僅9 km,而東部縣(市)區(qū)平均能見(jiàn)度在12 km左右,在市區(qū)和正定縣交界處有一大值中心,平均能見(jiàn)度為13 km,變化最大的是贊皇縣(圖4c)。

    3 能見(jiàn)度變化成因

    王喜全等[23]認(rèn)為,能見(jiàn)度變化是多方面因素綜合造成的結(jié)果,不僅與經(jīng)濟(jì)發(fā)展相伴隨的污染物排放增加有關(guān),近年來(lái)冷空氣活動(dòng)減弱也是一個(gè)不容忽視的自然因素。研究表明[24-25],石家莊地區(qū)屬煤煙型污染,PM2.5和PM10粒子的主要來(lái)源為燃煤。由河北省能源消耗總量及產(chǎn)業(yè)結(jié)構(gòu)可以看出(圖5a),1980年以來(lái),河北省能源消耗總量呈持續(xù)上升趨勢(shì),尤其在2000年之后增速明顯;煤炭在能源消耗總量的比例在90%附近波動(dòng),2000年以前相對(duì)較低,均在90%以下,且1985—1999年期間較穩(wěn)定,而后波動(dòng)增加,至2009年開(kāi)始持續(xù)減小。石家莊地區(qū)能見(jiàn)度迅速下降與同期煤炭比重增加密切相關(guān),2010年之后隨著能源消耗結(jié)構(gòu)的調(diào)整,煤炭比重明顯下降,這一時(shí)期的能見(jiàn)度也有所回升。

    圖5 1980—2013年河北省能源消耗總量和煤炭占能源消耗總量的比重(a)及石家莊3大產(chǎn)業(yè)所占比重(b)變化

    由圖5b可以看出,1980—1990年代中期石家莊第二產(chǎn)業(yè)的比重迅速增加,而后持續(xù)緩慢增加;第一產(chǎn)業(yè)比重在1990年代以前變化不大,而后持續(xù)減??;第三產(chǎn)業(yè)比重波動(dòng)不大。第二產(chǎn)業(yè)比重的持續(xù)增加和第一產(chǎn)業(yè)比重的持續(xù)減小,使得2000年之后工業(yè)源廢棄物和污染物排放效率明顯增多(表1)[26],對(duì)石家莊地區(qū)空氣質(zhì)量產(chǎn)生較大影響。能見(jiàn)度空間分布形勢(shì)的改變與城市化和工業(yè)化發(fā)展有很大關(guān)系。由于石家莊地區(qū)特殊的地形特征,西北部的山區(qū)地帶,在1990年代之前人煙稀少,工業(yè)化水平較低,而市區(qū)人口密度大,經(jīng)濟(jì)發(fā)展較快[5],因此在1998年之前山區(qū)能見(jiàn)度較高,而市區(qū)能見(jiàn)度偏低;1990年代之后隨著經(jīng)濟(jì)和人口的迅速發(fā)展,西北部縣市礦山開(kāi)采日益增加(圖6),且能源消耗增大,環(huán)境污染加重[27],形成了能見(jiàn)度東高西低的分布形勢(shì)。

    表1 2000年以來(lái)石家莊工業(yè)源廢棄物及污染物排放效率

    圖6 石家莊地區(qū)露天礦山區(qū)域分布

    按照《中華人民共和國(guó)國(guó)家標(biāo)準(zhǔn)(GB/T 20484-2006):冷空氣等級(jí)》,統(tǒng)計(jì)了1980—2013年石家莊地區(qū)發(fā)生中等強(qiáng)度以上冷空氣的次數(shù)。從空間分布來(lái)看(圖7a),石家莊地區(qū)中等強(qiáng)度以上冷空氣發(fā)生次數(shù)的空間分布以市區(qū)—平山一帶為中心分別向東北、西南、東南方向逐漸增加,東部辛集也存在一低值區(qū),總體來(lái)看,石家莊外圍地區(qū)冷空氣活動(dòng)相對(duì)頻繁。其中,新樂(lè)—正定—市區(qū)—贊皇一帶冷空氣發(fā)生次數(shù)的分布形勢(shì)與該地區(qū)能見(jiàn)度空間分布(圖4a)一致,市區(qū)為冷空氣活動(dòng)最差,平均能見(jiàn)度為低值中心,由市區(qū)分別向正定—新樂(lè)和贊皇方向冷空氣活動(dòng)加強(qiáng),平均能見(jiàn)度也呈增加趨勢(shì)。由圖7b可看出,1980年以來(lái)石家莊地區(qū)中等強(qiáng)度以上冷空氣發(fā)生次數(shù)整體呈顯著下降趨勢(shì)(通過(guò)α=0.001的顯著性檢驗(yàn)),與能見(jiàn)度的整體變化特征一致,表明冷空氣活動(dòng)減少不利于污染物的排放和稀釋,致使大氣能見(jiàn)度降低。

    圖7 1980—2013年石家莊地區(qū)中等強(qiáng)度以上冷空氣平均發(fā)生次數(shù)的空間分布(a,單位:次)及年際變化(b)

    綜上所述,石家莊地區(qū)能見(jiàn)度的時(shí)空變化與工業(yè)發(fā)展和污染物排放密切相關(guān),同時(shí)也受冷空氣活動(dòng)的影響,具體哪類因素占主要地位,還需通過(guò)數(shù)值模擬等方法進(jìn)行更加深入的分析。

    4 能見(jiàn)度與相對(duì)濕度和顆粒物濃度的相關(guān)模型

    研究表明,當(dāng)大氣能見(jiàn)度較差時(shí),大氣中的氣溶膠和相對(duì)濕度對(duì)能見(jiàn)度的影響比重較大[13]。相對(duì)濕度較低時(shí),粒徑<2 μm的細(xì)粒子對(duì)大氣能見(jiàn)度作用明顯[28],但隨著相對(duì)濕度的增大,氣溶膠粒子吸濕增長(zhǎng),使得粒徑>2 μm的粗粒子對(duì)能見(jiàn)度影響的貢獻(xiàn)逐漸增大[29]。因此,需同時(shí)考慮PM10、PM2.5和相對(duì)濕度對(duì)大氣能見(jiàn)度的影響,建立能見(jiàn)度預(yù)報(bào)模型。圖8給出2014年石家莊市逐時(shí)能見(jiàn)度與相對(duì)濕度和PM2.5濃度、PM10濃度的散點(diǎn)圖(樣本數(shù)共8 731個(gè)),并進(jìn)行了一元非線性擬合。發(fā)現(xiàn),能見(jiàn)度與相對(duì)濕度呈指數(shù)函數(shù)關(guān)系,而與PM2.5濃度、PM10濃度呈冪函數(shù)關(guān)系,回歸方程分別為:

    (1)

    式中:x1、x2、x3分別為相對(duì)濕度(%),PM2.5濃度(mg·m-3)和PM10濃度(mg·m-3),y1、y2、y3為能見(jiàn)度(km)。其中,能見(jiàn)度與相對(duì)濕度、PM2.5濃度及PM10濃度的相關(guān)系數(shù)R分別為-0.54、-0.71和-0.48,均通過(guò)α=0.001的顯著性檢驗(yàn)。可見(jiàn),能見(jiàn)度與PM2.5濃度的相關(guān)性更高。

    圖8 2014年石家莊市逐時(shí)能見(jiàn)度與相對(duì)濕度、PM2.5濃度、PM10濃度的擬合關(guān)系

    基于以上分析結(jié)果,利用1stOpt軟件對(duì)石家莊大氣能見(jiàn)度進(jìn)行多元非線性擬合,擬合公式如下:

    (2)

    式(2)中,y為大氣能見(jiàn)度(km),x1、x2分別表示PM2.5與PM10質(zhì)量濃度(mg·m-3),x3為相對(duì)濕度(%),a~g為系數(shù)。由于不同相對(duì)濕度條件下,PM2.5和PM10質(zhì)量濃度對(duì)能見(jiàn)度影響的貢獻(xiàn)不同,表2給出不同相對(duì)濕度條件下模型的擬合參數(shù)及均方差(RMSE)、相關(guān)系數(shù)(R),相關(guān)系數(shù)均通過(guò)α=0.001的顯著性檢驗(yàn)。

    表2 不同相對(duì)濕度下大氣能見(jiàn)度與相對(duì)濕度、PM2.5和PM10質(zhì)量濃度的非線性擬合系數(shù)

    為驗(yàn)證公式(2)大氣能見(jiàn)度擬合的合理性,利用2014年1月(代表冬季)和7月(代表夏季)相對(duì)濕度、PM2.5與PM10質(zhì)量濃度逐時(shí)數(shù)據(jù)(樣本數(shù)均為744個(gè)),根據(jù)表2中的系數(shù),對(duì)大氣能見(jiàn)度擬合,并結(jié)合同期逐時(shí)能見(jiàn)度資料,進(jìn)行擬合值與實(shí)測(cè)值的對(duì)比分析(圖9)。可以看出,2014年1月能見(jiàn)度擬合值與實(shí)測(cè)值隨時(shí)間變化吻合較好,擬合值大都比實(shí)測(cè)值偏高,且能見(jiàn)度越低誤差相對(duì)越小(圖9a、圖9b),擬合值與實(shí)測(cè)值相關(guān)系數(shù)達(dá)0.85(通過(guò)α=0.001的顯著性檢驗(yàn)),均方差為3.1 km,均方差與觀測(cè)值比值為0.79。7月能見(jiàn)度擬合值與實(shí)測(cè)值隨時(shí)間變化吻合也較好(圖9c、圖9d),擬合值大部分比實(shí)測(cè)值偏低,但能見(jiàn)度<10 km的擬合值大部分比實(shí)測(cè)值偏高,其相關(guān)系數(shù)為0.91(通過(guò)α=0.001的顯著性檢驗(yàn)),均方差為4.9 km,均方差與觀測(cè)值比值為0.41,誤差較1月略大。總體而言,(2)式對(duì)石家莊大氣能見(jiàn)度有一定的模擬能力,可利用PM10、PM2.5濃度及相對(duì)濕度對(duì)大氣能見(jiàn)度進(jìn)行預(yù)報(bào)預(yù)測(cè)。

    圖9 2014年1月(a,b)、7月(c,d)石家莊能見(jiàn)度的擬合值和實(shí)測(cè)值的時(shí)間演變(a,c)及二者散點(diǎn)圖(b,d)

    為提高能見(jiàn)度的業(yè)務(wù)預(yù)報(bào)水平,將上述建立的石家莊能見(jiàn)度多元非線性擬合公式應(yīng)用到業(yè)務(wù)預(yù)報(bào)工作中。PM2.5、PM10濃度采用本地建立的動(dòng)態(tài)多元回歸模型的預(yù)報(bào)值,該模型采用多元線性逐步回歸方法,基于前期30 d的地面氣象觀測(cè)資料和污染物濃度數(shù)據(jù)建立模型,對(duì)未來(lái)7 d污染物濃度進(jìn)行預(yù)測(cè)。該模型可根據(jù)不同天氣形勢(shì)或者季節(jié)變化不定時(shí)地修訂,使污染物預(yù)報(bào)值與實(shí)測(cè)值更接近。相對(duì)濕度采用T639預(yù)報(bào)的2 m相對(duì)濕度,PM2.5、PM10濃度和相對(duì)濕度每日預(yù)報(bào)時(shí)次為02:00、05:00、08:00、11:00、14:00、17:00、20:00、23:00,間隔為3 h。下面以2015年8月16—18日個(gè)例來(lái)評(píng)估公式(2)的預(yù)報(bào)準(zhǔn)確性。

    2015年8月16—18日石家莊市空氣質(zhì)量逐漸變差,能見(jiàn)度也逐漸隨之下降,其中16日為輕度污染,17—18日為中度污染,18日下午14:00開(kāi)始地面風(fēng)速增大到1.7 m·s-1以上,18:00開(kāi)始出現(xiàn)降雨,空氣質(zhì)量好轉(zhuǎn),能見(jiàn)度也隨之升高。圖10是2015年8月16—18日石家莊能見(jiàn)度的預(yù)報(bào)值與實(shí)測(cè)值對(duì)比結(jié)果,可以看出預(yù)報(bào)值與實(shí)測(cè)值變化趨勢(shì)基本一致,二者相關(guān)系數(shù)達(dá)0.93(通過(guò)α=0.001的顯著性檢驗(yàn)),但大部分預(yù)報(bào)值比實(shí)測(cè)值偏低,尤其是能見(jiàn)度較高時(shí)段二者差異明顯偏大,總體上該模型對(duì)能見(jiàn)度具有一定的預(yù)報(bào)能力。

    圖10 2015年8月16—18日石家莊市逐3 h能見(jiàn)度預(yù)報(bào)值與實(shí)測(cè)值對(duì)比

    5 結(jié) 論

    (1)1980—2013年,石家莊地區(qū)年均能見(jiàn)度為11.9 km,總體上以-1.0 km·(10 a)-1的速率呈下降趨勢(shì)。四季平均能見(jiàn)度也均呈下降趨勢(shì),其中夏季下降趨勢(shì)最明顯,春季下降趨勢(shì)最小。5月平均能見(jiàn)度為全年最大值(13.6 km),12月最小(10.4 km)。1998年前后能見(jiàn)度變化較大,1998年之前各月平均能見(jiàn)度均>11 km,之后明顯下降。

    (2)石家莊地區(qū)平均能見(jiàn)度在正定縣和新樂(lè)縣交界、贊皇縣分別存在2個(gè)高值中心,而在鹿泉市和市區(qū)交界、趙縣分別有2個(gè)低值中心。從春季到冬季平均能見(jiàn)度逐漸降低,高值中心范圍逐漸縮小,低值中心范圍逐漸擴(kuò)大。

    (3)石家莊地區(qū)相對(duì)濕度對(duì)能見(jiàn)度的影響較大,二者存在顯著指數(shù)函數(shù)關(guān)系,而PM2.5和PM10濃度與能見(jiàn)度存在顯著冪函數(shù)關(guān)系。因此,構(gòu)建了能見(jiàn)度與相對(duì)濕度和PM2.5、PM10濃度的多元非線性擬合模型,并獲得不同相對(duì)濕度下模型的擬合系數(shù)。經(jīng)檢驗(yàn),該模型能夠反映石家莊能見(jiàn)度的變化規(guī)律,對(duì)能見(jiàn)度具有一定的預(yù)報(bào)能力。

    [1] 秦世廣,石廣玉,陳林,等. 利用地面水平能見(jiàn)度估算并分析中國(guó)地區(qū)氣溶膠光學(xué)厚度長(zhǎng)期變化特征[J]. 大氣科學(xué),2010,34(2):449-456.

    [2] 白愛(ài)娟,鐘文婷,華蘭,等. 成都市大氣能見(jiàn)度變化特征及影響因子研究[J]. 中國(guó)環(huán)境監(jiān)測(cè),2014,30(2):21-25.

    [3] 周書(shū)華,倪長(zhǎng)健,劉陪川,等. 成都市1980~2010年能見(jiàn)度的變化趨勢(shì)分析[J]. 成都信息工程學(xué)院學(xué)報(bào),2014,29(1):91-96.

    [4] 張浩,石春娥,謝偉,等. 安徽省1955-2005年城市大氣能見(jiàn)度變化趨勢(shì)[J]. 氣象科學(xué),2008,28(5):515-520.

    [5] 范引琪,李春強(qiáng). 1980-2013年京、津、冀地區(qū)大氣能見(jiàn)度變化趨勢(shì)研究[J]. 高原氣象,2008,27(6):1392-1400.

    [6] Doyle M, Dorling S. Visibility trends in the UK 1950-1997[J]. Atmos Environ, 2002,36(19):3161-3172.

    [7] Ghim Y S, Moon K C, Lee S, et al. Visibility trends in Korea during the past two decades[J]. Journal of the Air & Waste Management Association, 2005,55(1):73-82.

    [8] 張智,陳玉華,周紅. 寧夏霧日和霾日的變化趨勢(shì)分析[J]. 干旱氣象,2013,31(4):714-719.

    [9] 江文華,劉德,陳勇航,等. 1980-2012年重慶地區(qū)霾日時(shí)空變化特征[J]. 干旱氣象,2015,33(4):602-606.

    [10] 邵振平. 鄭州機(jī)場(chǎng)能見(jiàn)度變化特征及霧的成因分析[J]. 氣象與環(huán)境科學(xué),2014,37(1):75-82.

    [11] 脫宇峰,潘攀,張麗娟,等. 濟(jì)源市能見(jiàn)度特征及其影響因素分析[J]. 氣象與環(huán)境科學(xué),2012,35(1):56-61.

    [12] 周巧蘭,束炯,王堅(jiān)捍. 上海地區(qū)低能見(jiàn)度特征分析[J]. 氣象與環(huán)境科學(xué),2009,32(2):26-29.

    [13] 劉和平,代佩玲,張青珍,等. 鄭州市大氣能見(jiàn)度變化特征及與空氣污染的關(guān)系[J]. 氣象與環(huán)境科學(xué),2008,31(4):44-46.

    [14] 田密,陳軍,張劍秋. 川西機(jī)場(chǎng)一次平流霧影響能見(jiàn)度過(guò)程的成因分析[J]. 高原山地氣象研究,2013,33(3):77-82.

    [15] 趙清越,劉小渝,李昕翼. 成都市民航氣象行業(yè)標(biāo)準(zhǔn)下的低能見(jiàn)度時(shí)空分布特征研究[J]. 高原山地氣象研究,2014,34(1):57-61.

    [16] 李霞,王勝利. 1980—2007年新疆地區(qū)能見(jiàn)度的變化特征及其影響因子研究[J]. 沙漠與綠洲氣象,2012,6(3):14-20.

    [17] 付桂琴,張迎新,張慶紅,等. 河北省低能見(jiàn)度事件特征分析[J]. 氣象,2013,39(8):1042-1049.

    [18] 劉寧微,馬雁軍,王揚(yáng)鋒. 遼寧中部城市群夏季大氣能見(jiàn)度的觀測(cè)研究[J]. 氣象學(xué)報(bào),2012,70(4):815-820.

    [19] 吳兌,鄧雪嬌,畢雪巖,等. 細(xì)粒子污染形成灰霾天氣導(dǎo)致廣州地區(qū)能見(jiàn)度下降[J]. 熱帶氣象學(xué)報(bào),2007,23(1):1-6.

    [20] 蔡子穎,劉愛(ài)霞,韓素芹,等. 天津低能見(jiàn)度特征初探[J]. 氣象,2014,40(1):114-118.

    [21] 張素敏,王贊紅,張榮英,等. 石家莊市大氣能見(jiàn)度變化特征及其與大氣顆粒物碳成分的關(guān)系[J]. 河北師范大學(xué)學(xué)報(bào),2008,32(6):825-833.

    [22] 李良玉,任愛(ài)玲,梁銀英. 石家莊市顆粒物質(zhì)量濃度及粒度隨高度變化規(guī)律[J]. 環(huán)境科學(xué)與管理,2007,32(1):61-64.

    [23] 王喜全,孫明生,楊婷,等. 京津冀平原地區(qū)灰霾天氣的年代變化[J]. 氣候與環(huán)境研究,2013,18(2):165-170.

    [24] 周靜博,張濤,李治國(guó),等. 利用SPAMS構(gòu)建石家莊市PM2.5固定排放源成分譜庫(kù)[J]. 河北工業(yè)科技,2015,32(5):443-450.

    [25] 戴春嶺,齊堃,康蘇花,等. 石家莊市大氣PM10單顆粒形貌成分分析與來(lái)源識(shí)別[J]. 河北工業(yè)科技,2015,32(5):451-457.

    [26] 牛子寧. 河北省污染物排放的區(qū)域比較研究[D]. 石家莊:河北師范大學(xué).2014.

    [27] 常中陽(yáng),徐振川,許權(quán)貴,等. 石家莊地區(qū)燃煤污染與防治分析[J]. 能源與環(huán)境,2001(7):30-34.

    [28] 宋明,韓素芹,張敏,等. 天津大氣能見(jiàn)度與相對(duì)濕度和PM10及PM2.5的關(guān)系[J]. 氣象與環(huán)境學(xué)報(bào),2013,29(2):34-41.

    [29] 林俊,劉衛(wèi),李燕,等. 大氣氣溶膠粒徑分布特征與氣象條件的相關(guān)性分析[J]. 氣象與環(huán)境學(xué)報(bào),2009,25(1):1-5.

    Variation Characteristics of Atmospheric Visibility and Their Relationship with Relative Humidity and Particle Concentration in Shijiazhuang of Hebei

    WANG Xiaomin, HAN Juncai, CHEN Jing, QIAN Weimiao, YUE Yanxia

    (ShijiazhuangMeteorologicalBureauofHebeiProvince,Shijiazhuang050081,China)

    Based on the visibility of 12 meteorological stations in Shijiazhuang of Hebei Province from 1980 to 2015 and the hourly relative humidity and PM2.5, PM10concentration from 2014 to 2015, the spatial and temporal distributions of visibility were analyzed in Shijiazhuang during 1980-2013, firstly. Then, the relationships between visibility and relative humidity as well as PM2.5, PM10concentration were studied, and on this basis, the multiple nonlinear forecast model of visibility was further established. The results are as follows: (1) The annual average visibility significantly decreased with rate of -1.0 km·(10 a)-1in Shijiazhuang since 1980, and the decreasing trend of average visibility in summer was most obvious, while it was smallest in spring. (2) By comparison of the visibility before and after 1998, the average visibility decreased by 15.3% in Shijiazhuang of Hebei. In space, there were two high centers in the north-central and south-central of Shijiazhuang and two low centers in the urban district of Shijiazhuang and Zhaoxian before 1998, while the visibility gradually decreased from east to west after 1998. (3) The relation of visibility with relative humidity was an exponential function, while the relationships of visibility with PM2.5, PM10concentration were power functions. Therefore, the multiple nonlinear fitting model of visibility was established based on relative humidity, PM2.5and PM10concentration, the model was able to simulate and forecast the visibility in Shijiazhuang to some extent by testing of the observation.

    visibility; relative humidity; particle concentration; temporal and spatial features; nonlinear fitting

    10.11755/j.issn.1006-7639(2016)-04-0648

    2016-01-14;改回日期:2016-04-04

    石家莊市科學(xué)技術(shù)研究與發(fā)展計(jì)劃“重污染天氣前兆‘強(qiáng)信號(hào)’及預(yù)警調(diào)控氣象評(píng)估技術(shù)研究”(151550083A)和石家莊科技計(jì)劃項(xiàng)目“石家莊霧霾天氣監(jiān)測(cè)預(yù)警技術(shù)研究”(131550363A)共同資助

    王曉敏(1986- ),女,碩士,工程師,主要從事環(huán)境氣象預(yù)報(bào)及應(yīng)用方面研究. E-mail:imisstar@hotmail.com

    1006-7639(2016)-04-0648-08DOI:10.11755/j.issn.1006-7639(2016)-04-0648

    P427.2

    A

    王曉敏,韓軍彩,陳靜,等.石家莊地區(qū)能見(jiàn)度變化特征及其與相對(duì)濕度和顆粒物濃度的關(guān)系[J].干旱氣象,2016,34(4):648-655, [WANG Xiaomin, HAN Juncai, CHEN Jing, et al. Variation Characteristics of Atmospheric Visibility and Their Relationship with Relative Humidity and Particle Concentration in Shijiazhuang of Hebei[J]. Journal of Arid Meteorology, 2016, 34(4):648-655],

    猜你喜歡
    能見(jiàn)度實(shí)測(cè)值石家莊
    石家莊曉進(jìn)機(jī)械制造科技有限公司
    肉類研究(2022年7期)2022-08-05 04:47:20
    ±800kV直流輸電工程合成電場(chǎng)夏季實(shí)測(cè)值與預(yù)測(cè)值比對(duì)分析
    常用高溫軸承鋼的高溫硬度實(shí)測(cè)值與計(jì)算值的對(duì)比分析
    哈爾濱軸承(2020年1期)2020-11-03 09:16:22
    2005—2017年白云機(jī)場(chǎng)能見(jiàn)度變化特征及其與影響因子關(guān)系研究
    市售純牛奶和巴氏殺菌乳營(yíng)養(yǎng)成分分析
    一種基于實(shí)測(cè)值理論計(jì)算的導(dǎo)航臺(tái)電磁干擾分析方法
    電子制作(2018年23期)2018-12-26 01:01:22
    低能見(jiàn)度下高速公路主動(dòng)誘導(dǎo)技術(shù)的應(yīng)用
    人民幣緣何誕生在石家莊
    前向散射能見(jiàn)度儀的常見(jiàn)異?,F(xiàn)象處理及日常維護(hù)
    前向散射能見(jiàn)度儀故障實(shí)例分析
    国产午夜精品论理片| 大型黄色视频在线免费观看| 免费电影在线观看免费观看| 丰满人妻一区二区三区视频av| 成人av一区二区三区在线看| 在线看三级毛片| 露出奶头的视频| 国产久久久一区二区三区| 亚洲国产精品成人综合色| 免费av不卡在线播放| 精品人妻一区二区三区麻豆 | 成人高潮视频无遮挡免费网站| 久久精品国产亚洲av天美| 午夜福利在线观看吧| 在线观看一区二区三区| 12—13女人毛片做爰片一| 国产精品亚洲美女久久久| 尾随美女入室| 99视频精品全部免费 在线| 久久草成人影院| 午夜激情欧美在线| 少妇的逼好多水| 亚洲精品一区av在线观看| 国产毛片a区久久久久| 亚洲人成网站高清观看| 国产伦精品一区二区三区四那| 久久久精品欧美日韩精品| 久久久久久伊人网av| 999久久久精品免费观看国产| 中国美白少妇内射xxxbb| 级片在线观看| 国产亚洲欧美98| 美女cb高潮喷水在线观看| 一个人观看的视频www高清免费观看| 午夜a级毛片| 一区二区三区高清视频在线| 亚洲无线观看免费| 午夜久久久久精精品| 亚洲va在线va天堂va国产| 乱人视频在线观看| 日日撸夜夜添| 人妻夜夜爽99麻豆av| 久久精品国产亚洲av涩爱 | 国产欧美日韩精品一区二区| 免费大片18禁| 久久久久久久久久黄片| 精品人妻1区二区| 欧美日韩黄片免| 亚洲午夜理论影院| av黄色大香蕉| 无人区码免费观看不卡| 欧美xxxx性猛交bbbb| 亚洲av日韩精品久久久久久密| 天堂动漫精品| 婷婷六月久久综合丁香| 欧美又色又爽又黄视频| 日韩精品青青久久久久久| 老熟妇乱子伦视频在线观看| 免费人成在线观看视频色| 成年免费大片在线观看| 国模一区二区三区四区视频| av在线天堂中文字幕| 深夜精品福利| 一区福利在线观看| 免费在线观看影片大全网站| 欧美日韩国产亚洲二区| 最近最新中文字幕大全电影3| 如何舔出高潮| 成人性生交大片免费视频hd| 国产高清视频在线观看网站| 黄片wwwwww| 国产真实伦视频高清在线观看 | 久久午夜福利片| 永久网站在线| 国产亚洲精品久久久久久毛片| 亚洲国产精品合色在线| av福利片在线观看| 亚洲中文字幕一区二区三区有码在线看| 不卡视频在线观看欧美| 日韩欧美在线乱码| 一级黄色大片毛片| 大又大粗又爽又黄少妇毛片口| 99热网站在线观看| 国产爱豆传媒在线观看| 国产蜜桃级精品一区二区三区| 18禁在线播放成人免费| 欧美激情在线99| 一夜夜www| 深爱激情五月婷婷| 婷婷色综合大香蕉| 麻豆久久精品国产亚洲av| 国产伦一二天堂av在线观看| 久久久久国产精品人妻aⅴ院| 在线a可以看的网站| 亚洲中文日韩欧美视频| 国产一区二区亚洲精品在线观看| 国产精品人妻久久久影院| 色精品久久人妻99蜜桃| 免费在线观看成人毛片| 日本一本二区三区精品| 国产毛片a区久久久久| 国产高潮美女av| 国产精品免费一区二区三区在线| 免费观看人在逋| 18禁在线播放成人免费| 国产精品爽爽va在线观看网站| 亚洲狠狠婷婷综合久久图片| a级毛片a级免费在线| 少妇的逼好多水| 久久久久久大精品| 免费不卡的大黄色大毛片视频在线观看 | 久久精品影院6| 午夜福利高清视频| 国产精品综合久久久久久久免费| 久久久久久久久久成人| 久久精品综合一区二区三区| 欧美一区二区亚洲| 1024手机看黄色片| 亚洲欧美精品综合久久99| 久久人妻av系列| 美女高潮的动态| 麻豆av噜噜一区二区三区| 老司机福利观看| 又爽又黄无遮挡网站| 两性午夜刺激爽爽歪歪视频在线观看| 免费观看精品视频网站| 欧美高清成人免费视频www| 亚洲无线观看免费| 女同久久另类99精品国产91| 夜夜看夜夜爽夜夜摸| 国产爱豆传媒在线观看| 国产爱豆传媒在线观看| 在线天堂最新版资源| 日本-黄色视频高清免费观看| 成人国产一区最新在线观看| 婷婷六月久久综合丁香| 成人午夜高清在线视频| 1024手机看黄色片| 亚洲美女搞黄在线观看 | 亚洲精品影视一区二区三区av| 日本免费a在线| 制服丝袜大香蕉在线| 天美传媒精品一区二区| 制服丝袜大香蕉在线| 亚洲av一区综合| 我要搜黄色片| 69人妻影院| 亚洲最大成人手机在线| 精品99又大又爽又粗少妇毛片 | 久99久视频精品免费| 中亚洲国语对白在线视频| 亚洲七黄色美女视频| 国产精品嫩草影院av在线观看 | 国产伦精品一区二区三区四那| 中文亚洲av片在线观看爽| 精品人妻偷拍中文字幕| 国产亚洲精品久久久com| 校园人妻丝袜中文字幕| 免费观看的影片在线观看| 国产亚洲av嫩草精品影院| 亚洲久久久久久中文字幕| 中文字幕av在线有码专区| 亚洲成人久久性| 丰满的人妻完整版| 亚洲av二区三区四区| 成人av一区二区三区在线看| 久久中文看片网| 中国美女看黄片| 亚洲人成伊人成综合网2020| 又黄又爽又刺激的免费视频.| 观看免费一级毛片| 国产色婷婷99| 国产高清三级在线| 午夜视频国产福利| 伊人久久精品亚洲午夜| 日韩中文字幕欧美一区二区| 我的老师免费观看完整版| av视频在线观看入口| 欧美激情在线99| 日日干狠狠操夜夜爽| 久久久色成人| 亚洲四区av| 久久99热6这里只有精品| 亚洲欧美精品综合久久99| 一卡2卡三卡四卡精品乱码亚洲| 一级av片app| 日本在线视频免费播放| 国产精品人妻久久久久久| 亚洲国产日韩欧美精品在线观看| 免费看光身美女| 变态另类成人亚洲欧美熟女| 99久久久亚洲精品蜜臀av| 男女做爰动态图高潮gif福利片| 亚洲av日韩精品久久久久久密| 国模一区二区三区四区视频| 国产精品美女特级片免费视频播放器| 成年女人毛片免费观看观看9| 国产综合懂色| 老司机深夜福利视频在线观看| 成人特级黄色片久久久久久久| 久久人人精品亚洲av| 国产精品自产拍在线观看55亚洲| 99久久精品国产国产毛片| 国产久久久一区二区三区| 美女高潮喷水抽搐中文字幕| 亚洲精品粉嫩美女一区| 免费人成视频x8x8入口观看| 国内精品宾馆在线| 国内久久婷婷六月综合欲色啪| 两性午夜刺激爽爽歪歪视频在线观看| 女同久久另类99精品国产91| 蜜桃亚洲精品一区二区三区| 成人鲁丝片一二三区免费| 久久久久久久亚洲中文字幕| 久久久久久久久久成人| 午夜影院日韩av| 最近视频中文字幕2019在线8| 国产成人av教育| 赤兔流量卡办理| av天堂在线播放| 亚洲狠狠婷婷综合久久图片| 午夜激情欧美在线| 免费av观看视频| 日本免费一区二区三区高清不卡| 午夜福利在线在线| 国产91精品成人一区二区三区| 一区二区三区四区激情视频 | 色综合婷婷激情| 国产视频一区二区在线看| 又粗又爽又猛毛片免费看| 日本 av在线| 变态另类丝袜制服| 最近最新免费中文字幕在线| 免费看日本二区| 中文字幕av在线有码专区| 麻豆成人av在线观看| 亚洲专区中文字幕在线| 97超视频在线观看视频| 亚洲黑人精品在线| 一级av片app| 国产精品乱码一区二三区的特点| 午夜免费成人在线视频| 99在线人妻在线中文字幕| 精品欧美国产一区二区三| 国产精品久久久久久久久免| 国产探花在线观看一区二区| 中文亚洲av片在线观看爽| 草草在线视频免费看| 99久久无色码亚洲精品果冻| av黄色大香蕉| 亚洲中文字幕一区二区三区有码在线看| 久久久久久久久大av| 久久久精品欧美日韩精品| 久9热在线精品视频| 波多野结衣巨乳人妻| 久久久午夜欧美精品| 日本爱情动作片www.在线观看 | 国产一级毛片七仙女欲春2| 成人av一区二区三区在线看| 国产精品人妻久久久影院| 久久九九热精品免费| 国产在线精品亚洲第一网站| 亚洲精品一卡2卡三卡4卡5卡| av国产免费在线观看| 日本a在线网址| 日韩欧美在线二视频| 女同久久另类99精品国产91| 美女高潮喷水抽搐中文字幕| 成年女人毛片免费观看观看9| 欧美精品啪啪一区二区三区| 精品久久久久久久久亚洲 | 成人三级黄色视频| 俺也久久电影网| 一级av片app| 国语自产精品视频在线第100页| 男女做爰动态图高潮gif福利片| 春色校园在线视频观看| 欧美高清性xxxxhd video| av.在线天堂| 精品久久久久久久久av| 最近在线观看免费完整版| 国产精品不卡视频一区二区| 欧美又色又爽又黄视频| 亚洲美女黄片视频| 国产av在哪里看| 一夜夜www| 色av中文字幕| 久久精品91蜜桃| 久久久久久国产a免费观看| 麻豆av噜噜一区二区三区| 日本黄色视频三级网站网址| 色综合站精品国产| 精品久久久久久成人av| 一区二区三区高清视频在线| 一级毛片久久久久久久久女| 精品人妻一区二区三区麻豆 | 亚洲av成人精品一区久久| 久久精品人妻少妇| 一个人观看的视频www高清免费观看| 又黄又爽又刺激的免费视频.| 99国产极品粉嫩在线观看| 亚洲精品影视一区二区三区av| 可以在线观看毛片的网站| 搡老岳熟女国产| 日本 欧美在线| 亚洲中文日韩欧美视频| 免费av毛片视频| 国内精品久久久久久久电影| 亚洲av五月六月丁香网| 国产精品,欧美在线| 麻豆成人午夜福利视频| av专区在线播放| 国产高清不卡午夜福利| 亚洲18禁久久av| 午夜久久久久精精品| 亚洲中文日韩欧美视频| 91av网一区二区| 久久久色成人| 成人国产综合亚洲| 十八禁网站免费在线| 91精品国产九色| 欧美最新免费一区二区三区| 亚洲精品粉嫩美女一区| 女的被弄到高潮叫床怎么办 | 夜夜爽天天搞| 亚洲三级黄色毛片| 国产一区二区三区av在线 | 男女下面进入的视频免费午夜| 精品欧美国产一区二区三| 亚洲性久久影院| 99热网站在线观看| 成人av在线播放网站| 日韩精品中文字幕看吧| АⅤ资源中文在线天堂| 精品99又大又爽又粗少妇毛片 | 能在线免费观看的黄片| 欧美激情国产日韩精品一区| 88av欧美| 欧美成人免费av一区二区三区| 亚洲中文字幕日韩| 舔av片在线| 国内久久婷婷六月综合欲色啪| 大型黄色视频在线免费观看| 免费人成视频x8x8入口观看| 日韩一区二区视频免费看| 少妇被粗大猛烈的视频| 精品人妻熟女av久视频| 欧美+日韩+精品| 亚洲av免费在线观看| 成人鲁丝片一二三区免费| 日本成人三级电影网站| 免费不卡的大黄色大毛片视频在线观看 | 中文在线观看免费www的网站| 国产爱豆传媒在线观看| 在线观看午夜福利视频| 中国美白少妇内射xxxbb| 国产伦精品一区二区三区视频9| 一本一本综合久久| av在线亚洲专区| 精品久久久久久久久久免费视频| 国产精品电影一区二区三区| 此物有八面人人有两片| 欧美日韩乱码在线| 午夜久久久久精精品| 午夜福利高清视频| 久久精品影院6| 久久中文看片网| 国产精品一区www在线观看 | 丰满人妻一区二区三区视频av| 亚洲最大成人手机在线| 搡老熟女国产l中国老女人| 美女黄网站色视频| 国产日本99.免费观看| 三级男女做爰猛烈吃奶摸视频| 久久精品综合一区二区三区| 欧美国产日韩亚洲一区| 男女啪啪激烈高潮av片| 嫩草影院入口| 联通29元200g的流量卡| 国产精品无大码| 久久久久国内视频| 久久精品国产亚洲av涩爱 | 久久精品国产鲁丝片午夜精品 | 国产精华一区二区三区| 尾随美女入室| 欧美黑人欧美精品刺激| 免费在线观看成人毛片| 国产精品1区2区在线观看.| 国产一区二区激情短视频| 亚洲精品乱码久久久v下载方式| 国内少妇人妻偷人精品xxx网站| 美女xxoo啪啪120秒动态图| 99热精品在线国产| 18禁黄网站禁片免费观看直播| 一区二区三区免费毛片| 最近最新中文字幕大全电影3| 国产乱人视频| 一区二区三区免费毛片| av天堂中文字幕网| 亚洲中文字幕一区二区三区有码在线看| av.在线天堂| 国内毛片毛片毛片毛片毛片| 久久久久精品国产欧美久久久| 日本与韩国留学比较| 美女cb高潮喷水在线观看| 亚洲欧美日韩东京热| АⅤ资源中文在线天堂| 国产日本99.免费观看| 欧美丝袜亚洲另类 | 日本色播在线视频| 成年人黄色毛片网站| 日韩国内少妇激情av| 露出奶头的视频| 中文字幕av在线有码专区| 亚洲内射少妇av| 老女人水多毛片| 99久久久亚洲精品蜜臀av| 极品教师在线免费播放| 如何舔出高潮| 国产高潮美女av| 性插视频无遮挡在线免费观看| 亚洲成av人片在线播放无| 亚洲av成人av| 精品99又大又爽又粗少妇毛片 | 51国产日韩欧美| 夜夜看夜夜爽夜夜摸| 一区二区三区四区激情视频 | 99精品久久久久人妻精品| 国产老妇女一区| 校园人妻丝袜中文字幕| 亚洲,欧美,日韩| 欧美日韩中文字幕国产精品一区二区三区| 亚洲男人的天堂狠狠| 久久久午夜欧美精品| av女优亚洲男人天堂| 亚洲电影在线观看av| www.www免费av| 国内精品一区二区在线观看| 久久热精品热| 看十八女毛片水多多多| 色哟哟·www| 别揉我奶头 嗯啊视频| 亚洲专区中文字幕在线| 国产精华一区二区三区| 亚洲精品一卡2卡三卡4卡5卡| 欧美成人性av电影在线观看| 亚洲精品一卡2卡三卡4卡5卡| 久久国内精品自在自线图片| 此物有八面人人有两片| 成年免费大片在线观看| 日韩精品青青久久久久久| 欧美色视频一区免费| 亚洲人成网站在线播放欧美日韩| 嫩草影院新地址| 国产精品综合久久久久久久免费| 最后的刺客免费高清国语| 999久久久精品免费观看国产| 午夜爱爱视频在线播放| 久久国产乱子免费精品| 日韩欧美国产在线观看| 久久久久久久午夜电影| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品久久视频播放| 久久久久久九九精品二区国产| 在线看三级毛片| 18禁黄网站禁片免费观看直播| 99久久成人亚洲精品观看| 免费高清视频大片| 91在线观看av| 99久久无色码亚洲精品果冻| 欧美日韩黄片免| 亚洲欧美激情综合另类| 国产主播在线观看一区二区| 大又大粗又爽又黄少妇毛片口| 亚洲性夜色夜夜综合| 白带黄色成豆腐渣| 国内精品宾馆在线| av中文乱码字幕在线| 老司机福利观看| 在线免费观看的www视频| 中文字幕高清在线视频| 色播亚洲综合网| 国内精品美女久久久久久| 国产日本99.免费观看| 成人国产麻豆网| 亚洲在线自拍视频| 听说在线观看完整版免费高清| 成人高潮视频无遮挡免费网站| 97热精品久久久久久| ponron亚洲| 欧美中文日本在线观看视频| 久久99热这里只有精品18| 黄色欧美视频在线观看| 琪琪午夜伦伦电影理论片6080| 国产欧美日韩精品一区二区| 天堂√8在线中文| 国产精品一区二区免费欧美| 男女下面进入的视频免费午夜| 啦啦啦啦在线视频资源| 91久久精品电影网| 亚洲五月天丁香| 亚洲,欧美,日韩| 听说在线观看完整版免费高清| 真人做人爱边吃奶动态| 久久久久久久久大av| 直男gayav资源| 亚洲美女视频黄频| 亚洲美女搞黄在线观看 | 亚洲精品一区av在线观看| 麻豆av噜噜一区二区三区| 亚洲欧美日韩无卡精品| 免费看av在线观看网站| 国产精品久久久久久精品电影| a级一级毛片免费在线观看| 狠狠狠狠99中文字幕| 亚洲中文字幕日韩| 国产亚洲欧美98| 高清日韩中文字幕在线| 在线观看舔阴道视频| 97热精品久久久久久| 97碰自拍视频| 亚洲图色成人| 国产精品日韩av在线免费观看| 少妇高潮的动态图| 能在线免费观看的黄片| 成人国产麻豆网| 99久久精品一区二区三区| 麻豆国产97在线/欧美| 亚洲中文日韩欧美视频| 偷拍熟女少妇极品色| 欧美日韩国产亚洲二区| 制服丝袜大香蕉在线| 熟妇人妻久久中文字幕3abv| 69人妻影院| 国产白丝娇喘喷水9色精品| 国产精品久久久久久亚洲av鲁大| 亚洲黑人精品在线| 看十八女毛片水多多多| 色精品久久人妻99蜜桃| 亚洲一区二区三区色噜噜| 18禁黄网站禁片午夜丰满| 性色avwww在线观看| 日本成人三级电影网站| 日本一二三区视频观看| 午夜福利在线观看免费完整高清在 | 日韩欧美精品免费久久| 永久网站在线| 99精品久久久久人妻精品| 在线观看av片永久免费下载| 国产中年淑女户外野战色| 精品人妻1区二区| 女人被狂操c到高潮| 3wmmmm亚洲av在线观看| 日本三级黄在线观看| 久久精品国产清高在天天线| 老师上课跳d突然被开到最大视频| 精品午夜福利在线看| 成人欧美大片| 91久久精品国产一区二区三区| 国产伦在线观看视频一区| 久久亚洲真实| 日日啪夜夜撸| 又黄又爽又刺激的免费视频.| 制服丝袜大香蕉在线| 国产免费一级a男人的天堂| 精品日产1卡2卡| 琪琪午夜伦伦电影理论片6080| 国产高清不卡午夜福利| 国产视频内射| 亚洲国产欧洲综合997久久,| 成人国产综合亚洲| 热99re8久久精品国产| 色综合亚洲欧美另类图片| 亚洲无线观看免费| 一本一本综合久久| 男女下面进入的视频免费午夜| 久久久久久久久中文| 一本久久中文字幕| 搡老岳熟女国产| 国产精品野战在线观看| 国产精品久久视频播放| 午夜福利18| 啦啦啦韩国在线观看视频| 啪啪无遮挡十八禁网站| 亚洲av中文av极速乱 | 久久久久久久午夜电影| 国产视频一区二区在线看| 免费看日本二区| 亚洲国产欧美人成| 精品不卡国产一区二区三区| 啦啦啦啦在线视频资源| 亚洲内射少妇av| 变态另类成人亚洲欧美熟女| 欧美3d第一页| 女生性感内裤真人,穿戴方法视频| 国产高清不卡午夜福利| 久久精品国产99精品国产亚洲性色| 亚洲国产欧洲综合997久久,| 国产白丝娇喘喷水9色精品| 国产视频内射| 99在线人妻在线中文字幕| .国产精品久久| 最新在线观看一区二区三区| 国产亚洲91精品色在线| .国产精品久久| 欧美中文日本在线观看视频| 久久久久久国产a免费观看| 久久香蕉精品热| 欧美成人免费av一区二区三区| 露出奶头的视频| 身体一侧抽搐| 久久精品国产鲁丝片午夜精品 | 成人特级黄色片久久久久久久| 人妻久久中文字幕网| 日韩强制内射视频|