日本地震對(duì)世界原油價(jià)格影響的干預(yù)分析 重慶工商大學(xué) 摘要 :近年來(lái)世界經(jīng)濟(jì)形勢(shì)不穩(wěn)定 ,原油價(jià)格的變動(dòng)很大程度上影響著各國(guó)經(jīng)濟(jì) 制定宏觀政策有極大的作用 型中利用時(shí)間序列的過(guò)去值預(yù)測(cè)未來(lái)值在統(tǒng)計(jì)預(yù)測(cè)分析中被廣泛運(yùn)用 紐約商品交易所交易的低硫輕質(zhì)原油的一種 ,是美國(guó)石油工業(yè)的基準(zhǔn)原油 日本發(fā)生 地震 ,在此期間 ,布倫特價(jià)格起伏 ,本文就是利用 型 就 日本地震對(duì)世界原油價(jià)格的影響進(jìn)行干預(yù)分析 震日（ 2011 年 3 月 11 日）前 60 日的 格和布倫特價(jià)格賦權(quán)平均后作為 世界原油價(jià)格 建立 型 ,預(yù)測(cè)出地震后 60 日的 世界原油價(jià)格 ,再與實(shí)際后 60 日的 格和布倫特價(jià)格賦權(quán)平均后的價(jià)格進(jìn)行比較 ,得出數(shù)據(jù) ,再進(jìn)行深入分析日本地震對(duì)原油價(jià)格的干預(yù) . 關(guān)鍵詞： 型 干預(yù)分析 模型 日本地震 世界原油價(jià)格 一、問(wèn)題的提出 2011 年 3 月 11 日日本發(fā)生 地震 ,這場(chǎng)地震引起全世界的關(guān)注 ,各國(guó)也針對(duì)日本地震對(duì)經(jīng)濟(jì)形勢(shì)的影響作出相應(yīng)的預(yù)測(cè)分析 世界原油價(jià)格作為現(xiàn)今測(cè)量經(jīng)濟(jì)形勢(shì)一大重要指標(biāo)并且與日本地震的關(guān)聯(lián)度最強(qiáng) ,被經(jīng)濟(jì)學(xué)家和普通百姓爭(zhēng)相關(guān)注 中國(guó)石化協(xié)會(huì)專門召開(kāi)了一次針對(duì)日本地震可能給石化行業(yè)造成何種影響的評(píng)估會(huì) ,并邀請(qǐng)了一些日本石化企業(yè)參加 ,給各方關(guān)注者解決了許多疑問(wèn) 大原油進(jìn)口國(guó) ,此次約有三分之一的產(chǎn)能受到影響 ,短期的需求銳減將打壓國(guó)際油價(jià) ,而地震導(dǎo)致的 核電站關(guān)閉估計(jì)會(huì)推高其他能源在震后重建中的需求 ,國(guó)際油價(jià)可能先降后升 日本地震的消息會(huì)成為國(guó)際油價(jià)最有力的主導(dǎo)因素：很多投資者 擔(dān)心嚴(yán)重的地震災(zāi)害令日本經(jīng) 濟(jì)的活力顯著喪失 ,對(duì)原油的需求也相應(yīng)下降 ,其造成的需求縮水幅度 ,足以對(duì)沖地緣政治局勢(shì)動(dòng)蕩造成的供應(yīng)短缺風(fēng)險(xiǎn) 界原油價(jià)格 帶來(lái)何種影響 ,以及其帶來(lái)的影響程度 ,我們從統(tǒng)計(jì)學(xué)的角度 就日本地震對(duì)原油價(jià)格的干預(yù)進(jìn)行深 1 層次研究 ,依靠準(zhǔn)確的數(shù)據(jù) ,進(jìn)行干預(yù)分析 ,得到更科學(xué)的分析結(jié)果 . 二、研究現(xiàn)狀及存在的問(wèn)題 近段時(shí)間以來(lái) ,國(guó)內(nèi)外學(xué)者和相關(guān)部門采用了多種方法對(duì)日本地震帶來(lái)的經(jīng)濟(jì)影響進(jìn)行了大量的探索 ,其中也不乏對(duì) 世界原油價(jià)格 影響的研究 于選取的指標(biāo)、研究的時(shí)點(diǎn)以及采取的方法不同 ,得出的結(jié)論也不 盡相同 . 從國(guó)內(nèi)學(xué)者的研究?jī)?nèi)容可以看出 ,對(duì)日本地震給石油價(jià)格帶來(lái)的影響的探索集中在對(duì)某幾個(gè)因素的研究上、缺乏綜合考慮 ,或以定性分析為主、缺少統(tǒng)計(jì)模型的有力支撐 ,或僅對(duì)影響的成因進(jìn)行分析、沒(méi)有對(duì)這種影響的未來(lái)走向進(jìn)行定量判斷 .（ 馬登科 , 張昕 , 世界原油價(jià)格決定的四個(gè)新特點(diǎn) ） 三、日本地震對(duì)國(guó)際油價(jià)的影響 日本地震引發(fā)海嘯 ,造成巨大損失 ,對(duì)石油需求的長(zhǎng)期和短期影響有不同的表現(xiàn) . 短期內(nèi) ,由于日本部分煉油廠和其他石油行業(yè)相關(guān)的基礎(chǔ)設(shè)施被地震和海嘯破壞 ,許多工廠停產(chǎn) ,短期經(jīng)濟(jì)的暫時(shí)停滯可能將使日本石油需求放緩 ,對(duì)于日本石油需求增長(zhǎng)的預(yù)期下降 ,往往會(huì)有 “應(yīng)聲而降 ”的作用 ,使得國(guó)際原油價(jià)格快速下挫 地震災(zāi)害對(duì)石油的需求也不一定減少 東京電力公司位于福島縣的核電在受到重創(chuàng) ,同時(shí)日本境內(nèi)所有核電站全部關(guān)閉 ,日本將失去近 1/3 的電力供應(yīng) ,火力發(fā)電和水力發(fā)電的功能也部分受損 電力供應(yīng)難以滿足需求 日本可能會(huì)增加火力發(fā)電 ,而柴油發(fā)電機(jī)由于其靈活簡(jiǎn)便 ,往往在自然災(zāi)害中起重要作用 日本對(duì)柴油的需求可能增加 ,對(duì)煤炭和油氣發(fā)電的需求也將增 加 . 中長(zhǎng)期 ,日本災(zāi)后重建過(guò)程需要增加 ,則可能對(duì)國(guó)際油價(jià)產(chǎn)生上升壓力 導(dǎo)致能源需求增加 核事故引發(fā)的核安全問(wèn)題再次引起國(guó)際社會(huì)對(duì)核電站安全的擔(dān)憂 ,這將可能影響全球核能的應(yīng)用 ,至少減緩核電發(fā)展進(jìn)程 ,一些核項(xiàng)目可能被推遲或取消 ,各國(guó)對(duì)核電項(xiàng)目將更加審慎 中長(zhǎng)期各國(guó)對(duì)非核能源 ,特別是對(duì)油氣的依賴可能會(huì)加強(qiáng) ,從而也或?qū)е聡?guó)際油氣價(jià)格的上浮 .（短期和長(zhǎng)期表現(xiàn)不同 日本地震影響油價(jià)）（鐘晶晶 ,日本強(qiáng)震致全球石油價(jià)格下挫） 四、模型建立前得準(zhǔn)備 （一）數(shù)據(jù)來(lái)源的說(shuō)明 2 布 倫特即 出產(chǎn)于北大西洋北海布倫特地區(qū)的輕質(zhì)低硫原油 ,是基準(zhǔn)品質(zhì) ,而且被廣泛交易 敦 洲際交易所和 美國(guó) 商品交易所均有他的期貨交易 ,是市場(chǎng)油價(jià)的標(biāo)桿 . 美國(guó) 西德 克 薩斯 輕質(zhì)原油 ,是北美地區(qū)較為通用的一類原油 目前 油已經(jīng)成為全球 原油定價(jià)的基準(zhǔn) 就短期來(lái)看 ,美元一樣 ,支配地位不可撼動(dòng) . 油多數(shù)時(shí)間內(nèi)跟隨 油走勢(shì) ,少數(shù)時(shí)間出現(xiàn)背離或影響 價(jià) 還有 油 ,亞洲 非洲 油 ,都是地區(qū)性基準(zhǔn)原油 ,但因影響力較小 ,且多數(shù)時(shí)間跟隨 油走勢(shì) ,故而較少作為參考 . 2011年 3月 11 日日本大地震對(duì)世界經(jīng)濟(jì)的沖擊 ,大宗商品 價(jià)格出現(xiàn)調(diào)整 ,國(guó)際原油受影響強(qiáng)勁反彈 以 2011 年 3 月 11 日為分界線 ,前后分別選取庫(kù)倫特和 60 個(gè)數(shù)據(jù)（國(guó)際石油網(wǎng)） ,然而 ,世界原油三分之二以上的交易量 ,是以同樣輕質(zhì)低硫的北海布倫特原油為基準(zhǔn)油作價(jià) 行賦權(quán)以求得的平均值作為 世界原油價(jià)格 ,即 = 3 + 4 T 原 油 價(jià) 格 布 倫 特 期 貨 價(jià) 格 期 貨 價(jià) 格（二）模型的介紹 1、 型 型全稱為自回歸移動(dòng)平均模型 (記 是由 博克思 ( 詹金斯 ( 70 年代初提出的一著名 時(shí)間序列預(yù)測(cè)方法 ,所以又稱為 型、博克思 其中 p,d,q）稱為差分自回歸移動(dòng)平均模型 ,自回歸 , p 為自回歸項(xiàng) ; 移動(dòng)平均 ,q 為移動(dòng)平均項(xiàng)數(shù) ,d 為時(shí)間序列成為平穩(wěn)時(shí)所做的差分次數(shù) . 型的基本思想是：將預(yù)測(cè)對(duì)象隨時(shí)間推移而形成的數(shù)據(jù)序列視為一個(gè)隨機(jī)序列 ,用一定的 數(shù)學(xué)模型 來(lái)近似描述這個(gè)序列 現(xiàn)代統(tǒng)計(jì)方法、 計(jì)量經(jīng)濟(jì)模型 在某種程度上已經(jīng)能夠幫助企業(yè)對(duì)未來(lái)進(jìn)行預(yù)測(cè) . 型預(yù)測(cè)的基本程序是： 1、根據(jù)時(shí)間序列 圖 、自相關(guān)函數(shù)和偏自相關(guān)函數(shù)圖以位根檢驗(yàn)其 方差 、趨勢(shì)及其季節(jié)性變化規(guī)律 ,對(duì)序列的平穩(wěn)性進(jìn)行識(shí)別 經(jīng)濟(jì)運(yùn)行的時(shí)間序列都不是平穩(wěn)序列 . 2、對(duì)非平穩(wěn)序列進(jìn)行平穩(wěn)化處理 并存在一定的增長(zhǎng)或下降趨勢(shì) ,則需要對(duì)數(shù)據(jù)進(jìn)行差分處理 ,如果數(shù)據(jù)存在異方差 ,則需對(duì)數(shù)據(jù) 3 進(jìn)行技術(shù)處理 ,直到處理后的數(shù)據(jù)的自 相關(guān)函數(shù)值和偏相關(guān)函數(shù)值無(wú)顯著地異于零 . 3、根據(jù)時(shí)間序列模型的識(shí)別規(guī)則 ,建立相應(yīng)的模型 而自相關(guān)函數(shù)是拖尾的 ,可斷定序列適合 型；若平穩(wěn)序列的偏相關(guān)函數(shù)是拖尾的 ,而自相關(guān)函數(shù)是截尾的 ,則可斷定序列適合 型；若平穩(wěn)序列的偏相關(guān)函數(shù)和自相關(guān)函數(shù)均是拖尾的 ,則序列適合 型 . 4、進(jìn)行 參數(shù)估計(jì) ,檢驗(yàn)是否具有統(tǒng)計(jì)意義 . 5、進(jìn)行 假設(shè)檢驗(yàn) ,診斷殘差序列是否為白噪聲 . 6、利用已通過(guò)檢驗(yàn)的模型進(jìn)行 預(yù)測(cè)分析 . 圖 1 2、干預(yù)分析模型 時(shí)間序列經(jīng)常會(huì)受到特殊事件及態(tài)勢(shì)的影響 ,諸如國(guó)內(nèi) ,國(guó)際經(jīng)濟(jì)政策或經(jīng)濟(jì)規(guī)則的變更 ,以及罷工、促銷之類事件的影響等等 ,學(xué)者們稱這類外部事件為干預(yù) 當(dāng)不存在影響時(shí) ,這兩種模型沒(méi)有本質(zhì)上的區(qū)別 . 研究干預(yù)分析模型的目的：從定量分析的角度來(lái)評(píng)估政策干預(yù)或突發(fā)事件對(duì)經(jīng)濟(jì)環(huán)境和經(jīng)濟(jì)過(guò)程的具體影響 干預(yù)分析模型是與時(shí)間序列模型結(jié)合起來(lái)進(jìn)行研究的 在干預(yù)事件發(fā)生后序列是否存在任 何事實(shí)上的變化？若有影響 ,其影響程度又如何？這就是干預(yù)分析模型所要解決的問(wèn)題 . 干預(yù)變量的形式 ：干預(yù)分析模型的基本變量是干預(yù)變量 ,有兩種常見(jiàn)的干預(yù)變量 . 一種是持續(xù)性的干預(yù)變量 ,表示 T 時(shí)刻發(fā)生以后 , 一直有影響 ,這時(shí)可以用階躍函數(shù)表示 , 4 形式是： 0,1, 干 預(yù) 事 件 發(fā) 生 之 前 （ ）干 預(yù) 事 件 發(fā) 生 之 后 （ ）第二種是短暫性的干預(yù)變量 ,表示在某時(shí)刻發(fā)生 , 僅對(duì)該時(shí)刻有影響 , 用單位脈沖函數(shù)表示 ,形式是： 1,0, 干 預(yù) 事 件 發(fā) 生 時(shí) （ ）其 它 時(shí) 間 （ ）干預(yù)事件的形式 ：干預(yù)事件雖然多種多樣 ,但按其影響的形式 ,歸納起來(lái)基本上有四種類型： a. 干預(yù)事件的影響突然開(kāi)始 ,長(zhǎng)期持續(xù)下去 從某一時(shí)刻 T 開(kāi)始 ,但影響的程度是未知的 ,即因變量的大小是未知的 表示干預(yù)影響強(qiáng)度的未知參數(shù) 平穩(wěn)時(shí)可以通過(guò)差分化為平穩(wěn)序列 ,則干預(yù)模型可調(diào)整為： (1 ) S,其中 B 為后移算子 如 b 個(gè)時(shí)期 ,那么干預(yù)模型可進(jìn)一步調(diào)整為 ： S b. 干預(yù)事件的影響逐漸開(kāi)始 ,長(zhǎng)期持續(xù)下去 并不能立刻產(chǎn)生完全的影響 ,而是隨著時(shí)間的推移 ,逐漸地感到這種影響的存在 , 0 11 更一般的模型是 ： 1 , 0 11 . . tt c. 干預(yù)事件突然開(kāi)始 ,產(chǎn)生暫時(shí)的影響 0 11 b 時(shí) ,干預(yù)的影響只存在一個(gè)時(shí)期 ,當(dāng) 1 時(shí) ,干預(yù)的影響將長(zhǎng)期存在 . d. 干預(yù)事件逐漸開(kāi)始 ,產(chǎn)生暫時(shí)的影響 在某個(gè)時(shí)刻到達(dá)高峰 ,然后又逐漸減弱以至消失 011 . . . 綜合上述 ,不管經(jīng)濟(jì)系統(tǒng)受到的干預(yù)影響多么復(fù)雜 ,都可以用上述四種形式或者是它們的某種組合來(lái)表達(dá) 也可以用這種組合去擬合多個(gè)干預(yù)事件所產(chǎn)生的影響 . 干預(yù)分析模型建模的思路和具體建 模步驟包括以下幾點(diǎn)： 1、利用干預(yù)影響產(chǎn)生前 的 數(shù) 5 據(jù)建立一個(gè)單變量的時(shí)間序列模型 ,然后利用此模型進(jìn)行外推 ,得到的預(yù)測(cè)值作為不受干預(yù)影響的數(shù)值 實(shí)際值減去預(yù)測(cè)值 ,得到受干預(yù)影響的具體結(jié)果 ,利用這些結(jié)果求估干預(yù)影響的參數(shù) 用排除干預(yù)影響后的全部數(shù)據(jù) ,識(shí)別與估計(jì)出一個(gè)單變量的時(shí)間序列模型 出總的干預(yù)分析模型 .（徐國(guó)祥 ,統(tǒng)計(jì)預(yù)測(cè)與決策 干預(yù)分析模型預(yù)測(cè)法） 五、 世界原油價(jià)格 型的構(gòu)建、檢驗(yàn)及預(yù)測(cè) （一）序列的檢驗(yàn)及處理 對(duì)于建立 型要求時(shí)間序列數(shù)據(jù) 是 平穩(wěn)的（所謂平穩(wěn)時(shí)間序列是 指隨機(jī)時(shí)間序列 的統(tǒng)計(jì)特征不隨時(shí)間的 推移而變化） 觀察其走勢(shì)和 然該時(shí)間序列不是平穩(wěn)序列 ,而是隨著時(shí)間的增加而不斷增大的 ,而且波動(dòng)幅度較大的序列 . 世界石油價(jià)格時(shí)間序列圖8590951001051101151 4 7 10 13 16 19 22 25 28 31 34 37 40 43 46 49 52 55 58世界石油價(jià)格圖 2 原始序列相關(guān)圖 圖 3 為了得到平穩(wěn) 時(shí)間序列 ,對(duì)該時(shí)間序列通過(guò)做差分進(jìn)行平穩(wěn)化處理 ,經(jīng)平穩(wěn)化處理后的二階差分的時(shí)間序列相關(guān)圖 ,可以看出新的 時(shí)間 序列已經(jīng)平穩(wěn)化了 ,為了保險(xiǎn)還需要對(duì)二階差分序列進(jìn)行 位根檢驗(yàn) 變量二階差分序列的 位根檢驗(yàn)對(duì)應(yīng)的統(tǒng)計(jì)量 的 P 值均小于 5%,表示在 5%的顯著水平下 ,二階差分序列存在單位根的原 假設(shè)被拒絕 ,即二階差分序列是平穩(wěn)的 ,可以對(duì)其建立 型 . 6 圖 4 二階差分序列相關(guān)圖 圖 5 % % 0% 位根檢驗(yàn) 表 1 世界原油價(jià)格二階差分后的序列圖 4 2 4 6 1 3 5 7 9 11 13 15 17 19 21 23 25 27 29 31 33 35 37 39 41 43 45 47 49 51 53 55 57 二階差分序列 7 （二）模型的識(shí)別與定階 根據(jù)圖 5,我們可以看到其自相關(guān)系數(shù)拖尾 ,偏自相關(guān)系數(shù) 2 階截尾 ,由此可以初步判斷p,d,q）模型中的 2p , 2d , 0q ,即使用 2,2,0） 通過(guò)直接觀察來(lái)確定模型的階數(shù)顯然是不夠準(zhǔn)確的 我們可以利用階數(shù)外擴(kuò)法及最小則來(lái)識(shí)別模型的類型和階數(shù) 件的反復(fù)計(jì)算比較 ,獲得 最后 的模型形式為 4,2,0） . 模型比較 表 2 （三）模型的參數(shù)估計(jì)和假設(shè)檢驗(yàn) 從表 3 可以看出 3）是不顯著的 ,因此將其從模型中剔除 ,根據(jù) 出數(shù)據(jù)建立 4,2,0）模型 ,即： 2 4 21 0 . 6 6 3 0 1 1 0 . 4 5 4 1 7 6 0 . 2 3 1 9 6 5 ( 1 ) B B x 從中可以看出參數(shù)的估計(jì)值都是符合統(tǒng)計(jì)意義的 從殘差序列相關(guān)圖可以看出自相關(guān)系數(shù)趨于 0,落入一個(gè)隨機(jī)區(qū)間 ,與 0 無(wú)顯著性差異 ,而起 P 值多大于 差為白噪聲序列 ,因此診斷該模型是可靠的 ,可預(yù)測(cè)的 . C R(1) ) R(3) R(4) 數(shù)估計(jì) 表 3 P 2 3 4 Q 0 0 0 8 殘差序列相關(guān)圖 圖 6 （四）模型的預(yù)測(cè) 根據(jù)式 2 4 21 0 . 6 6 3 0 1 1 0 . 4 5 4 1 7 6 0 . 2 3 1 9 6 5 ( 1 ) B B x 得到日本地震后 的 60 個(gè)世界原油價(jià)格 的預(yù)測(cè)值 ,世界原油價(jià)格 擬合及預(yù)測(cè)效果圖見(jiàn)圖 7 世界原油價(jià)格 擬合及預(yù)測(cè)效果圖 圖 7 9 六、日本地 震干預(yù)分析模型的建立和預(yù)測(cè) （一）分離出干預(yù)影響的具體數(shù)據(jù) ,求估干預(yù)模型的參數(shù) . 運(yùn)用模型估計(jì)出來(lái)的數(shù)據(jù)得到日本地震的效益值 t t tZ x x 由于日本地震的發(fā)生并不是立刻產(chǎn)生完全的影響 ,而是隨著時(shí)間的推移 ,逐漸的感到這種影響的存在 ,因此干預(yù)模型選取以下形式： 1 ,其中 , 0 , 2 0 1 1 3 1 11 , 2 0 1 1 3 1 1 年 月 日 以 前年 月 日 及 以 后對(duì)模型進(jìn)行參數(shù)估計(jì) ,得到的結(jié)果如下： t t - 1 = 0 . 7 7 6 0 8 9 , = 0 . 8 7 6 8 9 20 . 8 7 6 8 9 2 0 . 7 7 6 0 8 9其中 2 0 . 7 8 9 4 9 5 , F = 2 1 3 . 7 7 6 8 ( 0 . 0 0 0 0 0 0 顯 著 ）,模型系數(shù)的 t 檢驗(yàn)也是顯著的 ,說(shuō)明模型的擬合效果好 . （二）計(jì)算凈化序列 ,建立相應(yīng)的時(shí)間序列模型 凈化序列是指消除了日本地震影響的時(shí)間序列 ,是由原始序列 去干預(yù)影響值 到 ,即： 0 . 7 7 6 0 8 9y , 6 1 , 1 , 2 ? 1201 0 . 8 7 6 8 9 2 Tt t T ， ， 利用凈化序列 ,再次建立 型進(jìn)行擬合 ,經(jīng)過(guò)二階差 分后 ,凈化序列可以建立2,2,1） 通過(guò)檢驗(yàn)得到殘差序列為白噪聲序列 ,得到的模型如下： 221 0 . 3 7 0 7 1 9 B 0 . 3 6 9 8 6 2 ( 1 ) ( 1 0 . 9 8 8 2 7 6 B ) y 其中 2 0 . 6 8 0 2 3 9 , F = 7 9 . 4 2 0 7 0 ( 0 . 0 0 0 0 0 0 顯 著 ）,模型系數(shù)的 t 檢驗(yàn)也是顯著的 ,說(shuō)明模型的擬合效果好 . （三）組建干預(yù)分析模型并進(jìn)行預(yù)測(cè) 結(jié)合 擬合模型與干預(yù)參數(shù) ,的估計(jì)值 ,得到所求的干預(yù)分析預(yù)測(cè)模型：1 2 3 4 1 0 . 7 7 6 0 8 91 . 6 2 9 2 8 1 0 . 6 2 8 4 2 4 0 . 3 6 9 0 0 5 0 . 3 6 9 8 6 2 0 . 9 8 8 2 7 61 0 . 8 7 6 8 9 2 Tt t t t t t t tx y y y y 其中0 , 2 0 1 1 3 1 11 , 2 0 1 1 3 1 1年 月 日 以 前 (t 石油天然氣 綜合 分析綜述 附錄 附錄 1： 布倫特和 賦權(quán)平均后的 數(shù)據(jù) 序列號(hào) 日期 布倫特 界原油價(jià) 格 序列號(hào) 日期 布倫特 界原油價(jià)格 1 20101213 1 20110311 20101214 2 20110314 20101215 3 20110315 20101216 4 20110316 20101217 5 20110317 20101220 6 20110318 20101221 7 20110321 12 8 20101222 8 20110322 04 20101223 9 20110323 0 20101227 1 0 20110324 1 20101228 1 20110325 2 20101229 2 20110328 3 20101230 3 20110329 4 20101231 4 20110330 5 20110103 5 20110331 6 20110104 6 20110401 7 20110105 7 20110404 8 20110106 8 20110405 9 20110107 9 20110406 0 20110110 0 20110407 1 20110111 1 20110408 2 20110112 2 20110411 3 20110113 3 20110412 4 20110114 4 20110413 5 20110118 5 20110414 6 20110119 6 20110415 7 20110120 7 20110418 8 20110121 8 20110419 9 20110124 9 20110420 0 20110125 0 20110421 1 20110126 1 20110422 2 20110127 2 20110425 3 20110128 3 20110426 4 20110131 4 20110427 5 20110201 5 20110428 6 20110202 6 20110429 7 20110203 7 20110502 8 20110204 8 20110503 9 20110207 9 20110504 0 20110208 00 20110505 1 20110209 01 20110506 2 20110210 02 20110509 3 20110211 03 20110510 4 20110214 04 20110511 5 20110215 05 20110513 6 20110216 06 20110513 7 20110218 07 20110516 8 20110222 08 20110517 9 20110223 09 20110518 13 50 20110224 10 20110519 1 20110225 11 20110520 2 20110228 12 20110523 07 53 20110301 13 20110524 4 20110302 14 20110525 5 20110303 15 20110526 6 20110304 16 20110527 7 20110307 17