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

    天津市多模式氣溫集成預(yù)報方法

    2014-07-06 10:54:38吳振玲潘
    應(yīng)用氣象學(xué)報 2014年3期
    關(guān)鍵詞:最低氣溫遺傳算法站點

    吳振玲潘 璇 董 昊 徐 姝 汪 靖

    1)(天津市氣象臺,天津300074)2)(天津科技大學(xué),天津300222)

    天津市多模式氣溫集成預(yù)報方法

    吳振玲1)*潘 璇2)董 昊1)徐 姝1)汪 靖1)

    1)(天津市氣象臺,天津300074)2)(天津科技大學(xué),天津300222)

    在遺傳算法和粒子群算法的基礎(chǔ)上,采用權(quán)重分配方法開展基于混合演化算法的多模式氣溫集成預(yù)報方法研究。利用2012年5—10月中國氣象局GRAPES模式、北京市氣象局BJ-RUC模式、中國氣象局T639模式、天津市氣象局TJWRF模式24h預(yù)報時效的逐6h地面2m高度氣溫和35個天津區(qū)域自動氣象站點資料,通過逐日滾動建立集成預(yù)報模型,對混合演化算法的多模式氣溫集成預(yù)報方法進行了絕對誤差在2℃以內(nèi)的分級、分類及分站檢驗分析。結(jié)果表明:使用該方法建立的氣溫集成預(yù)報模型具有比較可靠的預(yù)報能力,預(yù)報誤差明顯小于任一成員,預(yù)報準確率高。按絕對誤差不大于2℃的檢驗標準,2012年35個站逐6h氣溫、最低氣溫、最高氣溫的集成預(yù)報平均準確率分別為76.34%,77.88%,78.00%。

    混合演化算法;滾動建模;多模式集成;氣溫預(yù)報

    引 言

    20世紀90年初,國外氣象學(xué)者已經(jīng)對數(shù)值天氣預(yù)報模式輸出的預(yù)報結(jié)果開展了集成預(yù)報方法研究,如Leslie等[1]采用統(tǒng)計回歸方法,使用模式預(yù)報結(jié)果和3個統(tǒng)計預(yù)報因子對澳大利亞東南部天氣系統(tǒng)進行集成組合預(yù)報,預(yù)報質(zhì)量顯著提高。隨著國內(nèi)運行的數(shù)值預(yù)報模式種類的迅速增加,面對海量的多模式數(shù)值預(yù)報產(chǎn)品,在短時間內(nèi)有效利用這些預(yù)報信息成為天氣預(yù)報業(yè)務(wù)的主要需求之一,也促使了多模式集成預(yù)報方法在國內(nèi)氣象部門快速發(fā)展。90年代末,多模式集成研究主要集中在降水預(yù)報方面,如金龍等[2]、魏鳳英[3]、周家斌等[4-6]、劉還珠等[7]、陳桂英[8]分別采用神經(jīng)網(wǎng)絡(luò)方法、區(qū)域動態(tài)權(quán)重方法、以汛期降水分布為對象的最小預(yù)測誤差平方和逐步算法、分區(qū)權(quán)重法研究了不同范圍、不同季節(jié)時段的降水集成預(yù)報,預(yù)報準確率優(yōu)于單一集成成員的預(yù)報結(jié)果。在現(xiàn)代天氣業(yè)務(wù)的精細化氣象要素預(yù)報要求下,氣溫集成預(yù)報方法應(yīng)用逐漸增多,如苗愛梅等[9]使用簡單的線性集成方法集成4個模式預(yù)報最低氣溫、最高氣溫;趙聲蓉[10]使用BP人工網(wǎng)絡(luò)神經(jīng)方法集成多種模式產(chǎn)品預(yù)報全國600多個站點的逐3h氣溫,采用同樣方法,李倩等[11]開展區(qū)域850hPa的溫度預(yù)報試驗,集成預(yù)報的平均絕對誤差均小于單模式預(yù)報結(jié)果;熊聰聰?shù)龋?2-13]嘗試分別利用遺傳算法、粒子群算法開展逐3h氣溫的集成預(yù)報試驗;牛保山等[14]利用最小二乘法集成各類上級指導(dǎo)產(chǎn)品,預(yù)報最高氣溫及最低氣溫;林春澤等[15]利用超級集合、多模式集合平均和消除偏差集合平均方法對地面氣溫進行集成預(yù)報。整體上,在精細化氣溫集成預(yù)報研究方面探討和試驗比較多。集成預(yù)報的關(guān)鍵在于構(gòu)造可靠的適應(yīng)度函數(shù),并選取較好的集成預(yù)報方案,以確保其結(jié)果在實際應(yīng)用中較為客觀準確[16],但由于天氣系統(tǒng)的復(fù)雜性,在高密度短期定時氣溫預(yù)報業(yè)務(wù)中還沒有形成比較成熟和完善的集成技術(shù),因此,繼續(xù)加強集成預(yù)報方法的探討仍是一項重要工作。

    天津市氣象局針對GRAPES模式、BJ-RUC模式、T639模式、TJWRF模式的逐6h氣溫預(yù)報產(chǎn)品,開展基于演化算法的多模式氣溫集成預(yù)報方法研究,利用權(quán)重分配方法將遺傳算法和粒子群算法結(jié)合在一起,嘗試在提高數(shù)值預(yù)報產(chǎn)品使用率的同時,進一步提高預(yù)報準確率和預(yù)報效率。

    1 資料與方法

    1.1 資 料

    本文使用2012年5—10月逐日地面氣溫資料,包括實況資料和數(shù)值預(yù)報產(chǎn)品資料。實況資料來自于天津233個區(qū)域自動氣象站中的35個代表站的逐日02:00(北京時,下同),08:00,14:00,20:00氣溫和最低氣溫、最高氣溫。數(shù)值預(yù)報產(chǎn)品選自中國氣象局GRAPES模式、北京市氣象局BJ-RUC模式、中國氣象局T639模式、天津市氣象局TJWRF模式2m高度的逐6h氣溫及最低氣溫、最高氣溫逐日預(yù)報數(shù)據(jù),模式起報時間均為20:00。

    由于中國氣象局GRAPES模式、北京市氣象局BJ-RUC模式、中國氣象局T639模式、天津市氣象局TJWRF模式產(chǎn)品預(yù)報分辨率不同,且輸出場為格點數(shù)據(jù),與天津區(qū)域自動氣象站不一致,因此,利用雙線性插值方法,使用站點周邊4個格點數(shù)據(jù),將格點預(yù)報數(shù)據(jù)轉(zhuǎn)換為區(qū)域站點預(yù)報數(shù)據(jù)。隨機數(shù)據(jù)檢驗表明,轉(zhuǎn)化的站點預(yù)報數(shù)據(jù)與原始格點數(shù)據(jù)具有相同的空間分布及預(yù)報結(jié)論。

    1.2 方 法

    氣溫具有明顯的季節(jié)性特征,尤其是在季節(jié)轉(zhuǎn)換時期,氣溫上下波動幅度相對較大,在保持氣溫連續(xù)性變化特征的基礎(chǔ)上解決季節(jié)影響,采用基于混合演化算法逐日滾動建立集成預(yù)報模型,建模資料長度保持30d。

    集成預(yù)報模型的訂正預(yù)報因子主要選取了上述4個模式的區(qū)域自動氣象站點氣溫預(yù)報數(shù)據(jù)和氣溫實況。

    2 混合演化方法簡介

    混合演化算法(簡稱混合算法)是利用遺傳算法(Genetic Algorithm,簡稱GA)與粒子群優(yōu)化算法(Particle Swarm Optimization,簡稱PSO)共同建立預(yù)報模型的方法。本研究針對同一站點、同一時刻的多組集成訓(xùn)練數(shù)據(jù)序列,即4個模式預(yù)報序列和實況序列,利用遺傳算法和粒子群算法分別建立該站點該時刻的集成預(yù)報模型,完成第1次集成預(yù)報,得到兩組集成預(yù)報數(shù)據(jù)序列,并將這兩組集成預(yù)報數(shù)據(jù)序列轉(zhuǎn)化為混合算法的集成訓(xùn)練數(shù)據(jù)序列,利用式(1)的權(quán)重分配方法,確定每個算法的模型預(yù)報權(quán)重(w),對遺傳算法和粒子群算法預(yù)報結(jié)果進行再集成,得到混合算法的預(yù)報模型(式(2))。

    式(1)中,EGA為遺傳算法模型預(yù)報平均誤差,EPSO為粒子群算法模型預(yù)報平均誤差。

    式(2)中,為兩算法融合后的預(yù)報結(jié)果,wGA和wPSO分別為遺傳算法和粒子群算法的模型預(yù)報權(quán)重,F(xiàn)GA和FPSO分別為遺傳算法和粒子群算法模型預(yù)報結(jié)果。

    2.1 建立氣溫集成預(yù)報模型

    遺傳算法和粒子群算法是在建立數(shù)學(xué)模型的基礎(chǔ)上,通過完成超高維、多局部極值的函數(shù)優(yōu)化計算,完成數(shù)學(xué)模型的求解方法。如果將多模式預(yù)報數(shù)據(jù)看作集成預(yù)報成員數(shù)據(jù),求解得到的結(jié)果就是集成預(yù)報成員的權(quán)重系數(shù)。

    考慮不同模式對不同站點和不同時刻的預(yù)報能力不同,采取了對同一種模式不同站點和不同時刻分別計算權(quán)重系數(shù)的方法,達到集中集合模式優(yōu)勢的目的[14],可建立氣溫集成預(yù)報數(shù)學(xué)模型:

    式(3)中,i為預(yù)報模式(集成預(yù)報成員),j為站點號,t為預(yù)報時刻,F(xiàn)i,j,t為氣溫在站點j第i個模式第t時刻的預(yù)報值,wi,j,t為在站點j第i個模式t時刻的權(quán)重系數(shù)j,t為氣溫站點j第t個時刻的集成預(yù)報值。

    根據(jù)前期的建模試驗,在wi,j,t滿足式(4)的約束條件下,預(yù)報模型最穩(wěn)定:

    2.2 求 解

    將多模式氣溫預(yù)報值看作算法空間中由單元按一定結(jié)構(gòu)組成的數(shù)組或個體,通過初始群體或初始化種群設(shè)定(即式(3)集成預(yù)報數(shù)學(xué)模型權(quán)重系數(shù)初始解的設(shè)定)、適應(yīng)度函數(shù)設(shè)計、遺傳算法或粒子群算法操作設(shè)計和運行參數(shù)設(shè)定等過程,經(jīng)過高次數(shù)的循環(huán)迭代,計算出精細化站點上某時刻不同模式氣溫預(yù)報的權(quán)重系數(shù),將權(quán)重系數(shù)代入數(shù)學(xué)模型,分別建立遺傳算法和粒子群算法的多模式氣溫集成預(yù)報模型。本研究循環(huán)迭代次數(shù)設(shè)定為5000次,遺傳算法求解和粒子群算法求解詳見文獻[12-13]。

    3 氣溫集成預(yù)報應(yīng)用與評估

    根據(jù)日期順序?qū)RAPES模式、BJ-RUC模式、T639模式、TJWRF模式的原始格點預(yù)報數(shù)據(jù)、站點轉(zhuǎn)換預(yù)報數(shù)據(jù)、天津區(qū)域自動氣象站觀測數(shù)據(jù)分別按照預(yù)報間隔時間(6h)錄入數(shù)據(jù)庫,將空缺的數(shù)據(jù)填寫為9999,便于程序檢查。從2012年5月開始分別對天津區(qū)域自動氣象站中的35個代表站點的每日02:00,08:00,14:00,20:00氣溫和最低氣溫、最高氣溫進行試驗預(yù)報。其中,由于模式提供產(chǎn)品的限制,最高氣溫、最低氣溫預(yù)報僅北京市氣象局BJRUC模式和天津市氣象局TJWRF模式參與集成。

    3.1 建模數(shù)據(jù)處理

    采用每日滾動建模方法,基礎(chǔ)數(shù)據(jù)更新的準確性是集成預(yù)報模型質(zhì)量保障。在數(shù)據(jù)更新過程中,用以下方法處理建模數(shù)據(jù):①當某個站點氣溫實況數(shù)據(jù)缺少時,按地理位置選擇距該站點最近的4個站點,以這4個站點的實況平均值替代該站點的數(shù)據(jù),若這4個站點不能滿足要求,則繼續(xù)尋找比較近的站點,直至滿足要求,求取平均值。②當某個模式氣溫預(yù)報數(shù)據(jù)缺少時,按照日期向前尋找存在的記錄,始終保持建模的30條數(shù)據(jù)。③若當日定時氣溫集成模式少于2個,該日氣溫集成預(yù)報模型不做滾動處理。

    3.2 定時氣溫集成預(yù)報分析

    3.2.1 預(yù)報的總體平均絕對誤差對比

    將參與氣溫集成預(yù)報的GRAPES模式、BJRUC模式、T639模式、TJWRF模式稱為集成模式成員。按氣象評分規(guī)則,氣溫誤差在±2℃之內(nèi)的為正確,這里根據(jù)式(5)分別計算4個模式及集成預(yù)報平均絕對誤差,檢驗其氣溫預(yù)報準確性,計算結(jié)果如圖1所示。

    圖1 2012年5—10月集成預(yù)報與GRAPES,BJ-RUC,T639,TJWRF模式預(yù)報總體平均絕對誤差Fig.1 Mean absolute error of temperature from consensus forecast and GRAPES,BJ-RUC,T639,TJWRF forecasts from May to October in 2012

    由圖1可知,無論是定時氣溫預(yù)報,還是最低氣溫、最高氣溫預(yù)報,2012年5—10月集成預(yù)報平均絕對誤差最小,為1.23~1.65℃,完全符合氣象評分正確規(guī)則。表明氣溫集成預(yù)報的準確率高于參與集成模式成員。

    定時氣溫集成預(yù)報的平均絕對誤差與GRAPES模式、BJ-RUC模式、T639模式、TJWRF模式預(yù)報中最大的平均絕對誤差的大小成正比。02:00,08:00,14:00,20:00定時氣溫預(yù)報的最大平均絕對誤差分別為2.96,3.78,7.67,4.25℃,對應(yīng)的集成預(yù)報平均絕對誤差分別為1.23,1.35,1.65,1.58℃。表明參與集成預(yù)報的模式預(yù)報準確率影響集成預(yù)報結(jié)果。

    4個集成模式成員預(yù)報平均絕對誤差不同。其中,定時氣溫及最低氣溫、最高氣溫預(yù)報中,TJWRF模式和BJ-RUC模式的平均絕對誤差相對穩(wěn)定,定時氣溫預(yù)報的平均絕對誤差變化幅度非常小,變化范圍為0.03~0.55℃,其平均絕對誤差也比較小,均在2.63℃之內(nèi)。GRAPES模式預(yù)報的平均絕對誤差為3~4℃,且白天誤差比夜間大。T639模式預(yù)報平均絕對誤差最不穩(wěn)定,距離起報時刻越近,平均絕對誤差越小,14:00,20:00氣溫平均絕對誤差分別超過4℃和7℃。

    3.2.2 預(yù)報穩(wěn)定性對比

    多模式集成預(yù)報結(jié)論直接受集成模式系統(tǒng)穩(wěn)定性影響,如果集成模式成員預(yù)報誤差變化過大,集成預(yù)報誤差也會不穩(wěn)定。由式(5)計算逐月整體平均絕對誤差,圖2為逐月集成預(yù)報與集成模式成員預(yù)報定時氣溫及最低氣溫、最高氣溫的平均絕對誤差對比。表1為集成預(yù)報與集成模式成員預(yù)報2012年5—10月整體平均絕對誤差。通過逐月平均絕對誤差分析,考察4個集成模式成員性能。

    圖2 2012年5—10月集成預(yù)報與GRAPES,BJ-RUC,T639,TJWRF模式預(yù)報氣溫平均絕對誤差Fig.2 Monthly mean absolute error of temperature from consensus forecast and GRAPES,BJ-RUC,T639,TJWRF forecasts

    2012年5—10月,多模式集成預(yù)報定時氣溫和最低氣溫、最高氣溫的平均絕對誤差變化非常小,預(yù)報誤差最小,其平均絕對誤差均在2℃之內(nèi),逐月平均絕對誤差為1.35~1.49℃(表1)。說明利用混合算法每日滾動建立集成預(yù)報模型對各個模式預(yù)報的軌跡描述比較準確,使模型具有預(yù)報穩(wěn)定性。春與夏、夏與秋季節(jié)轉(zhuǎn)換時期幾乎沒有影響多模式集成預(yù)報準確性。

    圖2顯示,盡管GRAPES模式、BJ-RUC模式、T639模式和TJWRF模式的定時氣溫或最低氣溫、最高氣溫預(yù)報平均絕對誤差大小不一致,如TJWRF模式、BJ-RUC模式平均絕對誤差為2~3℃,GRAPES模式為2~5℃,T639模式為2~9℃,但同一時刻氣溫預(yù)報或最低氣溫、最高氣溫預(yù)報的平均絕對誤差大小均具有較好的一致性,如預(yù)報誤差最大的T639模式,02:00,08:00誤差最小,不大于3℃,20:00誤差為4~5℃,14:00誤差最大,為7~9℃。

    表1 2012年5—10月集成預(yù)報與GRAPES,BJ-RUC,T639,TJWRF模式預(yù)報氣溫平均絕對誤差(單位:℃)Table 1 Monthly mean absolute error of temperaturefrom consensus forecast and GRAPES,BJ-RUC,T639,TJWRF forecasts from May to October in 2012(unit:℃)

    3.3 區(qū)域代表站點集成預(yù)報準確率

    按照氣溫預(yù)報評分標準,利用式(6)分別計算2012年5—10月35個站點氣溫集成預(yù)報的平均準確率,隨機、分級檢驗氣溫集成預(yù)報的準確率。即檢驗各站點絕對誤差Ea≤1℃,1℃<Ea≤2℃的預(yù)報準確率。圖3為2012年5—10月35個區(qū)域自動氣象站4次定時氣溫和最高氣溫、最低氣溫預(yù)報的平均準確率。

    圖3 2012年5—10月35個區(qū)域自動氣象站氣溫集成預(yù)報的平均準確率Fig.3 Mean accuracy of temperature at 35automatic meteorological stations from May to October in 2012

    式(6)中,R為預(yù)報準確率,Nr為預(yù)報正確次數(shù),Nf為總預(yù)報次數(shù)。檢驗標準是絕對誤差Ea≤1℃和1℃<Ea≤2℃的預(yù)報準確率。

    圖3顯示,各站氣溫集成預(yù)報絕對誤差Ea≤2℃的準確率基本高于70%。其中,所有站點4次定時氣溫集成預(yù)報的平均準確率達到76.34%,有30個站點的平均準確率為71.75%~82.13%,只有分布在塘沽區(qū)(3個站點)、大港區(qū)(1個站點)、漢沽區(qū)(1個站點)5個站點低于70%。所有站點最低氣溫和最高氣溫集成預(yù)報平均準確率分別達到77.88% 和78.00%,超過2012年人工預(yù)報的天津市13個站點的日最低氣溫、最高氣溫的平均預(yù)報準確率73.42%和74.30%①國家氣象信息中心.2012年1—12月全國城鎮(zhèn)天氣預(yù)報產(chǎn)品質(zhì)量評估報告.2013.。28個站點最低氣溫預(yù)報平均準確率為73.00%~83.33%,28個站點最高氣溫預(yù)報平均準確率為72.67%~86.67%。氣溫集成預(yù)報平均準確率在70%以下的站點多集中在塘沽保稅區(qū)、開發(fā)區(qū)五洲碼頭、大港沙井子、漢沽太平村等地。

    表2為4個定時氣溫集成預(yù)報Ea≤2℃的平均準確率,可以看出,02:00和08:00氣溫平均預(yù)報準確率較高,均超過80%,14:00和20:00的準確率下降,達到了70%以上。

    表2 2012年5—10月定時氣溫集成預(yù)報平均準確率Table 2 Mean accuracy of specified temperature consensus forecast from May to Octorber in 2012

    3.4 氣溫集成預(yù)報對天氣變化的敏感性

    天津北鄰山脈、東鄰渤海、西部與南部與河北平原相接,選取代表站點北部薊縣(54428)、南部靜海(54619)、西部西青(54527)、東部塘沽(54623)和天津(54517)5個站點,2012年利用9—10月各站逐日最低氣溫、最高氣溫集成預(yù)報和實況資料,檢驗氣溫集成預(yù)報對轉(zhuǎn)折性天氣過程的敏感程度。圖4為薊縣、靜海、西青、塘沽和天津5個站點最低氣溫、最高氣溫集成預(yù)報與實況對比。

    由圖4可知,各個單站的最高氣溫、最低氣溫集成預(yù)報與實況變化趨勢一致,表明集成預(yù)報模型對天氣變化的敏感性。最高氣溫集成預(yù)報變化與實況一致,最低氣溫變化略有不同,與平原相鄰西青站和靜海站變化基本一致,北部山區(qū)薊縣站和受熱島效應(yīng)影響最強的天津站,只有1~2次的降溫預(yù)報滯后1d,塘沽站距離渤海較近,其預(yù)報誤差相對較大。說明下墊面性質(zhì)對氣溫集成預(yù)報有一定影響,但通過逐日滾動建立模型技術(shù),可以最大程度地減少這種影響。

    圖4 2012年9—10月最低氣溫、最高氣溫集成預(yù)報與實況對比Fig.4 Comparison between consensus forecast and observation of minimin and maximum temperatures from September to October in 2012

    續(xù)圖4

    3.5 混合演化算法集成預(yù)報的優(yōu)勢

    采用2012年6月資料,將本文介紹按平均誤差計算權(quán)重的混合演化算法與利用遺傳算法、粒子群優(yōu)化算法、等權(quán)重混合算法的集成預(yù)報結(jié)果進行比較分析發(fā)現(xiàn),在相同的迭代計算,預(yù)報絕對誤差不超過2℃的集成預(yù)報準確率檢驗,遺傳算法的集成預(yù)報準確率高于粒子群優(yōu)化算法,按平均誤差計算權(quán)重的混合算法準確率高于按等權(quán)重的混合算法。表3為2012年6月平均最低氣溫集成預(yù)報準確率比較,由表3可知,按平均誤差計算權(quán)重的混合算法的集成預(yù)報準確率最高,比粒子群優(yōu)化算法高3%,比遺傳算法高1%,比等權(quán)重混合算法高3%。

    表3 2012年6月平均最低氣溫集成預(yù)報準確率(單位:%)Table 3 Mean accuracy of minimin temperature from consensus forecast in June 2012(unit:%)

    4 結(jié)論與討論

    1)混合演化算法建模集中了遺傳算法和粒子群優(yōu)化算法的優(yōu)勢。建模過程中,既可以利用遺傳算法比較準確地描述每一個集成模式成員氣溫預(yù)報的變化軌跡,又可以利用粒子群算法優(yōu)化權(quán)重系數(shù)的精確度,縮短計算時間,快速建模,適應(yīng)業(yè)務(wù)應(yīng)用需求。

    2)逐日滾動建立集成預(yù)報模型比較合理。滾動建模可以對實況和模式進行動態(tài)對比分析,并根據(jù)預(yù)報誤差調(diào)整權(quán)重系數(shù),更新預(yù)報模型。

    3)氣溫集成預(yù)報誤差小,可參考性強。定時氣溫及最低氣溫、最高氣溫集成預(yù)報平均絕對誤差均為1.23~1.65℃,明顯好于單一集成模式成員。35個區(qū)域自動氣象站最低氣溫、最高氣溫預(yù)報的絕對誤差不大于2℃的平均準確率超過2012年人工訂正預(yù)報的天津13個站點的24h預(yù)報準確率73.42%和74.30%。

    4)不同下墊面代表站點的逐日定時氣溫、最低氣溫、最高氣溫集成預(yù)報與實況變化基本一致。但下墊面性質(zhì)對最低氣溫集成預(yù)報有影響,尤其是對靠近沿海的自動氣象站影響較大,最低氣溫轉(zhuǎn)折點的預(yù)報滯后1d,且誤差略增。

    定時氣溫集成預(yù)報的平均絕對誤差與4個模式預(yù)報中最大平均絕對誤差大小成正比,本研究未考慮對集成模式成員進行偏差訂正,可能是產(chǎn)生氣溫集成預(yù)報誤差的原因之一。4個模式格點資料插值到站點時未考慮模式地形與測站地形的差異,可能是北部山區(qū)、沿海個別站點氣溫集成預(yù)報誤差較大的原因。另外,2012年5—10月的傍晚和夜間,沿海站點受海風影響較大,模式預(yù)報的天氣系統(tǒng)風向與實況差別較大,可能是導(dǎo)致沿海站點最低氣溫集成預(yù)報準確率下降的原因。

    [1] Leslie L M,Helland G L.Predicting regional forecast skill using single and ensemble forecast technique.Mon Wea Rev,1991,119(2):425-435.

    [2] 金龍,陳寧,林振山.基于人工神經(jīng)網(wǎng)絡(luò)的集成預(yù)報方法研究比較.氣象學(xué)報,1999,57(2):198-207.

    [3] 魏鳳英.全國夏季降水區(qū)域動態(tài)權(quán)重集成預(yù)報試驗.應(yīng)用氣象學(xué)報,1999,10(4):402-409.

    [4] 周家斌,彭瑤.制作汛期降水集成預(yù)報的分區(qū)權(quán)重法.應(yīng)用氣象學(xué)報,1999,10(4):428-435.

    [5] 周家斌,張海福.一種汛期降水分布的客觀集成預(yù)報方法.應(yīng)用氣象學(xué)報,2000,11(增刊):93-97.

    [6] 周家斌,張海福,楊桂英,等.制作汛期降水集成預(yù)報的分區(qū)權(quán)重法.應(yīng)用氣象學(xué)報,1999,10(4):428-435.

    [7] 劉還珠,郝為,林孔元,等.基于智能計算的多模型氣象綜合預(yù)報暴雨落區(qū)預(yù)報實用方法.北京:氣象出版社,2000:30-37.

    [8] 陳桂英.權(quán)重分布法集成預(yù)報試驗.應(yīng)用氣象學(xué)報,2000,11(增刊):51-57.

    [9] 苗愛梅,吳曉荃,董存仙,等.太原地區(qū)溫度集成預(yù)報方法研究.山西氣象,1997(3):2-6.

    [10] 趙聲蓉.多模式溫度集成預(yù)報.應(yīng)用氣象學(xué)報,2006,17(1):52-58.

    [11] 李倩,胡邦輝,王學(xué)忠,等.基于BP人工神經(jīng)網(wǎng)絡(luò)的區(qū)域溫度多模式集成預(yù)報試驗.干旱氣象,2011,29(2):231-235.

    [12] 熊聰聰,王靜,宋鵬,等.遺傳算法在多模式集成天氣預(yù)報中的應(yīng)用.天津科技大學(xué)學(xué)報,2008,23(4):80-84.

    [13] 熊聰聰,許淑祎,徐姝.基于粒子群算法的天氣滾動預(yù)報集成技術(shù).天津科技大學(xué)學(xué)報,2012,27(1):75-78.

    [14] 牛保山,姚繼先,郭文孝,等.縣站溫度集成預(yù)報方法.河南氣象,2004(1):11.

    [15] 林春澤,智協(xié)飛,韓艷,等.基于TIGGE資料的地面氣溫多模式超級集合預(yù)報.應(yīng)用氣象學(xué)報,2009,20(6):706-712.

    [16] 杞明輝,許美玲,程建剛,等.天氣預(yù)報集成技術(shù)和方法應(yīng)用研究.北京:氣象出版社,2006:77-87.

    Forecast Method of Multi-model Air Temperature Consensus in Tianjin

    Wu Zhenling1)Pan Xuan2)Dong Hao1)Xu Shu1)Wang Jing1)
    1)(Tianjin Municipal Meteorological Observatory,Tianjin300074)
    2)(Tianjin University of Science & Technology,Tianjin300222)

    Based on genetic algorithm and particle swarm optimization,multi-model air temperature consensus forecast technology(MMATCFT)of hybrid evolutionary algorithm(HEG)is studied.The main technical thought of this method is that two integrated forecast models are set up respectively by using the genetic algorithm and particle swarm optimization,and then the final mixed forecasting model is established by the weight distribution scheme,which is founded through comparing forecast mean errors between the two models.

    In order to eliminate the impact of seasonal temperature characteristics of Tianjin,the daily rolling integrated forecast model based on 30-day data is adopted in practical operation applications with hybrid evolutionary algorithm.Using 2mair temperature output data of four models of T639,GRAPES,TJWRF,BJ-RUC and observations of 35automatic weather stations(AWS)in villages and towns of Tianjin from May to October in 2012,the forecast test of MMTCFT is carried out.Then the experimentation result is evaluated using the way of classification and station-separation,according to the meteorological standard that absolute error of temperature forecast is within 2℃.T639,GRAPES,TJWRF and BJ-RUC are separately run by China National Meteorological Center,Tianjin Meteorological Bureau and Beijing Meteorological Bureau.The analysis shows that the temperature consensus forecast model is effective and reliable.The technical scheme of the consensus forecast based on rolling model is more rational.The forecast errors are obviously smaller than any model mentioned above,and the forecast accuracy is higher.The average forecast accuracy of 6htemperature,the daily maximum and minimum temperature in 35AWS is 76.34%,77.88%and 78.00%from May to October,respectively.

    genetic algorithm;rolling modeling;multi-model consensus;air temperature forecast

    吳振玲,潘璇,董昊,等.天津市多模式氣溫集成預(yù)報方法.應(yīng)用氣象學(xué)報,2014,25(3):293-301.

    2013-06-13收到,2014-02-27收到再改稿。

    中國氣象局氣象關(guān)鍵技術(shù)集成與應(yīng)用項目(CMAGJ2012M04)

    *email:bigwiselyyn@yahoo.com

    猜你喜歡
    最低氣溫遺傳算法站點
    基于Web站點的SQL注入分析與防范
    電子制作(2019年14期)2019-08-20 05:43:42
    北辰地區(qū)日最低氣溫特征及影響因素研究
    天津科技(2019年3期)2019-03-30 07:17:14
    67年來濟南最高和最低氣溫變化特征
    2017~2018年冬季西北地區(qū)某站點流感流行特征分析
    基于自適應(yīng)遺傳算法的CSAMT一維反演
    一種基于遺傳算法的聚類分析方法在DNA序列比較中的應(yīng)用
    首屆歐洲自行車共享站點協(xié)商會召開
    中國自行車(2017年1期)2017-04-16 02:53:52
    基于遺傳算法和LS-SVM的財務(wù)危機預(yù)測
    怕被人認出
    故事會(2016年21期)2016-11-10 21:15:15
    數(shù)學(xué)選擇題練習(xí)
    国产精品女同一区二区软件 | 国产三级在线视频| 午夜精品久久久久久毛片777| 国产三级在线视频| 免费在线观看影片大全网站| 九色国产91popny在线| 午夜免费激情av| 99热网站在线观看| 亚洲中文字幕日韩| 国产激情偷乱视频一区二区| 午夜爱爱视频在线播放| 久久久久久久精品吃奶| 狠狠狠狠99中文字幕| 国产精品一区二区三区四区久久| 婷婷六月久久综合丁香| 亚洲成人中文字幕在线播放| 一区二区三区高清视频在线| 在线a可以看的网站| 人妻制服诱惑在线中文字幕| 如何舔出高潮| 桃色一区二区三区在线观看| 午夜老司机福利剧场| а√天堂www在线а√下载| 婷婷色综合大香蕉| 黄色女人牲交| 美女大奶头视频| 村上凉子中文字幕在线| 国产精品av视频在线免费观看| 亚洲自偷自拍三级| 国内毛片毛片毛片毛片毛片| 欧美在线一区亚洲| 美女 人体艺术 gogo| 桃色一区二区三区在线观看| 韩国av在线不卡| 一本久久中文字幕| 国产人妻一区二区三区在| 国产精品嫩草影院av在线观看 | 男人和女人高潮做爰伦理| 丰满乱子伦码专区| 搡老妇女老女人老熟妇| 亚洲欧美清纯卡通| 在线看三级毛片| 亚洲图色成人| 99久久精品一区二区三区| 一个人看视频在线观看www免费| 欧美激情国产日韩精品一区| 啦啦啦韩国在线观看视频| 直男gayav资源| 国产久久久一区二区三区| 中出人妻视频一区二区| 女人十人毛片免费观看3o分钟| 亚洲av一区综合| 97超视频在线观看视频| eeuss影院久久| av在线观看视频网站免费| 久久国产乱子免费精品| 88av欧美| av在线天堂中文字幕| 久久午夜福利片| 韩国av在线不卡| 亚洲图色成人| 在线播放无遮挡| 久久香蕉精品热| 国产乱人视频| 永久网站在线| 国产高清不卡午夜福利| 夜夜看夜夜爽夜夜摸| 伦理电影大哥的女人| 成人永久免费在线观看视频| 免费无遮挡裸体视频| 老熟妇仑乱视频hdxx| 搡女人真爽免费视频火全软件 | 欧美日韩国产亚洲二区| 人人妻人人看人人澡| 欧美国产日韩亚洲一区| 大又大粗又爽又黄少妇毛片口| 特级一级黄色大片| 成年女人看的毛片在线观看| 97人妻精品一区二区三区麻豆| 国产av不卡久久| 欧美激情在线99| 色综合站精品国产| 少妇人妻一区二区三区视频| 少妇被粗大猛烈的视频| 亚洲国产精品sss在线观看| 日韩欧美 国产精品| 国产一级毛片七仙女欲春2| 久久人人精品亚洲av| 天堂动漫精品| 69av精品久久久久久| 狂野欧美白嫩少妇大欣赏| 黄色丝袜av网址大全| 国产午夜精品论理片| 999久久久精品免费观看国产| 欧美bdsm另类| 久久久久久久午夜电影| 国产国拍精品亚洲av在线观看| 长腿黑丝高跟| 欧美日韩乱码在线| 真人一进一出gif抽搐免费| 蜜桃亚洲精品一区二区三区| 欧美潮喷喷水| 精品不卡国产一区二区三区| 亚洲精品在线观看二区| 日本精品一区二区三区蜜桃| 色播亚洲综合网| 狂野欧美白嫩少妇大欣赏| 国产大屁股一区二区在线视频| 伊人久久精品亚洲午夜| 亚洲精品亚洲一区二区| 亚洲av.av天堂| 我要搜黄色片| a在线观看视频网站| 成年人黄色毛片网站| 免费看av在线观看网站| 欧美国产日韩亚洲一区| 天堂网av新在线| 白带黄色成豆腐渣| 午夜免费成人在线视频| АⅤ资源中文在线天堂| 欧美成人一区二区免费高清观看| 免费在线观看成人毛片| 两人在一起打扑克的视频| 在线看三级毛片| 日韩欧美国产在线观看| 中国美白少妇内射xxxbb| 内射极品少妇av片p| 国产精品亚洲美女久久久| 国产一级毛片七仙女欲春2| 国产伦精品一区二区三区视频9| 美女xxoo啪啪120秒动态图| 国产欧美日韩一区二区精品| 国产69精品久久久久777片| 成年女人毛片免费观看观看9| 干丝袜人妻中文字幕| 日本在线视频免费播放| netflix在线观看网站| 亚洲精华国产精华液的使用体验 | 男女之事视频高清在线观看| 99九九线精品视频在线观看视频| www.www免费av| 亚洲av不卡在线观看| 国产精品人妻久久久久久| 我的老师免费观看完整版| 色综合亚洲欧美另类图片| 男女边吃奶边做爰视频| 99九九线精品视频在线观看视频| 午夜福利在线观看吧| 国产精华一区二区三区| 日本在线视频免费播放| 国产精品国产高清国产av| 亚洲自拍偷在线| 亚洲18禁久久av| 欧美日韩乱码在线| 男插女下体视频免费在线播放| 桃色一区二区三区在线观看| 国产成人a区在线观看| 亚洲无线观看免费| 亚洲av二区三区四区| 成年女人看的毛片在线观看| 成人国产一区最新在线观看| 黄片wwwwww| 久久精品国产清高在天天线| 国产主播在线观看一区二区| АⅤ资源中文在线天堂| 变态另类丝袜制服| 我要搜黄色片| 亚洲国产精品成人综合色| 好男人在线观看高清免费视频| 欧美色欧美亚洲另类二区| 日本在线视频免费播放| 免费在线观看日本一区| 看十八女毛片水多多多| 中文亚洲av片在线观看爽| 国产午夜精品久久久久久一区二区三区 | 五月伊人婷婷丁香| 色噜噜av男人的天堂激情| 欧美日韩中文字幕国产精品一区二区三区| 中国美女看黄片| 老司机午夜福利在线观看视频| 久久久久精品国产欧美久久久| 午夜视频国产福利| 日韩大尺度精品在线看网址| 看免费成人av毛片| 亚洲美女黄片视频| 哪里可以看免费的av片| 欧美性猛交╳xxx乱大交人| 成人av在线播放网站| 亚洲av.av天堂| 国产一区二区激情短视频| 91久久精品国产一区二区成人| 免费黄网站久久成人精品| 欧美激情久久久久久爽电影| 国产高清不卡午夜福利| 身体一侧抽搐| 国产精品国产三级国产av玫瑰| 免费大片18禁| 亚洲av熟女| 韩国av在线不卡| 老师上课跳d突然被开到最大视频| 可以在线观看毛片的网站| 亚洲精品一卡2卡三卡4卡5卡| 日本 av在线| 黄色视频,在线免费观看| 欧美zozozo另类| 99在线视频只有这里精品首页| 精品人妻一区二区三区麻豆 | 国内精品宾馆在线| 欧美成人免费av一区二区三区| 天天一区二区日本电影三级| 午夜福利高清视频| 国产91精品成人一区二区三区| 99九九线精品视频在线观看视频| 毛片女人毛片| 色av中文字幕| 99久久精品一区二区三区| 精品一区二区三区人妻视频| 搡老岳熟女国产| 99在线人妻在线中文字幕| 99久久精品国产国产毛片| 成人二区视频| 亚洲美女搞黄在线观看 | 成人鲁丝片一二三区免费| 99久久成人亚洲精品观看| 日韩欧美在线乱码| 国产精品一区二区三区四区免费观看 | 欧美成人一区二区免费高清观看| 亚洲最大成人手机在线| 国产精品嫩草影院av在线观看 | 亚洲三级黄色毛片| 欧美+日韩+精品| 两个人的视频大全免费| 国产国拍精品亚洲av在线观看| 久久久国产成人免费| aaaaa片日本免费| 精品久久久久久久久久免费视频| 小说图片视频综合网站| 99久久中文字幕三级久久日本| 精品免费久久久久久久清纯| 老司机福利观看| 欧美最新免费一区二区三区| 我的女老师完整版在线观看| 精品欧美国产一区二区三| 国产又黄又爽又无遮挡在线| av在线亚洲专区| 夜夜看夜夜爽夜夜摸| 亚洲av.av天堂| 在现免费观看毛片| 联通29元200g的流量卡| 亚洲中文字幕日韩| 婷婷色综合大香蕉| or卡值多少钱| 成年人黄色毛片网站| 成人美女网站在线观看视频| 熟女电影av网| 免费高清视频大片| 中文字幕精品亚洲无线码一区| 国产伦一二天堂av在线观看| 精华霜和精华液先用哪个| 啪啪无遮挡十八禁网站| 观看免费一级毛片| 我的老师免费观看完整版| 狂野欧美激情性xxxx在线观看| 亚洲一区高清亚洲精品| 老熟妇乱子伦视频在线观看| 久久精品国产清高在天天线| 成人性生交大片免费视频hd| 国产欧美日韩精品一区二区| 日韩国内少妇激情av| 12—13女人毛片做爰片一| 亚洲av美国av| 色播亚洲综合网| 黄色丝袜av网址大全| 色哟哟哟哟哟哟| 亚洲av日韩精品久久久久久密| 成人鲁丝片一二三区免费| 亚洲午夜理论影院| 日本黄大片高清| 久久亚洲精品不卡| 国产伦人伦偷精品视频| 午夜福利成人在线免费观看| 在线观看66精品国产| 国产精品三级大全| 美女 人体艺术 gogo| 国产成人一区二区在线| 波多野结衣高清无吗| 男女边吃奶边做爰视频| 夜夜夜夜夜久久久久| 亚洲经典国产精华液单| 在线免费观看的www视频| 亚洲成人久久爱视频| 国产精品亚洲美女久久久| av国产免费在线观看| 99九九线精品视频在线观看视频| 一级黄片播放器| 免费在线观看日本一区| 波多野结衣高清无吗| 欧美激情在线99| 欧美最黄视频在线播放免费| 国产精品福利在线免费观看| 亚洲国产精品成人综合色| 欧美极品一区二区三区四区| 国产白丝娇喘喷水9色精品| 国产高清激情床上av| 精品一区二区三区视频在线| 欧美日韩国产亚洲二区| 少妇猛男粗大的猛烈进出视频 | 亚洲自拍偷在线| 国模一区二区三区四区视频| 日韩大尺度精品在线看网址| 日本黄色视频三级网站网址| 九九热线精品视视频播放| 国产精品久久久久久久久免| 免费观看在线日韩| 国产精品精品国产色婷婷| 欧美国产日韩亚洲一区| 日韩精品青青久久久久久| 亚洲七黄色美女视频| 亚洲av成人精品一区久久| 国产高清视频在线播放一区| 国产蜜桃级精品一区二区三区| 夜夜爽天天搞| 男女那种视频在线观看| 国产精品亚洲一级av第二区| 亚洲综合色惰| 精品欧美国产一区二区三| 午夜精品一区二区三区免费看| 国产精品亚洲一级av第二区| 亚洲av成人av| 久久精品国产99精品国产亚洲性色| 欧美日本视频| 18禁黄网站禁片免费观看直播| 午夜免费激情av| 欧美日韩国产亚洲二区| 国产精品久久久久久精品电影| 别揉我奶头~嗯~啊~动态视频| 国产免费男女视频| 88av欧美| 成年女人毛片免费观看观看9| 亚洲电影在线观看av| 一级毛片久久久久久久久女| 又爽又黄a免费视频| 亚洲精品亚洲一区二区| 精品久久久久久久末码| 小蜜桃在线观看免费完整版高清| 一区福利在线观看| 简卡轻食公司| h日本视频在线播放| 国产一区二区亚洲精品在线观看| 99热这里只有精品一区| 99久久久亚洲精品蜜臀av| 亚洲av熟女| 99久久精品一区二区三区| 国产主播在线观看一区二区| 欧美xxxx性猛交bbbb| 国产黄片美女视频| 欧美性猛交黑人性爽| 久久久久精品国产欧美久久久| 婷婷精品国产亚洲av在线| 最后的刺客免费高清国语| 久久久久免费精品人妻一区二区| 国产91精品成人一区二区三区| 蜜桃亚洲精品一区二区三区| 又粗又爽又猛毛片免费看| a在线观看视频网站| 免费av观看视频| 五月玫瑰六月丁香| 91精品国产九色| 两个人的视频大全免费| 亚洲黑人精品在线| 美女被艹到高潮喷水动态| 91麻豆精品激情在线观看国产| 如何舔出高潮| 2021天堂中文幕一二区在线观| 国产真实乱freesex| 国内精品一区二区在线观看| 高清日韩中文字幕在线| 少妇人妻精品综合一区二区 | 啦啦啦韩国在线观看视频| 五月伊人婷婷丁香| 日韩 亚洲 欧美在线| 欧美精品啪啪一区二区三区| 国产成年人精品一区二区| 亚洲人成网站在线播| 久久热精品热| 一个人观看的视频www高清免费观看| 欧美潮喷喷水| 一区福利在线观看| 天天躁日日操中文字幕| 欧美国产日韩亚洲一区| 大型黄色视频在线免费观看| 国产精品女同一区二区软件 | 日韩欧美精品v在线| 国产成人影院久久av| 在线播放无遮挡| 亚洲av美国av| 亚洲欧美日韩高清专用| 少妇的逼水好多| 久久精品国产自在天天线| 国产麻豆成人av免费视频| 在线观看免费视频日本深夜| 精品免费久久久久久久清纯| a级一级毛片免费在线观看| 一进一出好大好爽视频| 久久草成人影院| 又粗又爽又猛毛片免费看| 99久久精品国产国产毛片| 亚洲av五月六月丁香网| 极品教师在线免费播放| 中文字幕人妻熟人妻熟丝袜美| 3wmmmm亚洲av在线观看| 精品欧美国产一区二区三| 日韩欧美在线乱码| 美女免费视频网站| 久久午夜亚洲精品久久| av天堂中文字幕网| 综合色av麻豆| 人妻丰满熟妇av一区二区三区| 亚洲美女黄片视频| 丝袜美腿在线中文| 国产综合懂色| 999久久久精品免费观看国产| 日本免费a在线| 精品一区二区三区av网在线观看| 国产一区二区三区av在线 | 精品一区二区免费观看| 麻豆av噜噜一区二区三区| 日本免费一区二区三区高清不卡| 亚洲精品影视一区二区三区av| 欧美绝顶高潮抽搐喷水| 欧美性感艳星| 天堂av国产一区二区熟女人妻| 深夜精品福利| 精品日产1卡2卡| 亚洲一区二区三区色噜噜| 国产精品精品国产色婷婷| 国产在线精品亚洲第一网站| 色噜噜av男人的天堂激情| 欧美3d第一页| 国产三级中文精品| 亚洲国产色片| 亚洲电影在线观看av| 99国产精品一区二区蜜桃av| 久久久精品大字幕| 人人妻,人人澡人人爽秒播| 最近最新免费中文字幕在线| av在线亚洲专区| 国产成人福利小说| 亚洲乱码一区二区免费版| 日本免费一区二区三区高清不卡| 午夜免费激情av| 99久久九九国产精品国产免费| 男女之事视频高清在线观看| 大又大粗又爽又黄少妇毛片口| 热99在线观看视频| 最新在线观看一区二区三区| 日韩精品青青久久久久久| 成人性生交大片免费视频hd| 内地一区二区视频在线| 在线免费观看的www视频| 久久久国产成人精品二区| 色哟哟哟哟哟哟| 日韩欧美精品免费久久| 国产乱人伦免费视频| av视频在线观看入口| 国产伦在线观看视频一区| 男女边吃奶边做爰视频| 高清在线国产一区| 国产精品久久久久久亚洲av鲁大| 男插女下体视频免费在线播放| 国产白丝娇喘喷水9色精品| 女人十人毛片免费观看3o分钟| 欧美绝顶高潮抽搐喷水| 深爱激情五月婷婷| 18+在线观看网站| 亚洲内射少妇av| 欧美潮喷喷水| 免费高清视频大片| 91麻豆av在线| 日本免费一区二区三区高清不卡| 啪啪无遮挡十八禁网站| 91在线观看av| 久久精品综合一区二区三区| 午夜视频国产福利| 国产av不卡久久| 精品人妻一区二区三区麻豆 | 亚洲在线观看片| 色吧在线观看| 国产高清三级在线| 免费av毛片视频| 亚洲avbb在线观看| 精品国产三级普通话版| 久久精品国产99精品国产亚洲性色| 色哟哟哟哟哟哟| 国产精品久久久久久av不卡| 欧美日韩瑟瑟在线播放| 国产伦人伦偷精品视频| 老司机深夜福利视频在线观看| 亚洲av一区综合| 国产成人一区二区在线| 一级毛片久久久久久久久女| 国产伦在线观看视频一区| 一级毛片久久久久久久久女| 国产免费男女视频| 久久6这里有精品| 黄色配什么色好看| 日日干狠狠操夜夜爽| 少妇熟女aⅴ在线视频| 日本爱情动作片www.在线观看 | 日韩在线高清观看一区二区三区 | 日本 欧美在线| 色综合色国产| 91久久精品国产一区二区成人| 99久久九九国产精品国产免费| 日韩av在线大香蕉| 国产精品1区2区在线观看.| 欧美成人a在线观看| 小说图片视频综合网站| a在线观看视频网站| 欧美三级亚洲精品| 午夜精品在线福利| 1000部很黄的大片| 亚洲成人久久爱视频| 神马国产精品三级电影在线观看| aaaaa片日本免费| 久久人人精品亚洲av| 老司机午夜福利在线观看视频| eeuss影院久久| 日韩欧美在线二视频| 一进一出好大好爽视频| 真实男女啪啪啪动态图| 久久久久久久久久黄片| 亚洲三级黄色毛片| 色综合站精品国产| 韩国av一区二区三区四区| 欧美性感艳星| 国产免费男女视频| 色综合站精品国产| 成人综合一区亚洲| 精品久久国产蜜桃| 国产精品人妻久久久久久| 国产精品女同一区二区软件 | 偷拍熟女少妇极品色| 99精品在免费线老司机午夜| 深夜精品福利| 国产综合懂色| 国产蜜桃级精品一区二区三区| 老师上课跳d突然被开到最大视频| 亚洲综合色惰| 一卡2卡三卡四卡精品乱码亚洲| 国产午夜精品久久久久久一区二区三区 | 全区人妻精品视频| 午夜福利视频1000在线观看| 人妻夜夜爽99麻豆av| 狂野欧美白嫩少妇大欣赏| 午夜影院日韩av| 欧美极品一区二区三区四区| 一进一出好大好爽视频| 黄色欧美视频在线观看| 久久久久久久午夜电影| 亚洲,欧美,日韩| 欧美zozozo另类| 一个人看的www免费观看视频| 久久6这里有精品| 亚洲精品一卡2卡三卡4卡5卡| 欧美成人a在线观看| 欧美不卡视频在线免费观看| 久久久久久久精品吃奶| 免费av不卡在线播放| 成年人黄色毛片网站| 又粗又爽又猛毛片免费看| 久久久国产成人精品二区| 精品国产三级普通话版| 九色国产91popny在线| 国产av不卡久久| av在线亚洲专区| 夜夜爽天天搞| 国产一级毛片七仙女欲春2| av女优亚洲男人天堂| 午夜影院日韩av| xxxwww97欧美| 亚洲成人精品中文字幕电影| 婷婷丁香在线五月| 国产91精品成人一区二区三区| 可以在线观看毛片的网站| 亚洲性夜色夜夜综合| 国产精品一区二区三区四区免费观看 | 国产黄a三级三级三级人| 在线观看舔阴道视频| 99国产精品一区二区蜜桃av| 欧美三级亚洲精品| 非洲黑人性xxxx精品又粗又长| 精品人妻视频免费看| 色尼玛亚洲综合影院| 免费av毛片视频| 亚洲av不卡在线观看| 欧美+亚洲+日韩+国产| 麻豆av噜噜一区二区三区| 午夜a级毛片| 又黄又爽又刺激的免费视频.| 丰满乱子伦码专区| 欧美日韩瑟瑟在线播放| 午夜福利视频1000在线观看| 亚洲av中文av极速乱 | 色噜噜av男人的天堂激情| 国产高清三级在线| 淫秽高清视频在线观看| 亚洲国产精品合色在线| 久久国内精品自在自线图片| 欧美又色又爽又黄视频| 99热这里只有精品一区| 国产欧美日韩精品亚洲av| 91在线观看av| 亚洲精品影视一区二区三区av| 国产免费一级a男人的天堂|