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

    談?wù)劷y(tǒng)計(jì)預(yù)報(bào)的必要性

    2016-04-14 10:49:43陶祖鈺趙翠光陳敏北京大學(xué)物理學(xué)院大氣和海洋科學(xué)系北京0087國家氣象中心北京0008中國氣象局北京城市氣象研究所北京0008
    關(guān)鍵詞:統(tǒng)計(jì)天氣預(yù)報(bào)相關(guān)性

    陶祖鈺 趙翠光 陳敏( 北京大學(xué)物理學(xué)院大氣和海洋科學(xué)系,北京 0087; 國家氣象中心,北京 0008; 中國氣象局北京城市氣象研究所,北京 0008)

    ?

    談?wù)劷y(tǒng)計(jì)預(yù)報(bào)的必要性

    陶祖鈺1趙翠光2陳敏3
    (1 北京大學(xué)物理學(xué)院大氣和海洋科學(xué)系,北京 100871;2 國家氣象中心,北京 100081;3 中國氣象局北京城市氣象研究所,北京 100081)

    摘要:首先從認(rèn)識論的角度,即確定論和概率論,討論了天氣預(yù)報(bào)具有不確定性,在數(shù)值預(yù)報(bào)日益取代主觀經(jīng)驗(yàn)預(yù)報(bào)的今天,統(tǒng)計(jì)預(yù)報(bào)還有其不可取代的重要性。在方法論上,統(tǒng)計(jì)預(yù)報(bào)屬于客觀歸納法;在理論上,統(tǒng)計(jì)預(yù)報(bào)具有不受數(shù)值預(yù)報(bào)可預(yù)報(bào)性限制的特點(diǎn)。統(tǒng)計(jì)預(yù)報(bào)有相關(guān)法和相似法之分;小概率的極端天氣,可以用相似法進(jìn)行預(yù)報(bào)。集合預(yù)報(bào)是確定論方法和概率論方法的有機(jī)結(jié)合,是當(dāng)前天氣預(yù)報(bào)的發(fā)展方向。最后指出,統(tǒng)計(jì)預(yù)報(bào)本身同樣存在局限性,也必須經(jīng)受實(shí)際預(yù)報(bào)應(yīng)用效果的檢驗(yàn)。

    關(guān)鍵詞:天氣預(yù)報(bào),統(tǒng)計(jì),相關(guān)性,相似性,小概率極端天氣

    0 引言

    天氣預(yù)報(bào)有主觀預(yù)報(bào)(或經(jīng)驗(yàn)預(yù)報(bào))和客觀預(yù)報(bào)之分,又有確定性預(yù)報(bào)和不確定預(yù)報(bào)(或概率預(yù)報(bào))之分。數(shù)值模式做出的天氣預(yù)報(bào),既是客觀預(yù)報(bào),又是確定性預(yù)報(bào)。預(yù)報(bào)員的傳統(tǒng)預(yù)報(bào),是根據(jù)所掌握的理論知識和個人實(shí)踐經(jīng)驗(yàn)的主觀預(yù)報(bào),既不定量,又因人因時而變,有很大的不確定性。目前主觀經(jīng)驗(yàn)預(yù)報(bào)在天氣預(yù)報(bào)中的貢獻(xiàn)已在很大程度上被數(shù)值預(yù)報(bào)所取代。

    統(tǒng)計(jì)預(yù)報(bào)也可看作是經(jīng)驗(yàn)預(yù)報(bào),因?yàn)樗彩歉鶕?jù)以往的資料做出的預(yù)報(bào),只不過它避免了人的主觀影響,故也可看作是客觀預(yù)報(bào)。統(tǒng)計(jì)預(yù)報(bào)的結(jié)果雖然是定量的,但它是建立在概率統(tǒng)計(jì)基礎(chǔ)上的預(yù)報(bào),所以本質(zhì)上屬于概率預(yù)報(bào),而不屬于確定性預(yù)報(bào)。

    數(shù)值預(yù)報(bào)雖然被稱作20世紀(jì)人類最偉大的科學(xué)成果之一,并成為當(dāng)前天氣預(yù)報(bào)的基礎(chǔ)和主角,但根據(jù)混沌理論,作為確定性預(yù)報(bào)它有一個不可逾越的極限——2周的可預(yù)報(bào)性[1]。從自然辯證法的哲學(xué)觀點(diǎn)來看,偶然性和必然性都客觀存在,兩者都不可避免,天氣預(yù)報(bào)也不例外。自然辯證法還認(rèn)為,必然性寓于偶然性之中,即必然性通過無數(shù)個個別事例表現(xiàn)出來,只要樣本數(shù)足夠多。因此,統(tǒng)計(jì)預(yù)報(bào)應(yīng)該是天氣預(yù)報(bào)中一個不可缺少的成員。

    早年的統(tǒng)計(jì)預(yù)報(bào),由于觀測資料的年代短,樣本數(shù)不夠多,統(tǒng)計(jì)的對象也很不充分,所以,統(tǒng)計(jì)結(jié)果不穩(wěn)定,預(yù)報(bào)效果當(dāng)然也不理想。這就造成了在過去的數(shù)十年中統(tǒng)計(jì)預(yù)報(bào)的地位既不能與數(shù)值預(yù)報(bào)相比,也不能與主觀經(jīng)驗(yàn)預(yù)報(bào)效果相比的狀況。

    在數(shù)值預(yù)報(bào)逐漸逼近可預(yù)報(bào)性上限的今天,統(tǒng)計(jì)預(yù)報(bào)的不可缺少性也隨之顯現(xiàn)出來。本文將從科學(xué)研究的方法(演繹法和歸納法)及數(shù)值預(yù)報(bào)的現(xiàn)狀和極限談起,介紹相關(guān)法統(tǒng)計(jì)預(yù)報(bào)和相似法統(tǒng)計(jì)預(yù)報(bào)的特點(diǎn),討論統(tǒng)計(jì)預(yù)報(bào)超越確定論預(yù)報(bào)極限的可能性、大概率事件和小概率事件及極端天氣預(yù)報(bào)、集合預(yù)報(bào)是確定論和概率論的結(jié)合等問題,指出統(tǒng)計(jì)預(yù)報(bào)的必要性,強(qiáng)調(diào)集合預(yù)報(bào)是當(dāng)前天氣預(yù)報(bào)技術(shù)的一個重要發(fā)展方向。

    1 演繹法和歸納法

    科學(xué)研究方法分為演繹法和歸納法兩大類。數(shù)值天氣預(yù)報(bào)屬于演繹法,即根據(jù)已知的規(guī)律推演出未來的天氣。演繹法的基礎(chǔ)是分析,因?yàn)樗罁?jù)的規(guī)律是從對個別樣本(或事件)的分析中找出來的。雖然所分析的樣本數(shù)量是有限的或少量的,甚至可以只有一個樣本,但只要所得出的規(guī)律具有普遍性,就可以成為演繹法的基礎(chǔ),例如大家熟知的比薩斜塔自由落體實(shí)驗(yàn),在哲學(xué)上叫作從個別到一般;確定性預(yù)報(bào),則是從一般到個別。

    演繹法的特點(diǎn)是唯一性或確定性,不能模棱兩可。其數(shù)學(xué)形式表現(xiàn)為因變量唯一地由自變量決定。數(shù)值天氣預(yù)報(bào)就是根據(jù)一系列的方程,從初值通過時間積分計(jì)算出未來的天氣。如果沒有計(jì)算誤差,則理論上初值和預(yù)報(bào)兩者是一一對應(yīng)的,只不過演繹的方法不是推理而是數(shù)值積分。

    歸納法是指對已經(jīng)發(fā)生的大量的樣本(或事件)進(jìn)行分類歸納,根據(jù)總結(jié)出的規(guī)律來預(yù)報(bào)未來天氣。經(jīng)典天氣預(yù)報(bào)屬于歸納法,其基礎(chǔ)是分型和概念模型(例如鋒面氣旋模型),根據(jù)當(dāng)前形勢和概念模型的相似性來預(yù)報(bào)未來天氣。尋找歷史個例的相似性是經(jīng)典天氣預(yù)報(bào)(或主觀經(jīng)驗(yàn)預(yù)報(bào))的本質(zhì)[2]。但是主觀經(jīng)驗(yàn)的歸納不僅不夠定量,而且不同的人掌握科學(xué)知識的程度不同,實(shí)踐的經(jīng)驗(yàn)也不同,因此歸納的結(jié)果也不具有唯一性。

    統(tǒng)計(jì),是現(xiàn)代客觀歸納法的數(shù)學(xué)手段。相對于經(jīng)驗(yàn)歸納,統(tǒng)計(jì)是定量的、客觀的歸納。只要統(tǒng)計(jì)模型相同,對于同一組數(shù)據(jù)只有一個統(tǒng)計(jì)結(jié)果,所以現(xiàn)代的客觀歸納法具有一定程度的唯一性。當(dāng)然,如果模型或數(shù)據(jù)不完全相同,顯然,結(jié)果就不具有唯一性。

    統(tǒng)計(jì)歸納法本身雖然不受物理規(guī)律控制(沒有預(yù)設(shè)的控制規(guī)律),但如果統(tǒng)計(jì)方法足夠全面(或充分),統(tǒng)計(jì)樣本足夠大,則統(tǒng)計(jì)結(jié)果可以揭示客觀存在的物理規(guī)律。這方面比較成功的典型例子是Wallace 等[3]的“遙相關(guān)”研究[4]。他對北半球冬季500hPa月平均高度場并不事先設(shè)定關(guān)鍵區(qū),而是普遍地計(jì)算每一個格點(diǎn)與其他全部格點(diǎn)的相關(guān)系數(shù)。從最大正、負(fù)相關(guān)系數(shù)的空間分布?xì)w納出五種遙相關(guān)型,其中緯向的PNA型、EU型和EA型分別反映了太平洋—北美、歐洲—亞洲和東大西洋這三條主要的長波頻散路徑(圖1),具有非常清晰的天氣動力學(xué)意義[5-6],在天氣預(yù)報(bào)和氣候診斷的領(lǐng)域中得到了普遍的肯定和應(yīng)用。

    圖1 Wallace的遙相關(guān)型和長波頻散路徑[3](閉合圈為相關(guān)系數(shù)大值中心,實(shí)線和J分別為200hPa急流軸和急流中心)Fig.1 Wallace teleconnection patterns and long-wave dispersion path(closed loop is correlation coefficientmaximum center,Vector lines and J were 200 hPa jet axis and the jet center)

    2 數(shù)值預(yù)報(bào)的現(xiàn)狀和極限

    圖2 歐洲中期數(shù)值預(yù)報(bào)中心1981—2011年全球和南北半球500hPa高度場預(yù)報(bào)3,5,7和10天距平相關(guān)分布(來源:Uccellini LW美國宇航局/大氣海洋局衛(wèi)星資料同化聯(lián)合中心成立十周年紀(jì)念會報(bào)告幻燈片)Fig.2 The anomaly correlation coefficient of 500 hPa height forecast field in global and the northern and southern hemisphere at 3 days,5 days,7 days,10 days,Medium-Range Forecast Centre in Europe in 1981—2011

    從歐洲中期數(shù)值預(yù)報(bào)中心(ECMWF)對1981年以來數(shù)值預(yù)報(bào)水平的評估(圖2)可以看到,500hPa高度場3天預(yù)報(bào)的距平相關(guān)目前已高達(dá)97%,而且進(jìn)入21世紀(jì)后,3天預(yù)報(bào)相關(guān)系已接近水平線,看不出明顯的上升趨勢。它似乎表明,500hPa環(huán)流形勢的3天預(yù)報(bào)已達(dá)到幾乎完美的程度,很難再有所提高。其他數(shù)值模式與ECMWF相比沒有本質(zhì)上的差別,預(yù)報(bào)水平也大體相當(dāng)。但是,從降水和溫度等的單站要素的強(qiáng)度、時間的確定性預(yù)報(bào)來看,無論哪個模式都還有很多不能令人滿意之處。即使是看似相對比較容易的2m溫度預(yù)報(bào),也有相當(dāng)大的預(yù)報(bào)誤差(±2~4℃)。所以在預(yù)報(bào)溫度、降水等要素時,必須采用數(shù)值預(yù)報(bào)產(chǎn)品的統(tǒng)計(jì)釋用(如MOS方法(模式輸出統(tǒng)計(jì))和PP方法(完全預(yù)報(bào))等)。

    圖3 2012年7月21日06,09,12和15時(UTC)的觀測小時降水(單位:mm)及各起報(bào)時刻預(yù)報(bào)對應(yīng)時刻的降水分布[8](第1—6行分別表示CMORPH降水融合資料起報(bào)時間分別為20日21時,21日00,03,06,09和12時,第7行為21日15時實(shí)況;各列分別表示21日06,09,12和15時的降水分布)Fig.3 In July 21,2012 0600,0900,1200 and 1500/UTC observation hour precipitation(mm)and the corresponding time precipitation forecast map(1-6th row denote the reported time of CMORPH precipitation fusion data is at 2100 on the 20th,at 0000,0300,0600,0900 and 1200 on the 21th,the 7th row is real data at 1500 on the 21th; the columns represent precipitation distribution at 0000,0900,1200 and 1500 on the 21th)

    對于降水或強(qiáng)對流天氣的預(yù)報(bào)則更是難以令人滿意。例如2012年7月21日北京的特大暴雨過程,世界各國的模式大都在2~3天前就預(yù)報(bào)了這次強(qiáng)降水過程,但降水強(qiáng)度普遍偏小,降水時間普遍偏晚約6小時,都沒有報(bào)出最強(qiáng)的第一場鋒前暖區(qū)暴雨[7-8]。從圖3所給的北京快速循環(huán)中尺度數(shù)值模式(BJ-RUC)的預(yù)報(bào)檢驗(yàn)看[8],即使是同一個數(shù)值模式,不同初始時刻的預(yù)報(bào)結(jié)果也有明顯差別。更有甚者,預(yù)報(bào)水平并不一定隨著時效的縮短而變好,表現(xiàn)出時好時壞的現(xiàn)象。這種搖擺不定就是數(shù)值預(yù)報(bào)“不確定性”的典型表現(xiàn)。它說明雖然模式控制方程所代表的物理規(guī)律是確定的,但實(shí)際計(jì)算結(jié)果卻存在不確定性,因?yàn)槌踔岛兔恳徊讲罘钟?jì)算都不是100%的正確,總有一些誤差,即使非常微小。所以,不確定性是不可避免的,確定性預(yù)報(bào)的有效時間長度是有限的。

    在圖2中還可看到,5,7和10天的相關(guān)系數(shù)隨著預(yù)報(bào)時效的延長而下降,10天預(yù)報(bào)的相關(guān)系數(shù)已下降到不足0.5,作為確定性預(yù)報(bào)似可認(rèn)為已無參考價值。歐洲中期預(yù)報(bào)中心30年的總結(jié)說明,確定性預(yù)報(bào)的水平隨時效延長而下降是不可避免的,數(shù)值預(yù)報(bào)的實(shí)際能力不可能突破2周可預(yù)報(bào)性的理論極限。后面將會提到,要突破這個極限的辦法是采取確定論和概率論相結(jié)合的方法,即數(shù)值模式的集合預(yù)報(bào)。

    3 統(tǒng)計(jì)預(yù)報(bào)——超越確定性預(yù)報(bào)極限

    統(tǒng)計(jì)預(yù)報(bào)的性質(zhì)屬于歸納法,統(tǒng)計(jì)是歸納的數(shù)學(xué)手段。統(tǒng)計(jì)預(yù)報(bào)可分為相關(guān)法和相似法兩大類;相似法又可分為直接相似法和分型相似法兩類。

    3.1相關(guān)法統(tǒng)計(jì)預(yù)報(bào)

    相關(guān)法是用預(yù)報(bào)因子和預(yù)報(bào)量之間的直接相關(guān)性來做預(yù)報(bào),例如建立多因子回歸方程。在建立預(yù)報(bào)方程前,通常先要進(jìn)行因子篩選,即先根據(jù)經(jīng)驗(yàn)和科學(xué)知識選取一定數(shù)量的因子從中挑選出相關(guān)系數(shù)最大的作為預(yù)報(bào)因子。所以,統(tǒng)計(jì)方法雖然是定量的、客觀的,但也存在因人而異的問題。因?yàn)榭梢宰鳛轭A(yù)報(bào)因子的參數(shù)很多,選擇不同的因子,就會得到不同的統(tǒng)計(jì)結(jié)果(相對而言,采用不同的統(tǒng)計(jì)方法造成的差異則要小得多)。在長期天氣預(yù)報(bào)會商中經(jīng)常遇到這種情況,不同的統(tǒng)計(jì)預(yù)報(bào)得到不同的結(jié)果,因而難下預(yù)報(bào)決心。所以。統(tǒng)計(jì)預(yù)報(bào)也存在主觀性。

    理論上,如果因子的篩選能夠窮盡一切可能的因子,同時樣本數(shù)又足夠大,就可以克服統(tǒng)計(jì)預(yù)報(bào)的主觀性。當(dāng)前計(jì)算能力已強(qiáng)大到足已從海量的因子(如數(shù)百、數(shù)千甚至數(shù)萬個因子)中篩選預(yù)報(bào)因子,無需人為地圈定待選因子的范圍,那么,統(tǒng)計(jì)預(yù)報(bào)就有可能完全擺脫主觀預(yù)設(shè)的影響,而只反映客觀事實(shí)的內(nèi)在聯(lián)系。

    3.2相似法統(tǒng)計(jì)預(yù)報(bào)

    相似法統(tǒng)計(jì)預(yù)報(bào)又可分為兩類:直接相似預(yù)報(bào)和分型相似預(yù)報(bào)[2,9-10]。

    3.2.1直接相似預(yù)報(bào)

    如文獻(xiàn)[11]所述,直接相似預(yù)報(bào)就是“從眾多的歷史個例中找出一個與當(dāng)前形勢相似程度非常高的個例,則可以認(rèn)為引起當(dāng)前形勢演變的全部機(jī)制都基本上隱含在這個歷史個例中,因此可以直接根據(jù)歷史相似個例的演變作為當(dāng)前預(yù)報(bào)的基礎(chǔ)”,故稱為直接相似預(yù)報(bào)。在人們對天氣形勢演變和局地天氣產(chǎn)生的物理機(jī)制還不甚了解的時代,直接相似法是天氣預(yù)報(bào)的基本方法,故郭樹忠等[2]把它稱為經(jīng)典預(yù)報(bào)方法。他們定義了一個所謂相似系數(shù)(S),

    用它在眾多的歷史樣本中找出最相似的個例來制作當(dāng)前機(jī)場天空狀況的逐時預(yù)報(bào)。式中,Xi為當(dāng)前某一氣象要素的測站值,下標(biāo)i代表測站,Yi為同一要素的某一歷史個例,N為測站的總數(shù),和為全部測站的平均值,即

    理論上任意一個氣象要素場都可以用傅氏級數(shù)展開成諧波的集合。兩個波長相同的諧波,其相似性應(yīng)包括振幅、位相和其圍繞的平均值三個方面。可以證明[11],經(jīng)典的相關(guān)系數(shù)(C),

    由于對變量都做了歸一化處理,所以只能反映諧波的位相相似程度,同位相時相關(guān)系數(shù)為1,反位相時相關(guān)系數(shù)為-1。而S不僅反映位相的相似,還能反映振幅的相似程度。為了使S還能同時反映平均值的相似程度,陶祖鈺等[11]對S進(jìn)行了修改,即

    但是,如文獻(xiàn)[11]所指出的,“天氣形勢的相似應(yīng)包括天氣系統(tǒng)的位置、相對強(qiáng)度和絕對強(qiáng)度三個方面。相對強(qiáng)度相似是指高值系統(tǒng)和低值系統(tǒng)之間差異大小的相似。氣壓系統(tǒng)和溫度系統(tǒng)相對強(qiáng)度的大小與鋒區(qū)、急流等有直接的聯(lián)系,因此相對強(qiáng)度相似無疑是十分重要的。絕對強(qiáng)度相似是指天氣系統(tǒng)本身的數(shù)值是否接近。絕對強(qiáng)度和各年的氣候特征有關(guān)。例如,副熱帶高壓(以下簡稱副高)的強(qiáng)度和各年副高季節(jié)變化的早晚有關(guān),它是天氣形勢是否相似的一個重要的方面”。因此,試圖用一個相似系數(shù)來代表二維場的相似無疑是難以做到的。

    3.2.2分型相似預(yù)報(bào)

    天氣分型是經(jīng)典天氣預(yù)報(bào)中最常用的方法。分型相似法和相關(guān)法不同之處在于,前者不僅僅看到氣象變量本身,還看到變量的空間和時間分布,然后再按當(dāng)前實(shí)況所屬的類型或者最相似的歷史個例做出預(yù)報(bào)。天氣形勢的主觀經(jīng)驗(yàn)分型種類繁多,有的簡約,有的繁瑣;有的著眼高空,有的著眼低空或地面,如此等等不一而足。而且,分型的標(biāo)準(zhǔn)往往模糊不清,同一個樣本,不同的人會確定為不同的型。

    現(xiàn)代統(tǒng)計(jì)學(xué)中的聚類分析是一種可以客觀地按相似性進(jìn)行分類的方法。例如芬蘭Teuvo Kohonen于1981年提出的自組織特征映射(Self-organizing Feature Map,SOM)神經(jīng)網(wǎng)絡(luò)法[12]是一種通過反復(fù)調(diào)整識別圖像之間的相似程度(包括數(shù)值的相似和形狀的相似)從而達(dá)到聚類的目的。趙翠光用這一技術(shù)對2000—2009年6—8月華北地區(qū)107次區(qū)域暴雨前一天的環(huán)流形勢場(500hPa高度場、850hPa流場、850hPa溫濕場(θse))進(jìn)行了客觀分型[13]。從圖4給出的其中兩種型的降水分布和500hPa高度場形勢(省略850hPa形勢)可見,客觀分型結(jié)果的天氣學(xué)意義非常明顯。圖4a和4b為華北經(jīng)向型暴雨,其特點(diǎn)是副高偏北、河套低槽明顯;圖4c和4d為華北的緯向暴雨型,其特點(diǎn)是副高偏東、華北受低槽影響。

    現(xiàn)代統(tǒng)計(jì)學(xué)中分型的方法很多,例如EOF(Empirical Orthogonal Function)方法,即經(jīng)驗(yàn)正交函數(shù)分解方法,是一種既包括空間分布特征,也包括時間變化特征的客觀分型方法。文獻(xiàn)[13]中用旋轉(zhuǎn)主分量(REOF)方法[14-16]對華北地區(qū)夏季降水進(jìn)行分區(qū),客觀地反映出華北范圍內(nèi)不同區(qū)域降水的氣候特征。

    圖4 華北夏季區(qū)域性暴雨的客觀分型[13],副高偏北型暴雨典型場(河套低槽類)日平均暴雨實(shí)況(a,單位:mm)和500hPa高度場(b,單位:dagpm),副高偏東型暴雨典型場(華北低槽類)日平均暴雨實(shí)況(c)和500hPa高度場(d)Fig.4 Objective classification of Regional Rainstorm in North China in summer(13),Daily average observed precipitation about northerly subtropical rainstorm type(trough in Hetao)(a,mm)and 500hPa height(b,dagpm),daily average observed precipitation about easterly subtropical rainstorm type(trough in North China)(c,mm)and 500 hPa height(d,dagpm)

    圖像識別也是統(tǒng)計(jì)學(xué)應(yīng)用中的一個重要領(lǐng)域。廣義地看,尋找相似也屬于圖像識別領(lǐng)域。光流法是一種研究外界物體運(yùn)動的光影變化在人眼視網(wǎng)膜上所形成的光的流動的圖像識別技術(shù)[17]。它特別適用于圖像雷達(dá)回波系統(tǒng)的識別和追蹤,因此是在短時和臨近預(yù)報(bào)中很有應(yīng)用價值的一種統(tǒng)計(jì)學(xué)方法[18]。

    4 討論

    4.1預(yù)報(bào)時效

    統(tǒng)計(jì)預(yù)報(bào)的特點(diǎn)是不預(yù)設(shè)物理?xiàng)l件。無論時間間隔多長,只要前后兩組數(shù)據(jù)之間表現(xiàn)出有足夠的相關(guān)性,就可以把前一組數(shù)據(jù)當(dāng)作預(yù)報(bào)因子來做預(yù)報(bào),因此它沒有時效的限制。數(shù)值預(yù)報(bào)時效的理論極限是兩周,實(shí)際上超過7天高度距平、相關(guān)系數(shù)小于0.5表明預(yù)報(bào)接近隨機(jī)概率,參考價值不大[7]。因此10天以上的中長期預(yù)報(bào),在相當(dāng)長的時間內(nèi)還只能依靠統(tǒng)計(jì)預(yù)報(bào)。

    直接相似法統(tǒng)計(jì)預(yù)報(bào)的特點(diǎn)是只看結(jié)果,不問過程,不講“道理”。例如相似法統(tǒng)計(jì)預(yù)報(bào),只要現(xiàn)在相似,則未來也一定相似,無論3天、5天,還是10天、半月,甚至更長的時效,對相似法而言并無差異(注:當(dāng)然原則上隨著預(yù)報(bào)時效的延長相似的空間范圍應(yīng)該擴(kuò)大到全球)。雖然我們對統(tǒng)計(jì)出來的結(jié)果(如分型),可以做出一定的物理解釋,但統(tǒng)計(jì)過程并不受物理規(guī)律控制,因此可以說統(tǒng)計(jì)預(yù)報(bào)是不講“道理”的。這是統(tǒng)計(jì)預(yù)報(bào)和確定論的數(shù)值預(yù)報(bào)的最大不同。正因?yàn)榇耍y(tǒng)計(jì)預(yù)報(bào)可以不受數(shù)值預(yù)報(bào)時效的可預(yù)報(bào)性限制??傊?,從當(dāng)前數(shù)值預(yù)報(bào)的水平看,超過一周的延伸期預(yù)報(bào)就不能直接依靠數(shù)值預(yù)報(bào),而應(yīng)該采用統(tǒng)計(jì)預(yù)報(bào)方法來制作,只要樣本的數(shù)量足夠多。

    4.2轉(zhuǎn)折性天氣預(yù)報(bào)

    從氣象服務(wù)的角度看,轉(zhuǎn)折性天氣,如寒潮爆發(fā)、干旱、連續(xù)低溫陰雨的開始和結(jié)束,或夏季高溫伏旱的開始和結(jié)束等都是社會大眾最關(guān)心的,因此是天氣預(yù)報(bào)的重點(diǎn)。從預(yù)報(bào)的難度來看,連續(xù)數(shù)日不變的穩(wěn)定天氣最容易預(yù)報(bào);但是轉(zhuǎn)折性天氣,一旦啟動,變化迅速,導(dǎo)致預(yù)報(bào)失誤。此外,對于重大的天氣轉(zhuǎn)折,預(yù)報(bào)員很難下決心當(dāng)機(jī)立斷,也是造成預(yù)報(bào)失誤或影響預(yù)報(bào)服務(wù)效果的心理原因,而統(tǒng)計(jì)預(yù)報(bào)的客觀性可以為預(yù)報(bào)員下決心提供支持。因此,統(tǒng)計(jì)預(yù)報(bào)是一種比較適合轉(zhuǎn)折性天氣的預(yù)報(bào)方法。

    4.3極端天氣預(yù)報(bào)——大概率事件和小概率事件

    從統(tǒng)計(jì)學(xué)的角度看,天氣事件可簡單地分為兩大類:大概率事件和小概率事件。大概率事件的樣本多,統(tǒng)計(jì)關(guān)系穩(wěn)定,比較適用于相關(guān)法的統(tǒng)計(jì)預(yù)報(bào),例如晴雨預(yù)報(bào)、溫度預(yù)報(bào)等。

    極端天氣,如暴雨或大暴雨等,屬于小概率事件,由于樣本很少,無法建立穩(wěn)定的統(tǒng)計(jì)關(guān)系,因此不能直接使用相關(guān)性做預(yù)報(bào)。但是,有長期工作經(jīng)驗(yàn)的預(yù)報(bào)員仍可以根據(jù)親身經(jīng)歷的案例做出適當(dāng)?shù)念A(yù)報(bào),只要?dú)v史上曾經(jīng)有過類似的案例。

    相似法是對全體樣本進(jìn)行分類,既包括大概率事件,也包括小概率事件,所以相似法原則上也可用于極端天氣預(yù)報(bào)。例如,2012年7月21日北京的特大暴雨(簡稱“7?21”暴雨)對北京地區(qū)的雨量紀(jì)錄而言是百年,甚至是二百年一遇的極端暴雨事件。圖5出示了2012年7月21日08時的海平面氣壓場。按Zhao等[19]傾斜旋轉(zhuǎn)T模態(tài)主成分分析方法對1951—2012年夏季5704個海平面氣壓場客觀分型的結(jié)果看(圖6),“7?21”暴雨的環(huán)流形勢屬于暴雨概率最大的第六型,但其氣候概率也只有0.491%(圖6中環(huán)流型概率和暴雨概率的乘積),仍屬于小概率事件。從環(huán)流型看,此型海上副熱帶高壓偏強(qiáng),暖濕空氣與河套地區(qū)的冷空氣在華北對峙,是華北比較典型的經(jīng)向型暴雨形勢,在趙翠光[13]的分型中也同樣可以看到這種類型(圖4a和4b)。

    圖5 2012年7月21日08時的海平面氣壓場(單位:hPa)Fig.5 July 21,2013 0800 SLP(unit:hPa)

    與此相反的是第二型,環(huán)流特征是副高偏南,北方都在冷氣團(tuán)的控制下,所以出現(xiàn)此型暴雨的氣候概率更小,只有0.052%,可以說是真正的罕見暴雨。從分型結(jié)果中找出對這類罕見暴雨的歷史個例進(jìn)行分析,或許可為此類罕見暴雨的預(yù)報(bào)提供重要的線索。

    上述海平面氣壓場的客觀分型盡管分型時并未考慮降水的情況,但分型的結(jié)果卻能反映出環(huán)流形勢與暴雨之間的物理聯(lián)系。由此可見,只要樣本數(shù)足夠大,統(tǒng)計(jì)結(jié)果中就會反映出事物內(nèi)在的規(guī)律性。在自然辯證法中這就叫做必然性存在于偶然性之中(或普遍性存在于特殊性之中)。同時,即使是“小概率”事件,它也受普遍的內(nèi)在規(guī)律控制。像2008年金融危機(jī)和“9?11”事件,它們的發(fā)生也都有前兆可循[20],極端天氣預(yù)報(bào)也不例外。統(tǒng)計(jì)的本質(zhì)就是將所需要的信息從大量的、充滿噪音的信息中提煉出來。在信息爆炸時代(大數(shù)據(jù)時代)特別需要統(tǒng)計(jì)學(xué)。

    4.4集合預(yù)報(bào)——確定論和概率論的結(jié)合

    如前所述,數(shù)值預(yù)報(bào)屬于確定性預(yù)報(bào),但是數(shù)值模式預(yù)報(bào)的氣象要素(如溫度、風(fēng)和降水)存在相當(dāng)大的誤差。為了提高預(yù)報(bào)的準(zhǔn)確率,人們提出數(shù)值預(yù)報(bào)產(chǎn)品釋用方法,如模式輸出統(tǒng)計(jì)方法MOS(model output statistics)。它將模式輸出的物理量作為預(yù)報(bào)因子,通過建立與預(yù)報(bào)量之間的統(tǒng)計(jì)關(guān)系(如逐步回歸方程)來完成預(yù)報(bào)。所以MOS預(yù)報(bào)是一種確定論和概率論相結(jié)合的預(yù)報(bào)方法,是一種間接使用數(shù)值預(yù)報(bào)的方法。

    圖6 1951—2012年夏季9種海平面氣壓場客觀分型[19]、發(fā)生頻率(左上數(shù)字)與此型發(fā)生暴雨的頻率(右上數(shù)字)Fig.6 Objective classification of nine kinds SLP field in summer 1951-2012r[19],frequency of occurrence(upper left number)and the frequency of this type rainstorm occurs(upper right number)

    當(dāng)前數(shù)值預(yù)報(bào)和統(tǒng)計(jì)預(yù)報(bào)相結(jié)合的方法是集合預(yù)報(bào)。所謂“集合”,就是用數(shù)值模式制造出數(shù)十個乃至上百個不同的預(yù)報(bào)結(jié)果。其中的每一個成員,都與一定的模式過程和初值相對應(yīng),都是確定的預(yù)報(bào)。但是由于每個成員的模式過程或初值多存在一些微小的差異,所以預(yù)報(bào)結(jié)果各不相同,整個集合的預(yù)報(bào)有一定的離散度。換言之,集合預(yù)報(bào)是用確定論方法人為制造出來的不確定性。然后,再用統(tǒng)計(jì)方法從這個預(yù)報(bào)集合的離散特征中做出預(yù)報(bào)。集合預(yù)報(bào)不同于MOS預(yù)報(bào),仍是直接的數(shù)值預(yù)報(bào),統(tǒng)計(jì)只是歸納數(shù)值預(yù)報(bào)結(jié)果的方法。更重要的是,集合預(yù)報(bào)從方案設(shè)計(jì)開始就在概率論的指導(dǎo)之下進(jìn)行,所以集合預(yù)報(bào)是確定論和概率論方法的有機(jī)結(jié)合。顯然,集合預(yù)報(bào)的思路更符合客觀世界變化的根本規(guī)律。

    特別有意思的是,集合預(yù)報(bào)不僅可以提供大概率事件的預(yù)報(bào),也可提供極端事件的預(yù)報(bào)。圖7是從互聯(lián)網(wǎng)下載的各國數(shù)值模式對721暴雨過程的集合預(yù)報(bào)比較[7]。圖中給出的不是具體的降水量,而是24h極端降水的發(fā)生概率。雖然對圖中極端降水的定義及其發(fā)生概率的具體算法還沒有深入了解,但圖7所顯示的預(yù)報(bào)效果令人驚嘆。因此,開展集合預(yù)報(bào)的研究和應(yīng)用,在廣大預(yù)報(bào)員中普及集合預(yù)報(bào)的基本知識是很有必要的。

    本文第一作者早年曾做過少量的統(tǒng)計(jì)預(yù)報(bào)工作,所用的也都是比較簡單的統(tǒng)計(jì)方法。本文只是通過對統(tǒng)計(jì)預(yù)報(bào)基本概念的討論和幾個統(tǒng)計(jì)預(yù)報(bào)具體事例的介紹,試圖說明統(tǒng)計(jì)預(yù)報(bào)是天氣預(yù)報(bào)方法中不可缺失一個方面,值得引起大家的重視。

    篇后語:初稿完成后,作者又把《信號與噪音——為什么有的預(yù)報(bào)失敗,但有的成功》[20]和《赤裸裸的統(tǒng)計(jì)學(xué)——除去大數(shù)據(jù)的枯燥外衣,呈現(xiàn)真實(shí)的數(shù)據(jù)之美》[21]這兩本書找出來從頭至尾重新翻閱。感到除了強(qiáng)調(diào)統(tǒng)計(jì)預(yù)報(bào)的必要性,還必須指出統(tǒng)計(jì)預(yù)報(bào)的局限性,還有一些非常重要的話需要說,例如統(tǒng)計(jì)預(yù)報(bào)是否真的可以不講“道理”。為此補(bǔ)充介紹以下內(nèi)容。

    這兩本書都一方面指出“在大數(shù)據(jù)時代,統(tǒng)計(jì)的重要性怎么強(qiáng)調(diào)都不過分”,但同時也強(qiáng)調(diào)統(tǒng)計(jì)數(shù)字也會“撒謊”。如果數(shù)據(jù)不足,或者統(tǒng)計(jì)方法錯誤也會得出謬以千里的結(jié)論。書中舉了一個淺顯易懂但值得每一個預(yù)報(bào)員深思的例子:統(tǒng)計(jì)數(shù)據(jù)可能得出“工作時小憩的人更易死于癌癥”這樣的結(jié)論,但實(shí)際是離開辦公室休息的員工中有很多人在辦公樓入口處抽煙,是抽煙而非小憩引發(fā)了癌癥。可見,如果在數(shù)據(jù)處理過程中有先入為主的偏見,就可能犯重大的錯誤。只有對科學(xué)保持謙虛謹(jǐn)慎的態(tài)度,承認(rèn)客觀世界具有不可測性和不確定性,才能避免錯誤。

    圖7 各國模式2012年7月18日20時的第3~4天的24h雨量極端天氣概率預(yù)報(bào)[7](來源:http://tparc.mri-jma.go.jp/TIGGE/tigge_extreme_prob.html)Fig.7 Extreme weather probability forecast for 24h rainfall of 3-4th day at 2000 July 18th,in each country mode[7]

    書中舉了瑞典數(shù)學(xué)家、作家安德烈斯的一句話:用數(shù)據(jù)說謊容易,但是用數(shù)據(jù)說出真相卻很難。因此書的作者認(rèn)為最重要的是首先必須明確自己所持的觀點(diǎn)源于哪個科學(xué)原理,避免曲解統(tǒng)計(jì)結(jié)果。也就是說統(tǒng)計(jì)預(yù)報(bào)也必須講“道理”,不可盲目。

    書中還有一個重要的觀點(diǎn)是:“假設(shè)不能算是科學(xué),假設(shè)只有通過預(yù)測得到檢驗(yàn)才能成為科學(xué)?!币簿褪钦f,中長期預(yù)報(bào)和氣候研究中的統(tǒng)計(jì)結(jié)論,都要經(jīng)過預(yù)報(bào)檢驗(yàn)才能成立??傊?,“實(shí)踐才是檢驗(yàn)真理的唯一標(biāo)準(zhǔn)”。

    參考文獻(xiàn)

    [1]丑紀(jì)范.天氣和氣候的可預(yù)報(bào)性.氣象科技進(jìn)展,2011,1(2):11-14.

    [2]郭樹忠,羅寶玉.經(jīng)典預(yù)報(bào)自動化系統(tǒng).航空氣象,1989,(增刊1):11-13.

    [3]Wallace J M,Gultzler D S.Teleconections in the geopotential height field during the Northern Hemisphere winter.Mon Wea Rev,1981,109(4):784-812.

    [4]陶祖鈺,郭煒.冬季500 hPa侯平均高度場的遙相關(guān).氣象學(xué)報(bào),1989,47(4):486-492.

    [5]Hoskins B J,Karoly D J.The steady linear response of a spherical atmosphere to thermal and orographic forcing.J Atmos Sci,1981,38:1179-1196.

    [6]Simmons A J,Wallace J M,Branstator G M.Barotropic wave propagation and instability and atmospheric teleconnection patterns.J Atmos Sci,1983,40:1363-1392.

    [7]陶祖鈺,鄭永光.“7·21”北京特大暴雨的預(yù)報(bào)問題.暴雨災(zāi)害,2013,32(3):193-206.

    [8]姜曉曼,袁慧玲,等.北京“7.21”特大暴雨高分辨率模式分析場及預(yù)報(bào)分析.氣象學(xué)報(bào),2014,72(2):207-219.

    [9]李南聲.用網(wǎng)格點(diǎn)相關(guān)系數(shù)法選相似形勢場.氣象,1982,8(1):16-17.

    [10]趙鵬福,邱志生.環(huán)流型聚類分析問題的初步探討.氣象,1985,11(6):13-16.

    [11]陶祖鈺,湯桂生.尋找相似天氣形勢的客觀方法研究.大氣科學(xué)研究與應(yīng)用,1993,(5):81-88.

    [12]Kohonen T .Self-Org anization Maps.Berlin:Springer Verlag,1995.

    [13]趙翠光.華北夏季降水多尺度研究及客觀預(yù)報(bào).南京:南京信息工程大學(xué),2012.

    [14]施能.氣象統(tǒng)計(jì)預(yù)報(bào).北京:氣象出版社,2009.

    [15]黃嘉佑.氣象統(tǒng)計(jì)分析與預(yù)報(bào)方法.北京:氣象出版社,1990.

    [16]魏鳳英.現(xiàn)代氣候統(tǒng)計(jì)診斷與預(yù)測技術(shù).北京:氣象出版社,2007.

    [17]柳士俊,張蕾.光流法及其在氣象領(lǐng)域里的應(yīng)用.氣象科技進(jìn)展,2015,5(4):16-21.

    [18]韓雷,王洪慶,林隱靜.光流法在強(qiáng)對流天氣臨近預(yù)報(bào)中的應(yīng)用.北京大學(xué)學(xué)報(bào)(自然科學(xué)版),2008,44(5):751-755.

    [19]Zhao Y Y,Zhang Q H,Du Y,et al.Objective analysis of circulation extremes during the 21 July 2012 torrential rain in Beijing.Acta Meteor Sin,2013,27(5):626-635。

    [20]納特·希爾弗.信號與噪音.胡曉嬌,張新,朱辰辰,譯.北京:中信出版社,2013.

    [21]查爾斯·惠倫.赤裸裸的統(tǒng)計(jì)學(xué).曹檳,譯.北京:中信出版社,2013.

    The Necessity of Statistical Forecasts

    Tao Zuyu1,Zhao Cuiguang2,Chen Min3
    (1 Department of Atmospheric and Ocean Sciences,School of Physics,Peking University,Beijing 100871 2 National Meteorological Centre,Beijing 1000813 Institute of Urban Meteorology,China Meteorological Administration,Beijing 100081)

    Abstract:First,weather forecast uncertainty is discussed from the perspective of epistemology,namely,deterministic and probabilistic theory.Today,subjective forecast is increasingly replaced by the numerical weather prediction(NWP),but statistical forecasting is still irreplaceable.In terms of methodology,statistical forecasting is objective induction.Statistical prediction is not limited by predictability of NWP in theory.Statistical forecasting has a correlation or similar method; rare event related with extreme weather can be forecasted in similar method.Ensemble prediction method is a combination of deterministic and probabilistic methods,and is the development direction of the current weather forecast.Finally,there are limitations on statistical prediction itself,and statistical forecasting must be also verified by the effect in the forecast practice.

    Keywords:weather forecasts,statistics,correlation,similarity,rare event of extreme weather

    收稿日期:2014年3月3日;修回日期:2014年6月5日

    DOI:10.3969/j.issn.2095-1973.2016.01.001

    第一作者:陶祖鈺(1939—),Email:taozuyu@pku.edu.cn

    資助信息:公益性行業(yè)(氣象)科研專項(xiàng)(GYHY200906003,GYHY201206004,GYHY201306023);國家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(2013CB430106)

    猜你喜歡
    統(tǒng)計(jì)天氣預(yù)報(bào)相關(guān)性
    天氣預(yù)報(bào)員
    天氣預(yù)報(bào)的前世今生
    中期天氣預(yù)報(bào)
    2008—2015我國健美操科研論文的統(tǒng)計(jì)與分析
    人間(2016年27期)2016-11-11 17:33:19
    山東省交通運(yùn)輸投資計(jì)劃管理信息系統(tǒng)的設(shè)計(jì)
    我國居民收入與消費(fèi)關(guān)系的統(tǒng)計(jì)分析
    基于Kronecker信道的MIMO系統(tǒng)通信性能分析
    科技視界(2016年21期)2016-10-17 17:37:34
    市場經(jīng)濟(jì)背景下的會計(jì)統(tǒng)計(jì)發(fā)展探究
    商(2016年27期)2016-10-17 05:50:19
    小兒支氣管哮喘與小兒肺炎支原體感染相關(guān)性分析
    腦梗死與高同型半胱氨酸的相關(guān)性研究(2)
    国产黄频视频在线观看| 人成视频在线观看免费观看| 精品国产超薄肉色丝袜足j| 啦啦啦视频在线资源免费观看| 青春草亚洲视频在线观看| 久久久精品94久久精品| 最近中文字幕2019免费版| 亚洲av在线观看美女高潮| 天堂俺去俺来也www色官网| 日本欧美国产在线视频| 久久综合国产亚洲精品| 看非洲黑人一级黄片| 1024视频免费在线观看| 国产有黄有色有爽视频| 免费在线观看完整版高清| 青春草视频在线免费观看| 久久久精品区二区三区| 亚洲综合精品二区| 亚洲精品av麻豆狂野| 欧美日韩亚洲国产一区二区在线观看 | 男人操女人黄网站| 欧美精品一区二区大全| 你懂的网址亚洲精品在线观看| 九草在线视频观看| 国产毛片在线视频| 青春草视频在线免费观看| 人人妻,人人澡人人爽秒播 | 丰满乱子伦码专区| 最新在线观看一区二区三区 | 久久久欧美国产精品| 天天影视国产精品| 亚洲少妇的诱惑av| 国产亚洲精品第一综合不卡| 久久久久精品性色| 国产精品嫩草影院av在线观看| 黄片播放在线免费| 久久久久精品国产欧美久久久 | 一级毛片我不卡| 亚洲情色 制服丝袜| 久久综合国产亚洲精品| 精品第一国产精品| 免费av中文字幕在线| 精品国产乱码久久久久久小说| 亚洲精品自拍成人| 最近最新中文字幕大全免费视频 | 国产免费现黄频在线看| 日韩熟女老妇一区二区性免费视频| 亚洲在久久综合| 亚洲国产精品999| 日韩欧美一区视频在线观看| 精品亚洲成国产av| 免费黄网站久久成人精品| 国产在线一区二区三区精| 亚洲专区中文字幕在线 | 美女脱内裤让男人舔精品视频| 国产av码专区亚洲av| 国产熟女欧美一区二区| 亚洲欧美精品自产自拍| 最近最新中文字幕大全免费视频 | 中文字幕另类日韩欧美亚洲嫩草| 天天添夜夜摸| av.在线天堂| 中文字幕人妻丝袜一区二区 | 满18在线观看网站| 午夜免费观看性视频| tube8黄色片| 久久人人97超碰香蕉20202| 蜜桃国产av成人99| 欧美 日韩 精品 国产| 国产无遮挡羞羞视频在线观看| 在线免费观看不下载黄p国产| 久久久久精品性色| 亚洲国产精品成人久久小说| 国产 精品1| 丝袜脚勾引网站| 少妇人妻精品综合一区二区| 久久久久精品性色| 亚洲av日韩精品久久久久久密 | 久久人人爽人人片av| 最近2019中文字幕mv第一页| 水蜜桃什么品种好| 国产乱来视频区| 蜜桃国产av成人99| 亚洲专区中文字幕在线 | 成人手机av| 国产 一区精品| 国产国语露脸激情在线看| 国产麻豆69| 中国三级夫妇交换| 亚洲欧洲日产国产| 制服人妻中文乱码| 9色porny在线观看| 观看美女的网站| 中文字幕av电影在线播放| 天天躁日日躁夜夜躁夜夜| 久久久久精品性色| 亚洲av综合色区一区| 中文字幕亚洲精品专区| 久久久久久人妻| 日韩制服丝袜自拍偷拍| 国产深夜福利视频在线观看| 日韩中文字幕欧美一区二区 | 曰老女人黄片| 777久久人妻少妇嫩草av网站| 一二三四在线观看免费中文在| 热re99久久精品国产66热6| 女人爽到高潮嗷嗷叫在线视频| 亚洲,一卡二卡三卡| 无限看片的www在线观看| 国产精品久久久久成人av| 亚洲欧美一区二区三区国产| 国产精品一二三区在线看| 婷婷色av中文字幕| 亚洲少妇的诱惑av| 久久婷婷青草| 欧美黄色片欧美黄色片| 久久鲁丝午夜福利片| 你懂的网址亚洲精品在线观看| 久久99精品国语久久久| 日本色播在线视频| 少妇被粗大猛烈的视频| 精品国产一区二区久久| 久久久久视频综合| 久久久久精品国产欧美久久久 | 色婷婷久久久亚洲欧美| 久久人人爽人人片av| 欧美日韩亚洲高清精品| videos熟女内射| 午夜福利一区二区在线看| 黄片小视频在线播放| 日本黄色日本黄色录像| 亚洲人成77777在线视频| 中文字幕人妻丝袜制服| 久久天躁狠狠躁夜夜2o2o | 久久99精品国语久久久| 在线观看免费午夜福利视频| 中文精品一卡2卡3卡4更新| 国产激情久久老熟女| www.精华液| 中文字幕色久视频| 精品久久久久久电影网| 亚洲国产欧美在线一区| 欧美国产精品一级二级三级| 一区二区三区四区激情视频| 老司机影院成人| 一区二区三区激情视频| 久久精品国产a三级三级三级| 免费在线观看完整版高清| 黄频高清免费视频| 观看av在线不卡| 国产高清国产精品国产三级| 亚洲综合精品二区| 男人舔女人的私密视频| 国产淫语在线视频| 欧美日韩综合久久久久久| 最新的欧美精品一区二区| 精品久久久久久电影网| 青青草视频在线视频观看| 99九九在线精品视频| 麻豆精品久久久久久蜜桃| 天天影视国产精品| 免费在线观看完整版高清| 国产男女超爽视频在线观看| 亚洲精品日韩在线中文字幕| 桃花免费在线播放| 黑人欧美特级aaaaaa片| 热99国产精品久久久久久7| 久久久久久免费高清国产稀缺| 久久久国产欧美日韩av| 免费看不卡的av| 久久人妻熟女aⅴ| 欧美亚洲 丝袜 人妻 在线| 日韩伦理黄色片| 多毛熟女@视频| 丰满少妇做爰视频| 亚洲欧美一区二区三区国产| 日本91视频免费播放| 亚洲欧美成人精品一区二区| 色播在线永久视频| 中文字幕高清在线视频| 最近的中文字幕免费完整| 如何舔出高潮| 日韩视频在线欧美| 日韩,欧美,国产一区二区三区| 久久免费观看电影| 国产97色在线日韩免费| 成人18禁高潮啪啪吃奶动态图| 女人爽到高潮嗷嗷叫在线视频| 久久久久久人人人人人| 亚洲国产看品久久| 捣出白浆h1v1| 一级片免费观看大全| 日本一区二区免费在线视频| av网站在线播放免费| xxxhd国产人妻xxx| 久热爱精品视频在线9| 日日啪夜夜爽| 欧美黄色片欧美黄色片| 久久免费观看电影| 精品一区二区三区四区五区乱码 | 中文字幕另类日韩欧美亚洲嫩草| 国产精品国产三级专区第一集| 最新在线观看一区二区三区 | 欧美日韩成人在线一区二区| 日韩不卡一区二区三区视频在线| 国产激情久久老熟女| 午夜91福利影院| 蜜桃国产av成人99| 一级爰片在线观看| 欧美黑人欧美精品刺激| 亚洲精品日韩在线中文字幕| xxx大片免费视频| 久久久久精品国产欧美久久久 | 欧美 日韩 精品 国产| 日韩电影二区| 亚洲成色77777| 久久久久精品久久久久真实原创| 国产男女内射视频| 黄色视频在线播放观看不卡| 日日啪夜夜爽| 大香蕉久久网| 中文欧美无线码| 又粗又硬又长又爽又黄的视频| 啦啦啦在线免费观看视频4| 秋霞在线观看毛片| 国产色婷婷99| 国产免费福利视频在线观看| 少妇人妻久久综合中文| 91精品伊人久久大香线蕉| 国产无遮挡羞羞视频在线观看| 国产一区二区三区综合在线观看| 在线观看免费日韩欧美大片| 大片电影免费在线观看免费| 美国免费a级毛片| 一区二区三区激情视频| 男女无遮挡免费网站观看| 免费女性裸体啪啪无遮挡网站| 99香蕉大伊视频| 在线观看国产h片| 一级a爱视频在线免费观看| av在线老鸭窝| 国产成人精品在线电影| 亚洲美女黄色视频免费看| 亚洲国产中文字幕在线视频| 精品卡一卡二卡四卡免费| 丰满饥渴人妻一区二区三| 男人添女人高潮全过程视频| 51午夜福利影视在线观看| 久久久精品国产亚洲av高清涩受| 精品少妇久久久久久888优播| 一区二区三区激情视频| avwww免费| 高清黄色对白视频在线免费看| 女人精品久久久久毛片| 在线观看免费高清a一片| 久久久精品区二区三区| 美女中出高潮动态图| 日韩人妻精品一区2区三区| 中文精品一卡2卡3卡4更新| 日本欧美国产在线视频| 亚洲第一av免费看| 亚洲av日韩精品久久久久久密 | 99精国产麻豆久久婷婷| av福利片在线| 啦啦啦在线免费观看视频4| 热re99久久精品国产66热6| 人妻一区二区av| 两性夫妻黄色片| 亚洲av电影在线观看一区二区三区| 下体分泌物呈黄色| 久久精品熟女亚洲av麻豆精品| 90打野战视频偷拍视频| bbb黄色大片| 卡戴珊不雅视频在线播放| 蜜桃国产av成人99| 深夜精品福利| 大香蕉久久成人网| 国产成人啪精品午夜网站| 在现免费观看毛片| 一个人免费看片子| 免费观看人在逋| 狂野欧美激情性xxxx| 久久精品亚洲熟妇少妇任你| 一级爰片在线观看| 美女脱内裤让男人舔精品视频| 国产极品粉嫩免费观看在线| 久久人人爽av亚洲精品天堂| 亚洲精品在线美女| 深夜精品福利| 亚洲天堂av无毛| 丁香六月欧美| 美女视频免费永久观看网站| 久久鲁丝午夜福利片| 超色免费av| 久久久久人妻精品一区果冻| 美女脱内裤让男人舔精品视频| 日本一区二区免费在线视频| 久久97久久精品| 黑人巨大精品欧美一区二区蜜桃| 亚洲男人天堂网一区| 日本av免费视频播放| 最近中文字幕高清免费大全6| 一级黄片播放器| 少妇被粗大的猛进出69影院| 男女边摸边吃奶| svipshipincom国产片| 久热爱精品视频在线9| 亚洲国产精品999| 18在线观看网站| 欧美日韩视频高清一区二区三区二| 国产乱来视频区| 亚洲精品美女久久av网站| 一区在线观看完整版| 久久青草综合色| 久久人人97超碰香蕉20202| 亚洲精品美女久久久久99蜜臀 | 国产精品 国内视频| www.熟女人妻精品国产| 免费黄色在线免费观看| 国产不卡av网站在线观看| a级毛片在线看网站| 夫妻午夜视频| 亚洲精品第二区| 人人妻人人爽人人添夜夜欢视频| 日韩大码丰满熟妇| 成人漫画全彩无遮挡| 久久久久久免费高清国产稀缺| 久久久久久久久久久免费av| 色婷婷av一区二区三区视频| 超碰成人久久| 精品一品国产午夜福利视频| 十八禁人妻一区二区| 日本av手机在线免费观看| 97在线人人人人妻| 亚洲欧美一区二区三区黑人| 亚洲精品国产一区二区精华液| 亚洲欧美中文字幕日韩二区| 中文字幕人妻丝袜一区二区 | 美女国产高潮福利片在线看| 久久久欧美国产精品| 91国产中文字幕| 日韩成人av中文字幕在线观看| 丝瓜视频免费看黄片| 国产av码专区亚洲av| 一级片'在线观看视频| 久久久久国产精品人妻一区二区| 18禁国产床啪视频网站| 99热国产这里只有精品6| 91精品伊人久久大香线蕉| 精品亚洲成国产av| 亚洲精品aⅴ在线观看| 久久热在线av| 999久久久国产精品视频| 人人妻人人澡人人爽人人夜夜| 自拍欧美九色日韩亚洲蝌蚪91| 黄色视频在线播放观看不卡| 国产成人一区二区在线| av免费观看日本| 亚洲国产最新在线播放| 丝袜喷水一区| 天天躁夜夜躁狠狠久久av| 国产一区二区 视频在线| 黑人巨大精品欧美一区二区蜜桃| 看十八女毛片水多多多| 亚洲av男天堂| 下体分泌物呈黄色| 国产不卡av网站在线观看| 亚洲一卡2卡3卡4卡5卡精品中文| 热re99久久国产66热| 搡老岳熟女国产| 亚洲美女黄色视频免费看| av卡一久久| 日韩精品免费视频一区二区三区| 一区二区三区精品91| 91精品伊人久久大香线蕉| 少妇精品久久久久久久| a级毛片在线看网站| 观看美女的网站| videos熟女内射| 国产女主播在线喷水免费视频网站| 亚洲精品国产色婷婷电影| 晚上一个人看的免费电影| 免费观看人在逋| 免费av中文字幕在线| 国产乱来视频区| 成人亚洲欧美一区二区av| 中国三级夫妇交换| 亚洲欧美精品综合一区二区三区| 99国产精品免费福利视频| 99久久精品国产亚洲精品| 亚洲一区中文字幕在线| 人妻一区二区av| 精品国产一区二区久久| 欧美国产精品va在线观看不卡| 欧美老熟妇乱子伦牲交| 哪个播放器可以免费观看大片| 午夜日韩欧美国产| 秋霞伦理黄片| 日韩伦理黄色片| a 毛片基地| 1024香蕉在线观看| 久久99精品国语久久久| 两个人看的免费小视频| 老熟女久久久| 日韩一本色道免费dvd| 天天躁日日躁夜夜躁夜夜| 80岁老熟妇乱子伦牲交| 久久久精品国产亚洲av高清涩受| 亚洲欧美一区二区三区国产| 看免费成人av毛片| 狂野欧美激情性xxxx| 亚洲欧美日韩另类电影网站| 性高湖久久久久久久久免费观看| 色综合欧美亚洲国产小说| 另类精品久久| 一级毛片黄色毛片免费观看视频| 在线观看国产h片| 国产日韩欧美亚洲二区| a级毛片在线看网站| 久久韩国三级中文字幕| 国产xxxxx性猛交| 波多野结衣一区麻豆| 又黄又粗又硬又大视频| 欧美日本中文国产一区发布| 一级毛片 在线播放| 建设人人有责人人尽责人人享有的| 看免费成人av毛片| 美女扒开内裤让男人捅视频| 在线观看www视频免费| 人体艺术视频欧美日本| 欧美日韩视频高清一区二区三区二| 在线观看人妻少妇| 街头女战士在线观看网站| av不卡在线播放| 高清视频免费观看一区二区| 成年女人毛片免费观看观看9 | 国产麻豆69| 亚洲国产精品成人久久小说| 99国产精品免费福利视频| 亚洲精品乱久久久久久| 国产成人欧美| 高清黄色对白视频在线免费看| 亚洲 欧美一区二区三区| 天天影视国产精品| 操出白浆在线播放| 老司机影院成人| 亚洲av日韩精品久久久久久密 | 午夜福利网站1000一区二区三区| 亚洲,一卡二卡三卡| 一区福利在线观看| 国产成人91sexporn| 黄片播放在线免费| 啦啦啦视频在线资源免费观看| 日韩视频在线欧美| 亚洲欧美成人综合另类久久久| 狂野欧美激情性bbbbbb| 成年美女黄网站色视频大全免费| 亚洲国产最新在线播放| 欧美黄色片欧美黄色片| 日韩av免费高清视频| 另类精品久久| 日韩av不卡免费在线播放| 99九九在线精品视频| 日本黄色日本黄色录像| 高清av免费在线| 97在线人人人人妻| 午夜av观看不卡| 在线 av 中文字幕| 夫妻性生交免费视频一级片| 搡老乐熟女国产| 少妇精品久久久久久久| 成人黄色视频免费在线看| 老司机在亚洲福利影院| 免费高清在线观看日韩| 国产成人av激情在线播放| 国产又爽黄色视频| 婷婷色综合大香蕉| 欧美精品亚洲一区二区| 最近中文字幕2019免费版| 黑人欧美特级aaaaaa片| 亚洲精品久久成人aⅴ小说| 精品国产一区二区久久| 老汉色∧v一级毛片| 亚洲激情五月婷婷啪啪| 超碰成人久久| 国产男女超爽视频在线观看| 国产日韩一区二区三区精品不卡| 日韩一区二区视频免费看| 性少妇av在线| 欧美成人精品欧美一级黄| 十八禁人妻一区二区| 亚洲欧美中文字幕日韩二区| 黑丝袜美女国产一区| 久久99热这里只频精品6学生| kizo精华| 在线观看三级黄色| 高清不卡的av网站| 啦啦啦 在线观看视频| 高清不卡的av网站| 高清视频免费观看一区二区| 丰满迷人的少妇在线观看| 欧美黑人精品巨大| 爱豆传媒免费全集在线观看| 人人妻,人人澡人人爽秒播 | 青春草视频在线免费观看| 亚洲美女视频黄频| 999精品在线视频| 国产 一区精品| 日韩 亚洲 欧美在线| 亚洲一码二码三码区别大吗| 丰满少妇做爰视频| 黄色一级大片看看| 丰满乱子伦码专区| 亚洲国产看品久久| 日日撸夜夜添| 狠狠精品人妻久久久久久综合| 91aial.com中文字幕在线观看| 亚洲欧美中文字幕日韩二区| 欧美xxⅹ黑人| 99九九在线精品视频| 国产成人精品福利久久| 久久综合国产亚洲精品| kizo精华| 久久久久精品国产欧美久久久 | 国产成人精品久久二区二区91 | 久久久久久人人人人人| 亚洲国产欧美一区二区综合| 在现免费观看毛片| 亚洲av电影在线进入| 啦啦啦中文免费视频观看日本| 两个人免费观看高清视频| 爱豆传媒免费全集在线观看| 熟妇人妻不卡中文字幕| 国产精品偷伦视频观看了| 黄色视频在线播放观看不卡| 一级毛片黄色毛片免费观看视频| 亚洲,一卡二卡三卡| 久久亚洲国产成人精品v| 一边亲一边摸免费视频| 亚洲成人免费av在线播放| 欧美日韩一级在线毛片| 丝袜在线中文字幕| 亚洲一区二区三区欧美精品| 9热在线视频观看99| 狠狠婷婷综合久久久久久88av| 午夜日韩欧美国产| av有码第一页| 色视频在线一区二区三区| 亚洲成人国产一区在线观看 | 一本—道久久a久久精品蜜桃钙片| 18禁动态无遮挡网站| 汤姆久久久久久久影院中文字幕| 少妇 在线观看| 男女下面插进去视频免费观看| 不卡av一区二区三区| 亚洲欧美精品自产自拍| 国产精品久久久久成人av| 天堂中文最新版在线下载| 欧美 亚洲 国产 日韩一| 欧美激情高清一区二区三区 | 亚洲国产日韩一区二区| 一区二区三区四区激情视频| 欧美激情极品国产一区二区三区| 日本av手机在线免费观看| 精品少妇黑人巨大在线播放| 啦啦啦视频在线资源免费观看| 国产成人精品福利久久| 国产av国产精品国产| 国产欧美亚洲国产| 亚洲色图综合在线观看| 国产成人一区二区在线| 女性被躁到高潮视频| 在线天堂最新版资源| 午夜福利视频精品| 精品人妻熟女毛片av久久网站| 国产一区二区 视频在线| 宅男免费午夜| 成人影院久久| 久久天堂一区二区三区四区| 亚洲精品国产色婷婷电影| 亚洲视频免费观看视频| 丰满乱子伦码专区| 久久久精品免费免费高清| 男女国产视频网站| 一区福利在线观看| 国产男女内射视频| 一区福利在线观看| 亚洲精品国产一区二区精华液| 一区福利在线观看| 一级片'在线观看视频| 国产又爽黄色视频| 2018国产大陆天天弄谢| 国产xxxxx性猛交| 国产精品三级大全| 人人妻人人爽人人添夜夜欢视频| 国产成人av激情在线播放| 高清视频免费观看一区二区| 日韩成人av中文字幕在线观看| 国产一区二区三区综合在线观看| av国产精品久久久久影院| 久久精品人人爽人人爽视色| 99久久综合免费| 最近中文字幕2019免费版| 亚洲综合精品二区| 亚洲一码二码三码区别大吗| 亚洲四区av| 人成视频在线观看免费观看| 最近2019中文字幕mv第一页| 亚洲av国产av综合av卡| 激情视频va一区二区三区| 国产精品成人在线| 亚洲国产成人一精品久久久| 观看av在线不卡| 国产乱人偷精品视频| 女人精品久久久久毛片|