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

    基于經(jīng)驗模態(tài)分解的煤炭消費量組合預測

    2016-11-03 05:11:52鞠金艷祝榮欣
    黑龍江科技大學學報 2016年1期
    關(guān)鍵詞:消費量波動煤炭

    鞠金艷, 祝榮欣, 陳 玉

    (黑龍江科技大學 機械工程學院, 哈爾濱 150022)

    ?

    基于經(jīng)驗模態(tài)分解的煤炭消費量組合預測

    鞠金艷,祝榮欣,陳玉

    (黑龍江科技大學 機械工程學院, 哈爾濱 150022)

    為了明確煤炭消費量的發(fā)展變化規(guī)律和趨勢,科學引導煤炭行業(yè)健康、有序發(fā)展。針對煤炭消費量的時間序列變化具有增長性、波動性和非平穩(wěn)性的特點,采用經(jīng)驗模態(tài)分解法對1978~2014年煤炭消費量進行多層次分解,得到其發(fā)展變化的趨勢量和波動量。利用BP神經(jīng)網(wǎng)絡的非線性映射能力,分別對趨勢量和波動量進行預測,最終二者相加求和得到煤炭消費量的預測值。誤差分析表明,基于經(jīng)驗模態(tài)分解的煤炭消費量組合預測模型,擬合值的平均誤差為2.18%,預測值的平均誤差為1.24% 。該組合預測模型可以有效的提高煤炭消費量的預測精度,用該模型預測了2015~2020年煤炭消費量。預測結(jié)果表明,在未來幾年煤炭消費量將保持低速增長趨勢,到2020年將達到341 718.2萬t標準煤。

    煤炭消費量; 經(jīng)驗模態(tài)分解; 組合預測; BP網(wǎng)絡; 時間序列

    0 引 言

    煤炭是我國重要的基礎(chǔ)能源,經(jīng)濟發(fā)展、工業(yè)化和城鎮(zhèn)化進程都離不開煤炭資源的支持。在未來一段時期,煤炭仍然是能源消費的主體,預測煤炭消費量的動態(tài)發(fā)展變化趨勢,對科學測算煤炭開采規(guī)模,合理開發(fā)和利用煤炭資源,努力保持煤炭供需基本平衡,相關(guān)部門和機構(gòu)更合理地制定煤炭工業(yè)發(fā)展政策、措施,加強對煤炭總量的宏觀調(diào)控,確保煤炭行業(yè)平穩(wěn)健康發(fā)展等具有重要的意義[1-2]。國外對煤炭消費的研究主要集中于研究煤炭消費量的波動情況、價格的變化及其主要影響因素等,研究的方法主要有向量自回歸模型、小波分析、ADL空間分布模型等[3-5]。我國學者預測煤炭消費量相關(guān)的研究較多,主要分為三類:一類是對煤炭消費量時間序列進行預測,方法主要有指數(shù)平滑、彈性系數(shù)法、馬爾科夫預測模型、灰色GM(1,1)模型、支持向量機、神經(jīng)網(wǎng)絡預測模型、組合預測等[6-7];第二類是通過分析影響煤炭消費的主要因素,然后采用灰色GM(1,n)模型,多元回歸模型、自回歸模型或神經(jīng)網(wǎng)絡模型等,建立煤炭消費量和主要影響因素之間的關(guān)系模型,進而預測出煤炭消費量的變化情況[8-9];第三類是通過分析主要消耗煤炭資源的各部門的發(fā)展變化情況,得出煤炭消費量變化情況[10]。這些方法各有其優(yōu)點和不足,時間序列預測方法能很好的預測未來發(fā)展趨勢,但難以準確預測出煤炭消費量系統(tǒng)中的非線性關(guān)系,導致預測精度不高;根據(jù)主要影響因素的變化,能預測未來的煤炭消費情況,但需要事先預測出各影響因素的未來發(fā)展變化情況,因此操作起來復雜,而且影響因素的選取直接影響分析結(jié)果的差異。

    煤炭消費量受諸多因素的影響,發(fā)展變化既具有隨時間推移的增長性趨勢,又具有一定的波動性和非平穩(wěn)性的特點。因此,直接對煤炭消費量時間序列進行預測很難得到準確的預測結(jié)果。為提高預測結(jié)果可靠性,筆者采用經(jīng)驗模態(tài)分解法分解出煤炭消費量時間序列的變化趨勢量和波動量,分別對變化趨勢量和波動量進行預測,最終將二者預測結(jié)果相加求和得到煤炭消費量的預測值。

    1 建模思想與方法

    1.1建模思想

    煤炭消費量時間序列的增長性、波動性和非平穩(wěn)性的變化特點,決定了直接采用時間序列預測方法很難準確預測出其未來發(fā)展情況。Huang N.E.等提出的經(jīng)驗模態(tài)分解法(Empirical mode decomposition,EMD),是一種處理非平穩(wěn)、非線性信號的時頻信號分析法,對時間序列信號的分解具有客觀性和穩(wěn)定性的優(yōu)勢[11-12]。經(jīng)驗模態(tài)分解法通過對信號序列進行逐級線性化和平穩(wěn)化處理,分解出原信號序列的發(fā)展趨勢量和所有包含原信號不同時間尺度特征信息的波動分量,即本征模態(tài)函數(shù)。文獻[12]中詳細介紹了經(jīng)驗模態(tài)分解法的基本原理和計算過程,文中不再詳細介紹。采用經(jīng)驗模態(tài)分解法對煤炭消費量時間序列進行分解,得到煤炭消費的趨勢量和波動量,再選擇合適的方法分別對趨勢量和波動量進行預測,可以提高煤炭消費量預測結(jié)果的準確度[13-14]。

    BP神經(jīng)網(wǎng)絡具有較好的逼近非線性函數(shù)的能力,在預測領(lǐng)域應用得非常廣泛[15-18]。采用BP神經(jīng)網(wǎng)絡對分解出的煤炭消費趨勢量和波動量分別進行預測。BP神經(jīng)網(wǎng)絡由輸入層、若干個隱含層和輸出層構(gòu)成,如圖1所示。圖1是含有兩個隱含層I和J的BP網(wǎng)絡結(jié)構(gòu)。BP神經(jīng)網(wǎng)絡具有逼近非線性函數(shù)的能力,但有時會陷入局部極小值,影響泛化能力,進而影響預測精度。研究表明,直接采用BP神經(jīng)網(wǎng)絡對時間序列進行預測,必須采用滾動預測的方法,這種方法要求樣本量大而且預測精度不是很高[15,17]。傳統(tǒng)的時間序列預測方法能很好的擬合數(shù)據(jù)序列的發(fā)展趨勢,但不適合逼近復雜的非線性函數(shù)。因此,可以將BP神經(jīng)網(wǎng)絡和傳統(tǒng)的時間序列預測方法相結(jié)合,取長補短,采用組合預測的方法對時間序列進行預測。

    圖1 兩個隱含層的BP網(wǎng)絡結(jié)構(gòu)

    1.2建模方法

    基于經(jīng)驗模態(tài)分解法建立煤炭消費量組合預測模型,具體建模步驟如下:

    (1)采用經(jīng)驗模態(tài)分解法對1978~2014年煤炭消費量進行分解,得到煤炭消費趨勢量和波動量,然后選擇預測方法分別對趨勢量和波動量進行預測。

    (2)為了準確的預測出煤炭消費趨勢量的變化,采用傳統(tǒng)預測方法對煤炭消費趨勢量建立多種單一預測模型,并對建立的模型進行顯著性和擬合優(yōu)度檢驗,計算擬合值的平均絕對百分誤差,為了保證預測精度,一般誤差值應小于10%。

    (3)用建立的各單一預測模型擬合1978~2014年煤炭消費趨勢量的值,將其中1978~2011年煤炭消費趨勢量的擬合值作為BP網(wǎng)絡的輸入向量(P),煤炭消費趨勢量的實際值作為BP網(wǎng)絡的輸出向量(T),運用MATLAB7.6.0語言編寫程序,確定網(wǎng)絡參數(shù),建立組合預測模型中各單一預測模型擬合值與實際值之間的非線性映射關(guān)系。

    (4)判斷BP網(wǎng)絡的學習精度是否達到要求,并以2012~2014年煤炭消費趨勢量的數(shù)據(jù)作為檢驗樣本,來驗證網(wǎng)絡的泛化能力是否滿足要求,若滿足要求,運算停止;否則,重新設定網(wǎng)絡結(jié)構(gòu)的相關(guān)參數(shù),繼續(xù)訓練學習。

    (5)將各單一預測模型預測出的2015~2020年煤炭消費趨勢量的數(shù)值輸入訓練好的BP神經(jīng)網(wǎng)絡,仿真后得到相應的煤炭消費趨勢量預測值。

    (6)預測煤炭消費波動量的發(fā)展變化,由于其波動的非平穩(wěn)性,傳統(tǒng)預測方法無法對其進行準確預測,因此,采用BP神經(jīng)網(wǎng)絡進行滾動預測,確定滾動預測模型的相關(guān)參數(shù),實現(xiàn)非線性擬合,檢驗擬合精度和泛化能力,進而對煤炭消費波動量進行預測。

    (7)將1978~2014年煤炭消費趨勢量和波動量的預測值相加求和,得到總的煤炭消費量預測值,并計算預測值和實際值的誤差,若滿足要求,則用其對2015~2020年煤炭消費量進行預測。

    2 煤炭消費量的預測

    2.1數(shù)據(jù)的來源

    文中主要使用1978~2014年煤炭消費量的數(shù)據(jù),是折算成標準量(標準煤)后的數(shù)據(jù),該數(shù)據(jù)不能直接查到,但可通過查閱能源消費總量的統(tǒng)計數(shù)據(jù),并將其與煤炭消費量占能源消費總量的比重進行乘積,計算得到煤炭消費量的數(shù)據(jù)。1978~1985年數(shù)據(jù)來源為《新中國六十年統(tǒng)計資料匯編》,1986~2014年數(shù)據(jù)來源為《中國統(tǒng)計年鑒》。

    2.2煤炭消費量的經(jīng)驗模態(tài)分解

    煤炭消費量的發(fā)展變化受經(jīng)濟增長、煤炭價格、產(chǎn)業(yè)結(jié)構(gòu)、能源結(jié)構(gòu)和人口數(shù)量等多種因素的共同影響,煤炭消費量時間序列既有明顯的增長趨勢,又具有一定的波動性和非平穩(wěn)性的特點。采用經(jīng)驗模態(tài)分解法,應用Matlab軟件編程,對1978~2014年煤炭消費量時間序列進行分解,得到1個趨勢量和1個波動量,分解結(jié)果如圖2和表1所示。由圖2可知,煤炭消費趨勢量反映了各種影響因素所引起的煤炭消費量的長期增長趨勢,波動量主要反映了煤炭消費量總體的波動情況,短周期的波動較普遍發(fā)生。煤炭消費趨勢量的增長變化是煤炭消費總量變化的主導因素,煤炭消費波動量變化對煤炭消費總量的影響相對較小,因此,準確的預測煤炭消費趨勢量的增長變化情況,直接影響煤炭消費總量預測結(jié)果的準確性。

    圖2 煤炭消費量及其經(jīng)驗模態(tài)分解結(jié)果

    Fig. 2Coal consumption and its empirical mode decomposition

    2.3煤炭消費趨勢量的預測

    2.3.1單一預測模型的建立

    利用SPSS軟件對1978~2011年煤炭消費趨勢量時間序列數(shù)據(jù)進行曲線擬合,得到擬合精度較高的預測模型,分別為指數(shù)函數(shù)模型和一元線性回歸模型,如式(1)和(2)所示。指數(shù)函數(shù)模型的判定系數(shù)R2=0.979,一元線性回歸模型的判定系數(shù)R2=0.865,兩個模型都通過 檢驗,模型極顯著。用式(1)和(2)對我國1978~2011年煤炭消費趨勢量值進行預測,得到預測值。

    Y=35 346.154e0.057x,

    (1)

    Y=2 042.373+6 267.658x,

    (2)

    式中:Y——煤炭消費趨勢量值,萬t;

    x——時間變量,1978~2011年的對應取值分別為1~34。

    2.3.2 BP神經(jīng)網(wǎng)絡組合預測模型的建立

    (1)BP網(wǎng)絡輸入和輸出向量的確定。BP網(wǎng)絡輸入向量(P)是指數(shù)函數(shù)模型和一元線性回歸模型2種預測方法預測出的1978~2011年煤炭消費趨勢量的預測值,因此BP網(wǎng)絡輸入層有2個輸入節(jié)點,BP網(wǎng)絡的輸出向量(T)是煤炭消費趨勢量的實際值,因此,輸出層有1個輸出節(jié)點,從而建立了2個預測模型的預測值與相應實際值之間的非線性映射關(guān)系。

    表1 煤炭消費量分解及預測結(jié)果與誤差

    (2)隱含層及節(jié)點數(shù)的確定。逼近任何一個非線性函數(shù)可用含有一個隱含層的3層BP神經(jīng)網(wǎng)絡實現(xiàn)[13-15],文中選用1個隱含層的網(wǎng)絡。隱含層節(jié)點數(shù)采用逐漸增加的方法確定為3最理想,即BP網(wǎng)絡結(jié)構(gòu)為2-3-1。

    (3)BP網(wǎng)絡隱含層轉(zhuǎn)換函數(shù)為“Sigmoid”,輸出層轉(zhuǎn)換函數(shù)為“Purelin”,訓練函數(shù)采用預測精度和速度上都優(yōu)于其他算法的“Trainlm”,網(wǎng)絡訓練的相關(guān)參數(shù)為學習精度5×10-5,學習速率0.01,迭代步數(shù)2 000,對網(wǎng)絡進行訓練學習。

    (4)判斷BP網(wǎng)絡的學習精度是否達到要求,并以2012~2014年煤炭消費趨勢量的數(shù)據(jù)作為檢驗樣本,來驗證網(wǎng)絡的泛化能力是否滿足要求。若滿足要求,運算停止,得到相應節(jié)點的權(quán)值和閾值;否則,重新設定網(wǎng)絡結(jié)構(gòu)的相關(guān)參數(shù),繼續(xù)訓練學習。文中構(gòu)建的BP網(wǎng)絡在訓練精度MSE為4.15×10-5時,精度達到要求,訓練好的BP網(wǎng)絡的煤炭消費趨勢量的擬合值、預測值及其誤差值e如表1。由表1可知,BP神經(jīng)網(wǎng)絡組合預測模型得到的1978~2011年煤炭消費趨勢量擬合值的平均誤差為2.7%, 2012~2014年煤炭消費趨勢量預測值的平均誤差為1.2%,可見,組合預測模型既有較好的擬合能力和預測能力,可以用此模型對2015~2020年煤炭消費趨勢量進行預測。

    (5)將指數(shù)函數(shù)模型和一元線性回歸模型預測出的2015~2020年煤炭消費趨勢量的值輸入訓練好的BP神經(jīng)網(wǎng)絡預測模型,仿真后得到2015~2020年煤炭消費趨勢量的預測值,預測結(jié)果見表2。

    表22015~2020年煤炭消費量預測值

    Table 2Predicted values of coal consumption from 2015 to 2020 萬t

    年份煤炭消費趨勢量煤炭消費波動量煤炭消費總量2015293185.01702.7294887.72016297667.04250.1301917.12017302867.05192.5308059.52018311073.03690.9314763.92019325120.0-1688.9323431.12020343844.0-2125.8341718.2

    2.4煤炭消費波動量的預測

    由圖2可知,1978~2014年煤炭消費波動量時間序列的波動性較強,采用曲線回歸、自回歸模型和灰色預測等方法很難對其進行準確預測,因此,利用BP神經(jīng)網(wǎng)絡具有逼近非線性函數(shù)的特點,對其進行預測。采用BP神經(jīng)網(wǎng)絡對時間序列進行預測必須采用滾動預測的方法[13],即用前幾年的數(shù)據(jù)預測下一年的數(shù)據(jù)。

    BP神經(jīng)網(wǎng)絡的煤炭消費波動量預測模型的建模過程與預測趨勢量的建模過程基本相同。根據(jù)圖2的煤炭消費波動量變化規(guī)律及文獻[19]中采用BP 濾波對 1953~2012 年煤炭消費需求波動曲線的變化特征分析可知,煤炭消費波動量時間序列的波動周期為3~6 a,因此將BP網(wǎng)絡輸入層節(jié)點數(shù)分別確定為3個、4個和5個進行對比試驗,最終確定BP網(wǎng)絡輸入層節(jié)點為4個最理想,輸入向量(P)是1978~2010年煤炭消費波動量的實際值,隨著時間序列逐年滾動,每4年為一組作為輸入向量, BP網(wǎng)絡的輸出向量(T)是煤炭消費波動量輸入向量的下一年的實際值,即應為1982~2011年煤炭消費波動量的實際值,因此輸出層有1個輸出節(jié)點。例如:將1978~1981年煤炭消費波動量的實際值輸入網(wǎng)絡,可預測出1982年煤炭消費波動量的值,逐年滾動,共有30組訓練樣本。BP網(wǎng)絡的隱含層數(shù)選擇1,采用逐漸增加節(jié)點的方法確定隱含層節(jié)點數(shù)為4,BP網(wǎng)絡結(jié)構(gòu)為4-4-1。BP網(wǎng)絡各層函數(shù)的選擇與趨勢量預測相同,網(wǎng)絡訓練的參數(shù)為學習精度10-6,學習速率0.01,迭代步數(shù)3 000,對網(wǎng)絡進行訓練。以2012~2014年煤炭消費波動量的數(shù)據(jù)作為檢驗樣本,驗證網(wǎng)絡的泛化能力是否滿足要求,用訓練好的BP網(wǎng)絡預測煤炭消費波動量,擬合值、預測值及其相對誤差值如表1。

    由表1可知,BP神經(jīng)網(wǎng)絡預測得到的1982~2011年煤炭消費波動量擬合值的平均誤差為6.6%,2012~2014年煤炭消費波動量檢驗樣本預測值的平均誤差為3.7%,可見BP網(wǎng)絡預測的煤炭消費波動量擬合效果一般,短期預測效果很好。由于煤炭消費波動量在整個煤炭消費總量中所占的比重較小,預測精度稍低不會對整個模型的預測效果產(chǎn)生大的影響。將2011~2014年煤炭消費波動量的值輸入訓練好的BP神經(jīng)網(wǎng)絡,編寫Matlab程序,采用滾動預測的方法對2015~2020年煤炭消費波動量進行預測,預測結(jié)果見表2。

    2.5煤炭消費量的預測

    采用經(jīng)驗模態(tài)分解法將煤炭消費量分解為煤炭消費趨勢量和波動量,并分別對二者進行了預測。煤炭消費趨勢量和波動量的預測結(jié)果相加求和即可得到煤炭消費總量的預測值,見表1。由表1可知,1982~2011年煤炭消費量擬合值的平均誤差為2.18%,2012~2014年煤炭消費量檢驗樣本預測值的平均誤差為1.24%,因此組合模型精度較高。

    為了證明組合模型的預測效果更佳,采用BP神經(jīng)網(wǎng)絡對1978~2014年煤炭消費量時間序列數(shù)據(jù)建立滾動預測模型,將時間序列分為訓練樣本和檢驗樣本。BP神經(jīng)網(wǎng)絡輸入層節(jié)點為4個,網(wǎng)絡結(jié)構(gòu)為4-4-1,訓練函數(shù)為“Trainlm”。 網(wǎng)絡訓練的參數(shù)為學習精度0.007、學習速率0.01、迭代步數(shù)1 000。訓練樣本為時間序列組成的前30組數(shù)據(jù),檢驗樣本為最后3組數(shù)據(jù),通過對BP網(wǎng)絡訓練精度和檢驗精度的控制,得到BP網(wǎng)絡滿足泛化能力要求的權(quán)值。用訓練好的BP網(wǎng)絡對煤炭消費量進行預測,擬合值、預測值及其相對誤差值見表1最后兩列。由表1可知,BP網(wǎng)絡滾動預測得到的1982~2011年煤炭消費量擬合值的平均誤差為4.18%,2012~2014年煤炭消費量檢驗樣本預測值的平均誤差為3.45%,因此BP網(wǎng)絡模型精度低于文中提出的基于經(jīng)驗模態(tài)分解的組合預測模型。故該組合預測模型可用于對2015~2020年煤炭消費量進行預測,預測值見表2。

    由表2可知,在今后一段時期內(nèi),煤炭消費量還將保持低速增長態(tài)勢,煤炭仍將是經(jīng)濟社會發(fā)展和人民生活水平提高中的重要能源之一。近年來,我國經(jīng)濟低速增長,宏觀經(jīng)濟進入轉(zhuǎn)型期,所以能源需求增長幅度不大,煤炭消費的增長速度漸緩,加上各城市環(huán)境污染問題日益突出,為改善環(huán)境質(zhì)量,煤炭消費面臨越來越大的壓力,煤炭行業(yè)由快速發(fā)展時期轉(zhuǎn)入低速發(fā)展期?!笆濉睍r期,我國經(jīng)濟社會發(fā)展呈現(xiàn)新趨勢,煤炭工業(yè)發(fā)展機遇與挑戰(zhàn)并存。一方面,從宏觀經(jīng)濟發(fā)展趨勢看,隨著我國工業(yè)化、信息化、城鎮(zhèn)化、農(nóng)業(yè)現(xiàn)代化持續(xù)推進,能源需求仍將保持增長,煤炭作為我國能源的主體地位不會改變。另一方面,我國經(jīng)濟發(fā)展進入新常態(tài),經(jīng)濟結(jié)構(gòu)不斷優(yōu)化,能源結(jié)構(gòu)不斷優(yōu)化,非化石能源比重上升,替代煤炭作用增強。同時,隨著科技進步,煤炭利用方式、利用領(lǐng)域不斷拓展,煤炭將由燃料向燃料和原料并重轉(zhuǎn)變。

    3 結(jié) 論

    (1)從提高煤炭消費量預測精度的目的出發(fā),針對煤炭消費量時間序列具有增長性、波動性和非平穩(wěn)性的特點,采用經(jīng)驗模態(tài)分解法分解出煤炭消費量時間序列中的趨勢量和波動量,由分解結(jié)果可知煤炭消費趨勢量增長趨勢明顯,波動性較小,并且其發(fā)展變化是影響煤炭消費總量變化的主要因素,煤炭消費波動量的波動性較大,其發(fā)展變化對煤炭消費總量的影響較小。

    (2)煤炭消費趨勢量占煤炭消費總量的比重較大,為了較準確的對煤炭消費趨勢量進行預測,分別建立了煤炭消費趨勢量預測的指數(shù)函數(shù)模型和一元線性回歸模型,并利用BP神經(jīng)網(wǎng)絡的非線性映射能力,綜合利用指數(shù)函數(shù)模型和一元線性回歸模型提供的有效信息,確定組合預測模型中各單一預測模型的權(quán)重,并利用建立的組合預測模型對煤炭消費趨勢量進行了預測。通過預測精度分析可知,組合預測模型的擬合平均誤差為2.7%,預測平均誤差為1.2%,模型具有較好的預測效果。

    (3)煤炭消費波動量的波動性較大,用傳統(tǒng)預測方法很難對其進行準確預測,因此利用BP神經(jīng)網(wǎng)絡具有逼近非線性函數(shù)的特點,采用滾動預測的方法,BP網(wǎng)絡結(jié)構(gòu)為4-4-1,對其進行預測。預測結(jié)果表明,模型擬合值的平均誤差為6.6%,預測值的平均誤差為3.7%,模型的預測精度稍低,但滿足預測精度要求。由于每年的煤炭消費波動量占煤炭消費總量的比重較小,因此其對煤炭消費總量的預測結(jié)果影響不大。

    (4)提出的基于經(jīng)驗模態(tài)分解的煤炭消費量組合預測模型,對1982~2014年煤炭消費總量進行了預測,1982-2011年煤炭消費量擬合值的平均誤差為2.18%,2012~2014年煤炭消費量預測值的平均誤差為1.24%,模型預測精度較高。

    (5)研究結(jié)果為煤炭消費量預測提供了一種新的方法,應用該方法對2015~2020年煤炭消費量進行了預測。預測結(jié)果表明,2015~2020年煤炭消費量將保持低速增長的態(tài)勢,研究結(jié)果為有關(guān)部門掌握煤炭消費量的短期發(fā)展趨勢,頒布煤炭工業(yè)發(fā)展的相關(guān)政策、措施等提供參考。

    [1]高峰.我國煤炭需求總量分析與預測[J].煤炭經(jīng)濟研究,2014, 34(4): 10-14.

    [2]劉滿芝,高曉峰,屈傳智,等. 中國煤炭需求波動規(guī)律研究[J].資源科學, 2013, 35(4): 681-689.

    [3]MUDIT K, JYOTI K P. Modeling demand for coal in India: vector autoregressive models with cointegrated variables[J]. Energy, 2000, 25(2): 149-168.

    [4]GUDARZI FARAHANI YAZDAN, VARMAZYARI BEHZAD, MOSHTARIDOUST SHIVA. Consumption in Iran: past trends and future directions[J]. Economics and Management, 2012, 62(1): 12-17.

    [5]CRISTINA CATTANEO, MATTEO MANERS, ELISA SCARPA. Industrial coal demand in China: A provincial analysis[J]. Resource and Energy Economics, 2011, 33(1): 12-35.

    [6]孫涵, 付曉靈, 張先鋒. 基于支持向量回歸機的中國煤炭長期需求預測[J].中國地質(zhì)大學學報, 2011, 11(9) :15-18.

    [7]楊俊祥, 程盛芳.灰色-周期外延組合模型在煤炭需求預測中的應用[J]. 統(tǒng)計與決策, 2010, (13): 162-163.

    [8]樊愛宛, 潘中強, 王巍. 灰色 GM(1,N)模型在河南省煤炭需求預測中的應用[J].煤炭技術(shù),2011, 30(10): 7-9.

    [9]丁宏飛, 黃福玲, 吳建樂. 基于GA-SVR 的煤炭需求預測模型研究[J].西南民族大學學報, 2010, 36(3): 402-404.

    [10]林伯強, 毛東昕. 煤炭消費終端部門對煤炭需求的動態(tài)影響分析[J].中國地質(zhì)大學學報, 2014, 14(6): 1-12.

    [11]鞠金艷, 趙林, 王金峰. 農(nóng)機總動力增長波動影響因素分析[J]. 農(nóng)業(yè)工程學報, 2016, 32(2): 84-91.

    [12]HUANG N E, SHEN Z, LONG S R, et al. The empirical mode decomposition and the Hilbert spectrum for nonlinear and non-stationary time series analysis[J]. Proc R Soc Land A,1998, 454(3): 899-955.

    [13]郭書琴. 基于經(jīng)驗模態(tài)分解的短期電力負荷預測[D].北京: 華北電力大學, 2010.

    [14]蔡吉花, 張世軍, 楊麗. EMD-SVM在南京市月平均氣溫預測中的應用, 2014, 44(22): 103-111.

    [15]鞠金艷, 王金武. 黑龍江省農(nóng)業(yè)機械化作業(yè)水平預測方法[J].農(nóng)業(yè)工程學報, 2009, 25(5): 83-88.

    [16]楊娟麗, 徐梅, 王福林, 等. 基于BP神經(jīng)網(wǎng)絡的時間序列預測問題研究[J].數(shù)學的實踐與認識, 2013, 43(4): 158-164.

    [17]鞠金艷. 黑龍江省農(nóng)業(yè)機械化發(fā)展的系統(tǒng)分析與對策研究[D].哈爾濱: 東北農(nóng)業(yè)大學, 2011.

    [18]AMINIAN F, SUAREZ E D.Forecasting economic data with neural networks[J].Computational Economies, 2006, 28(1): 71-88.

    [19]張洪潮, 王澤江, 李曉利, 等. 中國煤炭消費需求波動規(guī)律及成因分析[J].中國人口資源與環(huán)境,2014, 24(1): 94-101.

    (編輯徐巖)

    Combined prediction method of coal consumption based on empirical mode decomposition

    JUJinyan,ZHURongxin,CHENYu

    (School of Mechanical Engineering, Heilongjiang University of Science & Technology, Harbin 150022, China)

    This paper aims to investigate the law behind how coal consumption tends to develop in an effort to guide the healthy and orderly development of coal industry. In response to the growth, volatility, and instability inherent in the time series of coal consumption, the research using empirical mode decomposition method consists of the multi- level decomposition of the time series of coal consumption from 1978 to 2014 to obtain its development trend quantity and fluctuation quantity; the prediction of the trend quantity and fluctuation quantity using the nonlinear mapping ability of BP neural network and ultimate achievement of the forecasting results of coal consumption by summing the two forecast results. The error analysis shows that combined prediction model of coal consumption based on empirical mode decomposition gives an average error of 2.18% for fitted values and the average error of 1.24% for predicted values. The model enables an effective improvement in the forecast accuracy of coal consumption and works better for predicting coal consumption between 2015 and 2020. The prediction show that the coal consumption tends to keep a moderate growth in the coming years, and is expected to reach 3.417182 billion tce by 2020.

    coal consumption; empirical mode decomposition; combination forecast; BP network; time series

    2015-12-06

    鞠金艷(1982-),女,黑龍江省哈爾濱人,講師,博士,研究方向:生產(chǎn)管理與系統(tǒng)工程,E-mail:ju_jinyan@163.com。

    10.3969/j.issn.2095-7262.2016.01.024

    F407.21

    2095-7262(2016)01-0110-07

    A

    猜你喜歡
    消費量波動煤炭
    消除煤炭:是時候?qū)⒚禾枯d入史冊了 精讀
    英語文摘(2021年3期)2021-07-22 06:30:12
    煤炭
    羊肉價回穩(wěn) 后期不會大幅波動
    微風里優(yōu)美地波動
    中國化肥信息(2019年3期)2019-04-25 01:56:16
    煤炭:去產(chǎn)能繼續(xù) 為煤炭正名
    能源(2018年4期)2018-01-15 22:25:25
    干濕法SO2排放波動對比及分析
    2015年全球廢鋼消費量同比下降
    上海金屬(2016年4期)2016-04-07 16:43:41
    2014年國際橡膠消費量增加6.7%
    煤炭的“未來”
    能源(2015年8期)2015-05-26 09:15:34
    天堂中文最新版在线下载| 精品人妻一区二区三区麻豆| 日韩熟女老妇一区二区性免费视频| av在线观看视频网站免费| 2018国产大陆天天弄谢| 成年人午夜在线观看视频| 国产在线一区二区三区精| 91精品三级在线观看| 男的添女的下面高潮视频| 亚洲av.av天堂| 在线观看免费高清a一片| 婷婷成人精品国产| 亚洲中文av在线| 日本与韩国留学比较| 免费av中文字幕在线| 2021少妇久久久久久久久久久| 超碰97精品在线观看| 亚洲,欧美精品.| 亚洲激情五月婷婷啪啪| 国产日韩一区二区三区精品不卡| 欧美丝袜亚洲另类| 十八禁网站网址无遮挡| 日本wwww免费看| 91精品伊人久久大香线蕉| 女人被躁到高潮嗷嗷叫费观| 性高湖久久久久久久久免费观看| 成人毛片60女人毛片免费| 国产精品一国产av| 亚洲av.av天堂| 美女内射精品一级片tv| 国产亚洲av片在线观看秒播厂| 女性生殖器流出的白浆| 十分钟在线观看高清视频www| 大片免费播放器 马上看| 国产精品久久久久久精品古装| 亚洲精品一区蜜桃| 国产xxxxx性猛交| 久久久精品免费免费高清| 18+在线观看网站| 满18在线观看网站| 五月伊人婷婷丁香| 日韩不卡一区二区三区视频在线| 亚洲一级一片aⅴ在线观看| 大码成人一级视频| 久久热在线av| 水蜜桃什么品种好| 国产在线一区二区三区精| 国产免费又黄又爽又色| 国产不卡av网站在线观看| 精品酒店卫生间| 亚洲精品自拍成人| 99热网站在线观看| 久久ye,这里只有精品| 久久青草综合色| 五月天丁香电影| 欧美日韩亚洲高清精品| 国产亚洲av片在线观看秒播厂| 一区二区日韩欧美中文字幕 | 国产熟女欧美一区二区| 汤姆久久久久久久影院中文字幕| www.色视频.com| 男女边吃奶边做爰视频| 91精品伊人久久大香线蕉| 国产成人aa在线观看| 999精品在线视频| 午夜影院在线不卡| 久久精品国产亚洲av涩爱| 午夜福利在线观看免费完整高清在| 观看av在线不卡| 亚洲欧美日韩卡通动漫| a级毛色黄片| 赤兔流量卡办理| 日韩av在线免费看完整版不卡| 夫妻午夜视频| 国产精品久久久久久久久免| 国产在线免费精品| 亚洲精品乱久久久久久| 欧美xxxx性猛交bbbb| 国产又色又爽无遮挡免| 免费观看在线日韩| 美女主播在线视频| av一本久久久久| 咕卡用的链子| 香蕉精品网在线| 久久这里有精品视频免费| 欧美最新免费一区二区三区| 国产福利在线免费观看视频| 成人毛片60女人毛片免费| 日本午夜av视频| 97在线视频观看| 大片免费播放器 马上看| 啦啦啦啦在线视频资源| 精品一区二区三卡| 国产精品国产三级专区第一集| av一本久久久久| 天堂俺去俺来也www色官网| 日韩成人伦理影院| 精品一区二区免费观看| 国产 精品1| 宅男免费午夜| 菩萨蛮人人尽说江南好唐韦庄| 一区二区三区乱码不卡18| av在线老鸭窝| 免费人妻精品一区二区三区视频| 亚洲高清免费不卡视频| 成人午夜精彩视频在线观看| 日本午夜av视频| 嫩草影院入口| 两个人免费观看高清视频| 视频在线观看一区二区三区| 在线观看免费高清a一片| 乱码一卡2卡4卡精品| 日本av手机在线免费观看| 韩国av在线不卡| 91精品三级在线观看| 免费播放大片免费观看视频在线观看| av电影中文网址| 精品视频人人做人人爽| av播播在线观看一区| 啦啦啦中文免费视频观看日本| 黑人巨大精品欧美一区二区蜜桃 | 国内精品宾馆在线| 亚洲丝袜综合中文字幕| 99香蕉大伊视频| 美女内射精品一级片tv| 亚洲图色成人| 人妻一区二区av| 九九在线视频观看精品| 国产精品国产三级专区第一集| xxx大片免费视频| 国产视频首页在线观看| 婷婷色综合www| 精品99又大又爽又粗少妇毛片| 国产精品国产av在线观看| av.在线天堂| 亚洲国产精品专区欧美| 黄色视频在线播放观看不卡| 九草在线视频观看| 精品久久国产蜜桃| 亚洲精品国产av蜜桃| 99精国产麻豆久久婷婷| 久久热在线av| 国产一区有黄有色的免费视频| 纵有疾风起免费观看全集完整版| 中文字幕av电影在线播放| av黄色大香蕉| 国产精品成人在线| 日韩免费高清中文字幕av| 男女无遮挡免费网站观看| 超色免费av| 久久久久久久精品精品| 国产精品一区www在线观看| 纯流量卡能插随身wifi吗| h视频一区二区三区| 亚洲人成77777在线视频| 国产精品成人在线| 亚洲精品国产色婷婷电影| 成年人午夜在线观看视频| 男女午夜视频在线观看 | 夜夜骑夜夜射夜夜干| 国产午夜精品一二区理论片| 久久精品国产自在天天线| 欧美日韩国产mv在线观看视频| 青春草视频在线免费观看| 男女边摸边吃奶| 免费久久久久久久精品成人欧美视频 | 亚洲欧洲精品一区二区精品久久久 | 18禁观看日本| 国产成人精品无人区| 黄片播放在线免费| 日本欧美视频一区| 哪个播放器可以免费观看大片| 王馨瑶露胸无遮挡在线观看| 一区二区日韩欧美中文字幕 | 亚洲婷婷狠狠爱综合网| 亚洲五月色婷婷综合| 中文字幕免费在线视频6| 一级片'在线观看视频| 日本av手机在线免费观看| 天堂俺去俺来也www色官网| 午夜免费观看性视频| 亚洲av电影在线进入| av播播在线观看一区| 自拍欧美九色日韩亚洲蝌蚪91| 午夜91福利影院| 在线亚洲精品国产二区图片欧美| 久久久久久久久久成人| 久久99精品国语久久久| 美国免费a级毛片| 两个人免费观看高清视频| 26uuu在线亚洲综合色| 欧美 亚洲 国产 日韩一| 国产精品不卡视频一区二区| 国产精品人妻久久久影院| 天堂俺去俺来也www色官网| 99视频精品全部免费 在线| 久久久精品94久久精品| 在线观看人妻少妇| 69精品国产乱码久久久| 欧美精品亚洲一区二区| 九九在线视频观看精品| 亚洲 欧美一区二区三区| 日韩视频在线欧美| 99热网站在线观看| 女的被弄到高潮叫床怎么办| 日韩一本色道免费dvd| 一二三四中文在线观看免费高清| 久久毛片免费看一区二区三区| 亚洲欧洲国产日韩| 超碰97精品在线观看| 久久人妻熟女aⅴ| 色视频在线一区二区三区| 国产色爽女视频免费观看| 超色免费av| 国语对白做爰xxxⅹ性视频网站| 亚洲欧美成人精品一区二区| 国产成人精品婷婷| 精品亚洲乱码少妇综合久久| 日本色播在线视频| 99热6这里只有精品| 大陆偷拍与自拍| 99久久中文字幕三级久久日本| 9191精品国产免费久久| 亚洲精品中文字幕在线视频| 亚洲人成网站在线观看播放| 日韩大片免费观看网站| 久久午夜福利片| 亚洲精品国产av蜜桃| 成人免费观看视频高清| 搡老乐熟女国产| 97在线人人人人妻| 国产欧美日韩综合在线一区二区| 韩国高清视频一区二区三区| 纵有疾风起免费观看全集完整版| 一本久久精品| 国产成人精品福利久久| 国产男女内射视频| 一级爰片在线观看| 国产欧美日韩综合在线一区二区| 韩国av在线不卡| 肉色欧美久久久久久久蜜桃| 久久99热这里只频精品6学生| 欧美97在线视频| 亚洲美女黄色视频免费看| 亚洲av免费高清在线观看| 黄色配什么色好看| 视频中文字幕在线观看| 蜜桃在线观看..| 一本色道久久久久久精品综合| 亚洲精品456在线播放app| 欧美精品一区二区大全| 国内精品宾馆在线| 全区人妻精品视频| 9191精品国产免费久久| 亚洲中文av在线| 久久久久人妻精品一区果冻| 国产色婷婷99| 我的女老师完整版在线观看| 亚洲美女搞黄在线观看| 中文字幕人妻熟女乱码| 少妇 在线观看| 亚洲第一区二区三区不卡| 免费观看性生交大片5| 亚洲丝袜综合中文字幕| 亚洲av福利一区| 欧美成人午夜免费资源| av福利片在线| 黄色 视频免费看| 日韩一本色道免费dvd| 午夜福利影视在线免费观看| 亚洲经典国产精华液单| 国产1区2区3区精品| 欧美激情 高清一区二区三区| 啦啦啦在线观看免费高清www| 中文字幕免费在线视频6| 最新的欧美精品一区二区| 国产亚洲精品久久久com| 成年人免费黄色播放视频| 男女边吃奶边做爰视频| 男人添女人高潮全过程视频| 爱豆传媒免费全集在线观看| 丝瓜视频免费看黄片| 国产精品一国产av| 国产无遮挡羞羞视频在线观看| 精品福利永久在线观看| 精品一区二区三区视频在线| 亚洲国产精品一区二区三区在线| 久久国内精品自在自线图片| 在线天堂最新版资源| av国产精品久久久久影院| kizo精华| 中文字幕最新亚洲高清| 蜜桃在线观看..| 高清av免费在线| videos熟女内射| av福利片在线| 三上悠亚av全集在线观看| freevideosex欧美| 国产精品一区www在线观看| 欧美日韩一区二区视频在线观看视频在线| 亚洲国产最新在线播放| 亚洲成色77777| 啦啦啦中文免费视频观看日本| 亚洲精品一二三| 2018国产大陆天天弄谢| 男女无遮挡免费网站观看| 成人国产av品久久久| 美女福利国产在线| 免费大片黄手机在线观看| 久久午夜综合久久蜜桃| 亚洲精华国产精华液的使用体验| 国产精品免费大片| 午夜福利乱码中文字幕| 夫妻午夜视频| 亚洲人成网站在线观看播放| 亚洲精品国产色婷婷电影| 免费人成在线观看视频色| 少妇的丰满在线观看| 久久精品国产a三级三级三级| 亚洲熟女精品中文字幕| 欧美亚洲日本最大视频资源| 哪个播放器可以免费观看大片| 亚洲精品成人av观看孕妇| 国产女主播在线喷水免费视频网站| 午夜福利影视在线免费观看| 肉色欧美久久久久久久蜜桃| 国内精品宾馆在线| 国产福利在线免费观看视频| a级毛片黄视频| 天堂中文最新版在线下载| 午夜老司机福利剧场| 久久精品人人爽人人爽视色| 嫩草影院入口| 99国产综合亚洲精品| 王馨瑶露胸无遮挡在线观看| 久久精品人人爽人人爽视色| 国产免费一区二区三区四区乱码| 美女脱内裤让男人舔精品视频| 美女内射精品一级片tv| 精品国产一区二区三区久久久樱花| 91精品国产国语对白视频| 成年人免费黄色播放视频| 天美传媒精品一区二区| 国产在线视频一区二区| 国产免费视频播放在线视频| 国产在线免费精品| 精品卡一卡二卡四卡免费| 成人亚洲精品一区在线观看| 午夜91福利影院| 夫妻性生交免费视频一级片| 日韩欧美精品免费久久| 插逼视频在线观看| 午夜福利影视在线免费观看| 寂寞人妻少妇视频99o| 高清欧美精品videossex| 99国产综合亚洲精品| 亚洲国产最新在线播放| 狠狠婷婷综合久久久久久88av| 好男人视频免费观看在线| 女性生殖器流出的白浆| 成人无遮挡网站| 波野结衣二区三区在线| 人妻人人澡人人爽人人| 国产免费福利视频在线观看| 欧美日韩av久久| 久久久久久久久久久久大奶| 欧美日韩av久久| 伊人亚洲综合成人网| 国产成人aa在线观看| 天天躁夜夜躁狠狠久久av| 国产精品99久久99久久久不卡 | 成人免费观看视频高清| 多毛熟女@视频| 日日啪夜夜爽| 一级片免费观看大全| 欧美日韩视频高清一区二区三区二| 日韩视频在线欧美| 国产极品粉嫩免费观看在线| 高清欧美精品videossex| 熟女av电影| 99久久中文字幕三级久久日本| 少妇人妻精品综合一区二区| 国产在线免费精品| 一级片'在线观看视频| 男人舔女人的私密视频| 亚洲激情五月婷婷啪啪| 欧美日韩国产mv在线观看视频| 欧美激情 高清一区二区三区| 一本—道久久a久久精品蜜桃钙片| 热re99久久精品国产66热6| 婷婷色麻豆天堂久久| 少妇 在线观看| 欧美亚洲 丝袜 人妻 在线| 激情视频va一区二区三区| 国产精品一国产av| 欧美成人午夜免费资源| 一本色道久久久久久精品综合| 在线观看一区二区三区激情| 亚洲内射少妇av| 成人午夜精彩视频在线观看| 久久人人爽av亚洲精品天堂| 精品人妻熟女毛片av久久网站| 亚洲色图综合在线观看| 久久女婷五月综合色啪小说| 2021少妇久久久久久久久久久| 爱豆传媒免费全集在线观看| 国产淫语在线视频| 成人18禁高潮啪啪吃奶动态图| 欧美+日韩+精品| 日日爽夜夜爽网站| 一级毛片黄色毛片免费观看视频| 9热在线视频观看99| 国产永久视频网站| 青春草国产在线视频| 一区在线观看完整版| 免费黄网站久久成人精品| 精品国产一区二区久久| 一个人免费看片子| 欧美少妇被猛烈插入视频| 亚洲国产精品999| 久久97久久精品| 亚洲人成77777在线视频| 婷婷色综合www| 亚洲精品aⅴ在线观看| 久久精品夜色国产| 日本色播在线视频| 新久久久久国产一级毛片| 自拍欧美九色日韩亚洲蝌蚪91| 精品人妻在线不人妻| 成人18禁高潮啪啪吃奶动态图| 久久久久精品人妻al黑| 精品久久久精品久久久| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产成人av激情在线播放| 新久久久久国产一级毛片| 亚洲精品,欧美精品| 97人妻天天添夜夜摸| 午夜福利在线观看免费完整高清在| av卡一久久| 亚洲性久久影院| 国产精品人妻久久久影院| 丁香六月天网| av播播在线观看一区| 女的被弄到高潮叫床怎么办| 一本—道久久a久久精品蜜桃钙片| 国产综合精华液| 国产又爽黄色视频| 我的女老师完整版在线观看| 男人舔女人的私密视频| 成人黄色视频免费在线看| 青春草国产在线视频| 久久免费观看电影| 大香蕉97超碰在线| 美女脱内裤让男人舔精品视频| 大片免费播放器 马上看| 精品酒店卫生间| 日韩成人av中文字幕在线观看| kizo精华| 精品一区在线观看国产| 精品一区二区三卡| 亚洲精品自拍成人| 亚洲精品乱久久久久久| 国产欧美日韩一区二区三区在线| 只有这里有精品99| 老司机亚洲免费影院| 亚洲av电影在线进入| 日本爱情动作片www.在线观看| 成人午夜精彩视频在线观看| 成人国语在线视频| 高清毛片免费看| 又大又黄又爽视频免费| 99久久人妻综合| 美女国产视频在线观看| 亚洲图色成人| 亚洲美女视频黄频| 久久人人97超碰香蕉20202| 最黄视频免费看| 晚上一个人看的免费电影| 久久人人97超碰香蕉20202| 制服诱惑二区| 99热全是精品| 欧美日韩视频精品一区| 妹子高潮喷水视频| 国产极品天堂在线| 日韩欧美精品免费久久| 亚洲国产色片| 国产精品一区二区在线观看99| av在线播放精品| 亚洲在久久综合| 国产有黄有色有爽视频| 国产视频首页在线观看| 亚洲丝袜综合中文字幕| 国产亚洲最大av| 免费看不卡的av| www.熟女人妻精品国产 | 日韩制服骚丝袜av| 精品亚洲成国产av| 777米奇影视久久| 国产免费又黄又爽又色| 有码 亚洲区| 自拍欧美九色日韩亚洲蝌蚪91| videos熟女内射| 极品人妻少妇av视频| 成年av动漫网址| 亚洲欧美成人综合另类久久久| 中文字幕亚洲精品专区| freevideosex欧美| a级毛色黄片| 国产精品无大码| 亚洲国产欧美日韩在线播放| 亚洲精品乱久久久久久| 亚洲,欧美,日韩| 亚洲精品国产av成人精品| 欧美亚洲 丝袜 人妻 在线| 九九爱精品视频在线观看| 有码 亚洲区| 国产一区二区激情短视频 | 久久这里有精品视频免费| 久久青草综合色| 在线看a的网站| 香蕉精品网在线| 又粗又硬又长又爽又黄的视频| 亚洲成av片中文字幕在线观看 | 亚洲五月色婷婷综合| 精品国产一区二区久久| 久久久欧美国产精品| 精品国产国语对白av| 婷婷色麻豆天堂久久| 97在线人人人人妻| 国产日韩一区二区三区精品不卡| 99精国产麻豆久久婷婷| 人妻系列 视频| 街头女战士在线观看网站| 老司机影院成人| 黄网站色视频无遮挡免费观看| 久久久亚洲精品成人影院| 最近的中文字幕免费完整| 18在线观看网站| 建设人人有责人人尽责人人享有的| 这个男人来自地球电影免费观看 | 成年美女黄网站色视频大全免费| 久久精品国产综合久久久 | 我要看黄色一级片免费的| 国产女主播在线喷水免费视频网站| 国产精品国产三级国产专区5o| 国产成人免费观看mmmm| 亚洲成av片中文字幕在线观看 | 欧美亚洲日本最大视频资源| 久久人人爽人人爽人人片va| 日韩熟女老妇一区二区性免费视频| 国产黄色视频一区二区在线观看| 久久这里只有精品19| 午夜免费男女啪啪视频观看| 一区二区三区乱码不卡18| 亚洲精品,欧美精品| 国产精品久久久av美女十八| 亚洲 欧美一区二区三区| av网站免费在线观看视频| 久久精品国产亚洲av天美| 久久午夜综合久久蜜桃| 天堂中文最新版在线下载| 女人被躁到高潮嗷嗷叫费观| 97在线视频观看| 亚洲人成77777在线视频| 建设人人有责人人尽责人人享有的| 免费不卡的大黄色大毛片视频在线观看| 亚洲精华国产精华液的使用体验| 又黄又粗又硬又大视频| av免费观看日本| 精品人妻在线不人妻| 看免费av毛片| 欧美日韩综合久久久久久| 久久久久精品性色| 免费播放大片免费观看视频在线观看| 热99久久久久精品小说推荐| 亚洲,一卡二卡三卡| 欧美日韩一区二区视频在线观看视频在线| 黄色视频在线播放观看不卡| 国产一区二区在线观看av| 欧美激情 高清一区二区三区| 亚洲av免费高清在线观看| 男人操女人黄网站| 亚洲成人一二三区av| 久久人人爽人人爽人人片va| 久久精品人人爽人人爽视色| 黄网站色视频无遮挡免费观看| 亚洲精品色激情综合| 大香蕉久久网| 国产精品国产三级国产专区5o| av又黄又爽大尺度在线免费看| 伦理电影大哥的女人| 国产日韩欧美视频二区| 亚洲综合色网址| 夫妻性生交免费视频一级片| 亚洲av综合色区一区| 午夜久久久在线观看| 成年人免费黄色播放视频| 婷婷色综合www| 最新的欧美精品一区二区| 欧美 日韩 精品 国产| 26uuu在线亚洲综合色| 丁香六月天网| 黑人猛操日本美女一级片| 精品少妇黑人巨大在线播放| 9191精品国产免费久久| 交换朋友夫妻互换小说| 激情视频va一区二区三区| www.色视频.com| 日韩一区二区视频免费看| 青春草亚洲视频在线观看| 久久鲁丝午夜福利片| 9191精品国产免费久久| 美女大奶头黄色视频| 亚洲精品456在线播放app| 亚洲综合色惰|