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

    基于負(fù)荷分解的短期負(fù)荷預(yù)測(cè)傳遞函數(shù)模型

    2012-10-08 06:33:14萬(wàn)志宏文福拴
    關(guān)鍵詞:最低氣溫傳遞函數(shù)算子

    劉 敏,萬(wàn)志宏,文福拴

    (1.廣州供電局,廣東 廣州 510620;2.華南理工大學(xué) 電力學(xué)院,廣東 廣州 510640;3.浙江大學(xué) 電氣工程學(xué)院,浙江 杭州 310027)

    0 引言

    短期負(fù)荷預(yù)測(cè)是電力系統(tǒng)規(guī)劃和運(yùn)行的基礎(chǔ),關(guān)系到整個(gè)電力系統(tǒng)運(yùn)行的安全性和經(jīng)濟(jì)性[1]。負(fù)荷需求受用電結(jié)構(gòu)、氣候變化、星期類型、需求側(cè)管理等多方面因素的影響,這使得負(fù)荷變化呈現(xiàn)出非平穩(wěn)的隨機(jī)過(guò)程。因此,要對(duì)負(fù)荷進(jìn)行準(zhǔn)確預(yù)測(cè)是相當(dāng)困難的。到目前為止,已提出來(lái)多種負(fù)荷預(yù)測(cè)方法,如時(shí)間序列分析法[2~4]、相似日法[5]、人工神經(jīng)元網(wǎng)絡(luò)[6~8]和支持向量機(jī)[9]等。

    在上述方法中,基于時(shí)間序列分析的預(yù)報(bào)模型在很多領(lǐng)域包括電力系統(tǒng)短期負(fù)荷預(yù)測(cè)中得到了比較成功的應(yīng)用。到目前為止,應(yīng)用于電力系統(tǒng)短期負(fù)荷預(yù)測(cè)的時(shí)間序列模型主要為累積式自回歸-滑動(dòng)平均模型 (Auto-Regressive Integrated Moving Average,ARIMA)。由于 ARIMA模型只涉及負(fù)荷一個(gè)變量,無(wú)法描述各影響因素對(duì)負(fù)荷變化的影響,因此這種方法的預(yù)測(cè)準(zhǔn)確性仍不如人意。為考慮各種相關(guān)因素對(duì)負(fù)荷變化的影響,文獻(xiàn)[10]采用了傳遞函數(shù)模型分別對(duì)雙休日和工作日的負(fù)荷進(jìn)行建模,預(yù)測(cè)精度因此有所提高,但這種方法存在進(jìn)一步改善的空間。

    在此背景下,以文獻(xiàn)[10]的工作為基礎(chǔ),本文做了進(jìn)一步的研究,以改善負(fù)荷預(yù)測(cè)的準(zhǔn)確性?;舅悸肥鞘紫劝沿?fù)荷分解為周期分量和非周期分量,然而對(duì)這兩個(gè)分量分別進(jìn)行預(yù)測(cè)。具體地講,首先基于歷史負(fù)荷數(shù)據(jù),采用傳統(tǒng)分解法提取出負(fù)荷變化的周期分量,并得到剩下的非周期分量。在此基礎(chǔ)上,采用逐步回歸法首先從影響負(fù)荷非周期分量的因素中篩選出主要因素,然后發(fā)展了預(yù)測(cè)非周期分量的傳遞函數(shù)模型。最后,用廣東電力系統(tǒng)實(shí)際負(fù)荷數(shù)據(jù)說(shuō)明了所發(fā)展的短期負(fù)荷預(yù)測(cè)模型的實(shí)際效果。

    1 負(fù)荷分解模型

    電力系統(tǒng)負(fù)荷一般呈現(xiàn)出以周、日為周期的變化規(guī)律。所以,可以將負(fù)荷分解為周期分量和非周期分量,以提高預(yù)測(cè)精度??紤]到每天同一時(shí)刻的負(fù)荷是等時(shí)間間隔的 (相隔1天),這樣每天同一時(shí)刻的負(fù)荷就組成了等時(shí)間間隔的隨機(jī)時(shí)間序列。因此,可以采用時(shí)間序列理論[11]對(duì)負(fù)荷進(jìn)行分析。時(shí)間序列分析主要包括兩種模型:a.乘法模型;b.加法模型。這里采用時(shí)間序列分析中的加法模型對(duì)負(fù)荷數(shù)據(jù)進(jìn)行分解。這樣,對(duì)于1天24 h 96個(gè)時(shí)段 (每15 min一個(gè)時(shí)段)就有96個(gè)加法模型,模型的具體形式如式(1):

    式中:Pk,t,Dk,t,Sk,t和 Ik,t分別為歷史樣本中第 k 天中的時(shí)段t的負(fù)荷序列、負(fù)荷序列的趨勢(shì)分量、負(fù)荷序列的周期分量和負(fù)荷序列的不規(guī)則分量;n為歷史樣本天數(shù)。

    各分量的分離過(guò)程如下。

    (1)對(duì)于具有明顯非線性趨勢(shì)的原始序列,可以采用指數(shù)、對(duì)數(shù)或者多項(xiàng)式函數(shù)去除原始序列中的非線性。負(fù)荷的趨勢(shì)分量主要受經(jīng)濟(jì)增長(zhǎng)影響,考慮到短期內(nèi)經(jīng)濟(jì)增長(zhǎng)速度保持不變,因此這里采用線性回歸法確定趨勢(shì)分量Dk,t。

    (2)用原始負(fù)荷序列減去步驟 (1)所獲得的趨勢(shì)分量 Dk,t得到序列{Yk,t},該序列中含有周期分量和不規(guī)則分量??梢圆捎檬?(2)獲得周期分量:

    式中:Z為Pk,t周期分量的周期;N為歷史數(shù)據(jù)樣本中所包含的周期數(shù)。

    在提取各分量時(shí),先將負(fù)荷數(shù)據(jù)按照各個(gè)時(shí)刻分類,形成96個(gè)負(fù)荷子序列。之后,依次對(duì)96個(gè)序列提取各分量。這樣,序列{Pk,t}是以周為周期的,即Z=7。

    (3)用原始負(fù)荷序列減去步驟 (1)和 (2)所得的趨勢(shì)分量 Dk,t和周期分量 Sk,t,求出原始負(fù)荷序列的不規(guī)則分量 Ik,t。

    通過(guò)上述步驟可以提取出序列{Pk,t}的周期分量Sk,t,這樣就可以預(yù)測(cè)出待預(yù)測(cè)日時(shí)段 t所對(duì)應(yīng)的周期分量~St。然后,就只需對(duì)序列{Pk,t}去除周期分量后得到的序列{Lk,t}建立預(yù)測(cè)模型,下面將對(duì)此展開(kāi)討論。

    2 數(shù)學(xué)模型

    2.1 變量篩選

    由于日最高溫度和日最低溫度是影響電力負(fù)荷的重要因素,而采用前述方法提取出的周期分量Sk,t只與星期類型有關(guān),與日最高氣溫和日最低氣溫?zé)o關(guān)。因此,時(shí)刻t的負(fù)荷序列{Pk,t}去除周期分量Sk,t后得到的序列{Lk,t}一般會(huì)與日最高氣溫和日最低氣溫相關(guān)。基于這樣的考慮,在建立序列{Lk,t}的傳遞函數(shù)模型時(shí),應(yīng)考慮以日最高氣溫和日最低氣溫作為模型的輸入量。傳遞函數(shù)模型的輸出受到隨機(jī)事件、日最高氣溫和日最低氣溫的共同影響。為簡(jiǎn)化問(wèn)題,可以首先對(duì)模型的候選輸入量進(jìn)行篩選,剔除相關(guān)性較小的輸入變量。這里采用逐步回歸法[12,13]對(duì)輸入變量進(jìn)行篩選。在逐步回歸方法中,一般把輸入變量稱為自變量,把Lk,t稱為因變量。逐步回歸法每次只引入或剔除一個(gè)自變量。一個(gè)自變量是否被引入或被剔除則取決于其偏回歸平方和的F變量篩選檢驗(yàn)結(jié)果或校正決定系數(shù)的大小。逐步回歸法的步驟如下[12,13]:

    (1)計(jì)算每個(gè)自變量與因變量的相關(guān)系數(shù)r,并根據(jù)其絕對(duì)值由大到小對(duì)自變量排序。選擇相關(guān)系數(shù)絕對(duì)值最大的那個(gè)自變量進(jìn)行一元線性回歸,并對(duì)回歸方程進(jìn)行顯著性檢驗(yàn)。若回歸效果不顯著,則說(shuō)明自變量對(duì)因變量的影響不顯著。若顯著,則進(jìn)行步驟 (2)。

    (2)引入變量。按相關(guān)系數(shù)絕對(duì)值由大到小逐個(gè)將自變量引入到回歸方程,建立新的回歸方程。對(duì)引入的變量進(jìn)行偏回歸平方和的F檢驗(yàn)。若顯著,則引入該變量,否則就不引入。

    (3)剔除變量。剔除變量的過(guò)程與引入變量的過(guò)程相反。

    (4)重復(fù)步驟 (2)和 (3),直至自變量篩選結(jié)束為止。

    2.2 傳遞函數(shù)模型

    傳遞函數(shù)模型可以描述自變量與因變量的關(guān)系,因此其適于描述Lk,t與日最高氣溫、日最低氣溫的關(guān)系。為描述方便,給出以下定義:

    定義 1 若算子B可以使得:BLk,t= Lk-1,t,BnLk,t=Lk-n,t,則稱 B 為延遲算子。

    由時(shí)間序列理論可知,傳遞函數(shù)模型適用于平穩(wěn)時(shí)間序列。因此,在建立描述Lk,t的傳遞函數(shù)模型前,必須對(duì)日最高氣溫序列{Tk,1}、日最低氣溫序列{Tk,2} 和序列{Lk,t} 進(jìn)行平穩(wěn)化處理??紤]到序列{Lk,t}是負(fù)荷序列{Pk,t}去除周期分量得到的,因此可以認(rèn)為序列{Lk,t}不含有周期趨勢(shì)?;谏鲜隹紤],為獲得平穩(wěn)時(shí)間序列,可以對(duì)日最高氣溫序列{Tk,1}、日最低氣溫序列{Tk,2} 和序列{Lk,t}分別進(jìn)行差分運(yùn)算。若在這三個(gè)序列中,存在非平穩(wěn)序列,則對(duì)該序列進(jìn)行一階差分運(yùn)算。若經(jīng)一階差分運(yùn)算后,得到的序列仍為非平穩(wěn)序列,則對(duì)其繼續(xù)進(jìn)行一階差分運(yùn)算,直至序列平穩(wěn)為止。

    設(shè) {Lk,t} 經(jīng)過(guò)平穩(wěn)化處理后得到序列{lk,t},日最高氣溫序列{Tk,1}和日最低氣溫序列{Tk,2}經(jīng)過(guò)平穩(wěn)化處理后分別得到序列{Xk,1} 和{Xk,2},則Lk,t的傳遞函數(shù)模型可描述如式 (3):

    式中:lk,,t為模型的輸入時(shí)間序列;Xk,i為模型的輸出時(shí)間序列;Q為經(jīng)逐步回歸篩選后得到的因變量的個(gè)數(shù);αk,t為零均值白噪聲序列;bi為第i個(gè)輸入序列的時(shí)間延遲;θt(B)=1 - θt,1B - θt,2B2- … -θt,qBq為q階滑動(dòng)平均算子;φt(B) = 1 -φt,1B -φt,2B2- … - φt,pBp為p階自回歸算子;Ωi(B) =1 - Ωi,1B - Ωi,2B2- … - Ωi,γBγ為 γ 階自變量時(shí)間序列;Ei(B)=1 - Ei,1B - Ei,2B2- …- Ei,μBμ為 μ 階自變量時(shí)間序列。

    傳遞函數(shù)模型的建立過(guò)程如下[14,15]。

    (1)平穩(wěn)性檢驗(yàn)。前已述及,傳遞函數(shù)模型只是對(duì)平穩(wěn)時(shí)間序列適用,這樣在模型建立之前就必須對(duì)輸入序列和輸出序列進(jìn)行平穩(wěn)性檢驗(yàn)。這里采用單位根檢驗(yàn)法[11]進(jìn)行平穩(wěn)性檢驗(yàn),若輸入序列和響應(yīng)序列均為平穩(wěn)序列,則進(jìn)行步驟(2),否則對(duì)非平穩(wěn)序列進(jìn)行平穩(wěn)化處理。

    (2)預(yù)白噪化。對(duì)經(jīng)過(guò)適當(dāng)差分后的輸入序列和輸出序列利用 Pandit-Wu建模方法[2]建立相應(yīng)的自回歸-滑動(dòng)平均模型 (Auto-Regressive Moving Average,ARMA),以產(chǎn)生白噪聲序列 {νk,t}和{εk,i}:

    式中:Wt(B)=1 - Wt,1B - Wt,2B2- … - Wt,υBυ為υ 階自回歸算子;Ut(B) =1 - Ut,1B - Ut,2B2- …- Ut,ωBω為ω階滑動(dòng)平均算子;Γi(B) = 1 -Γi,1B - Γi,2B2- … - Γi,ρBρ為 ρ階自回歸算子;Hi(B) =1 - Hi,1B - Hi,2B2- … - Hi,λBλ為 λ 階滑動(dòng)平均算子。

    (3) 建立預(yù)測(cè)模型??疾煨蛄?{νk,t}和{εk,i}互相關(guān)函數(shù),建立傳遞函數(shù)模型:

    式中:{ηk,t}為回歸殘差序列。

    (4)建立殘差序列{ηt}的ARMA模型:

    3 算例與分析

    這里采用廣東電力系統(tǒng)2009年5月8日至2009年7月22日的負(fù)荷數(shù)據(jù)以及廣東省日最高氣溫和日最低氣溫?cái)?shù)據(jù)作為樣本數(shù)據(jù),對(duì)2009年7月23日的負(fù)荷進(jìn)行預(yù)測(cè)。短期負(fù)荷預(yù)測(cè)模型建立過(guò)程如下。

    (1)負(fù)荷分解,剔除周期分量。采用第1節(jié)介紹的負(fù)荷分解模型對(duì)各時(shí)段負(fù)荷組成的序列{Pk,t} 進(jìn)行分解,剔除周期分量 Sk,t,獲得剔除周期分量后的負(fù)荷序列{Lk,t}。

    (2)采用逐步回歸方法篩選出影響 {Lk,t}的主要因素 /變量。以{Lk,t}為指標(biāo)、日最高氣溫和日最低氣溫為候選因素進(jìn)行逐步回歸篩選?,F(xiàn)在以2009年5月8日至7月22日0點(diǎn)的、剔除周期分量后的負(fù)荷數(shù)據(jù){Lk,1}為例來(lái)描述篩選過(guò)程。計(jì)算得知{Lk,1}與日最高氣溫的相關(guān)系數(shù)為0.699 9,{Lk,1}與日最低氣溫的相關(guān)系數(shù)為0.439 7。這樣,在逐步回歸中,就先以日最高氣溫Tk,1為自變量、Lk,1為因變量進(jìn)行一元線性回歸。表1列出了逐步篩選過(guò)程中的一些量和相關(guān)的檢驗(yàn)指標(biāo),其中 F1為檢驗(yàn) Tk,1對(duì) Lk,1影響是否顯著的統(tǒng)計(jì)量,F(xiàn)2為檢驗(yàn) Tk,2對(duì) Lk,1影響是否顯著的統(tǒng)計(jì)量。

    表1 逐步回歸Tab.1 The stepwise regression procedure

    從表1可以看出,在一元線性回歸、引入和剔除過(guò)程中,統(tǒng)計(jì)量均大于 Fα(這里取 α =0.05)。因此,經(jīng)過(guò)逐步回歸后,日最高氣溫 Tk,1和日最低氣溫Tk,2均被引入到傳遞函數(shù)模型之中。

    (3)建立傳遞函數(shù)模型并進(jìn)行預(yù)測(cè)。以步驟(2)篩選出來(lái)的主要影響因素/變量作為模型的輸入變量,建立 {Lk,t}的傳遞函數(shù)模型,并根據(jù)傳遞函數(shù)模型得到待預(yù)測(cè)日時(shí)段t剔除周期分量后的負(fù)荷數(shù)據(jù)的預(yù)測(cè)值。

    (4)待預(yù)測(cè)日的負(fù)荷預(yù)測(cè)。從負(fù)荷周期分量中得到待預(yù)測(cè)日時(shí)段t負(fù)荷的周期分量,以及待預(yù)測(cè)日時(shí)段t的負(fù)荷預(yù)測(cè)值。

    短期負(fù)荷預(yù)測(cè)流程如圖1所示。

    圖1 短期負(fù)荷預(yù)測(cè)流程Fig.1 The procedure of short-term load forecasting

    圖2為求得的2009年7月23日的負(fù)荷周期分量,圖3為該日去除周期分量后的負(fù)荷預(yù)測(cè)值~Lt,圖4為該日負(fù)荷預(yù)測(cè)值與實(shí)際值。

    圖2 負(fù)荷周期分量Fig.2 The periodic component of loads

    圖3 去除周期分量的負(fù)荷預(yù)測(cè)值Fig.3 The forecasted load with the periodic quantity removed

    圖4 預(yù)測(cè)負(fù)荷與實(shí)際負(fù)荷比較Fig.4 Comparisons of forecasted loads and actual loads

    為說(shuō)明本文所發(fā)展的預(yù)測(cè)模型的準(zhǔn)確性,這里將本文模型與 ARIMA模型、文獻(xiàn)[10]中的模型進(jìn)行比較。這三種模型在各時(shí)段的預(yù)測(cè)值與負(fù)荷實(shí)際值見(jiàn)附錄A中的表A1。

    為了比較這三種模型,這里定義了如下日負(fù)荷預(yù)測(cè)準(zhǔn)確率:

    式中:Ei為預(yù)測(cè)日第i個(gè)點(diǎn)的相對(duì)誤差;A為日負(fù)荷預(yù)測(cè)準(zhǔn)確率。

    從表2可以看出本文模型的日負(fù)荷預(yù)測(cè)準(zhǔn)確率高于另外兩種模型。

    附錄A

    為更全面地說(shuō)明所發(fā)展的負(fù)荷預(yù)測(cè)模型的準(zhǔn)確性,表3列出了三種模型在2009年7月24日至7月30日這一周內(nèi)的日負(fù)荷預(yù)測(cè)準(zhǔn)確率,從中可見(jiàn)本文模型明顯優(yōu)于其它兩種模型。

    表2 三種模型預(yù)測(cè)準(zhǔn)確率比較Tab.2 Comparisons of load forecasting accuracy indices of three load forecasting models

    表3 三種預(yù)測(cè)模型在一周內(nèi)的預(yù)測(cè)準(zhǔn)確率比較Tab.3 Comparisons of load forecasting accuracy indices of three load forecasting models for one-week period

    4 結(jié)論

    本文發(fā)展了一種基于負(fù)荷分解的短期負(fù)荷預(yù)測(cè)傳遞函數(shù)模型。采用廣東電力系統(tǒng)實(shí)際負(fù)荷數(shù)據(jù)所做短期負(fù)荷預(yù)測(cè)結(jié)果表明,所發(fā)展的方法能夠取得較現(xiàn)有的其它兩種方法更為準(zhǔn)確的負(fù)荷預(yù)測(cè)結(jié)果。所發(fā)展的負(fù)荷預(yù)測(cè)模型具有下述優(yōu)點(diǎn):

    (1)該模型能較好反映各種相關(guān)因素對(duì)負(fù)荷變化的影響。

    (2)采用負(fù)荷分解模型將負(fù)荷分解為周期和非周期分量,有效的提高了負(fù)荷預(yù)測(cè)的精度。

    (3)采用逐步回歸法篩選主要影響因素,也有助于減小負(fù)荷預(yù)測(cè)誤差。

    表A1 7月23日短期負(fù)荷預(yù)測(cè)結(jié)果Tab.A1 Results of load forecasting for July 23

    續(xù)表

    續(xù)表

    續(xù)表

    [1]羅滇生,姚建剛,何洪英,等.基于自適應(yīng)滾動(dòng)優(yōu)化的電力負(fù)荷多模型組合預(yù)測(cè)系統(tǒng)的研究與開(kāi)發(fā)[J].中國(guó)電機(jī)工程學(xué)報(bào),2003,23(5):58-61.

    [2]楊叔子,吳雅,軒建平,等.時(shí)間序列分析的工程應(yīng)用[M].武漢:華中科技大學(xué)出版社,2007.

    [3]BOX G E P,JENKINS G M.時(shí)間序列分析預(yù)測(cè)與控制[M].中國(guó)統(tǒng)計(jì)出版社,1997.

    [4]Huang S J,Shih K R.Short-term load forecasting via ARMA modelidentification including non-Gaussian process considerations[J].IEEE Trans on Power Systems,2003,18(2):673 -679.

    [5]莫維仁,張伯明,孫宏斌,等.短期負(fù)荷預(yù)測(cè)中選擇相似日的探討[J].清華大學(xué)學(xué)報(bào) (自然科學(xué)版),2004,44(1):106-109.

    [6]Hipper H S,Pefreira C E,Souza R C.Neural network for short-term load forecasting:a review and evaluation[J].IEEE Trans on Power Systems,2001,16(2):44-54.

    [7]HiroyukoI M.Deterministic annealing clustering for ANN-termload forecasting.IEEE Trans on Power Systems,2001,16(3):545 -551.

    [8]Taylor J W,Roberto B.Neural network load forecasting with weather ensemble predictions[J].IEEE Trans on Power Systems,2001,17(3):626 -632.

    [9]康重慶,夏清,劉梅.電力系統(tǒng)負(fù)荷預(yù)測(cè)[M].北京:中國(guó)電力出版社,2007.

    [10]李妮,江岳春,黃珊,等.基于累積式自回歸移動(dòng)平均傳遞函數(shù)模型的短期負(fù)荷預(yù)測(cè)[J].電網(wǎng)技術(shù),2009,33(8):93-97.

    [11]潘紅宇.金融時(shí)間序列模型[M].北京:對(duì)外經(jīng)濟(jì)貿(mào)易大學(xué)出版社,2008.

    [12]何曉群.現(xiàn)代統(tǒng)計(jì)分析方法與應(yīng)用[M].北京:中國(guó)人民大學(xué)出版社,1998.

    [13]方開(kāi)泰,金輝,陳慶云.實(shí)用回歸分析[M].北京:科學(xué)出版社,1988.

    [14]SAS/ETS User’s Guide,vision 5 edition[M].SAS Institute Inc.,Cary,NC,USA,1984.

    [15]高惠璇 編譯.SAS系統(tǒng) SAS/ETS軟件使用手冊(cè)[M].北京:中國(guó)統(tǒng)計(jì)出版社,1988.

    猜你喜歡
    最低氣溫傳遞函數(shù)算子
    擬微分算子在Hp(ω)上的有界性
    各向異性次Laplace算子和擬p-次Laplace算子的Picone恒等式及其應(yīng)用
    北辰地區(qū)日最低氣溫特征及影響因素研究
    天津科技(2019年3期)2019-03-30 07:17:14
    67年來(lái)濟(jì)南最高和最低氣溫變化特征
    一類Markov模算子半群與相應(yīng)的算子值Dirichlet型刻畫
    基于LabVIEW的火焰?zhèn)鬟f函數(shù)測(cè)量系統(tǒng)
    Roper-Suffridge延拓算子與Loewner鏈
    紅外成像系統(tǒng)的調(diào)制傳遞函數(shù)測(cè)試
    數(shù)學(xué)選擇題練習(xí)
    基于傳遞函數(shù)自我優(yōu)化的BP網(wǎng)絡(luò)算法改進(jìn)
    免费观看无遮挡的男女| 国产亚洲一区二区精品| 久久久久网色| 亚洲精品第二区| 亚洲va在线va天堂va国产| 中文字幕av电影在线播放| 日韩成人av中文字幕在线观看| 丰满乱子伦码专区| 国产国拍精品亚洲av在线观看| 亚洲精品色激情综合| 久久精品国产自在天天线| 男女边摸边吃奶| 男的添女的下面高潮视频| 色网站视频免费| av天堂中文字幕网| 国产av精品麻豆| 亚洲精品日韩在线中文字幕| 久久免费观看电影| 免费av不卡在线播放| 黄色配什么色好看| 伦理电影大哥的女人| 少妇高潮的动态图| 波野结衣二区三区在线| 国产一区二区三区综合在线观看 | 久久人妻熟女aⅴ| 在线天堂最新版资源| 国产黄片美女视频| 夫妻性生交免费视频一级片| 黄色日韩在线| 在线看a的网站| 人妻夜夜爽99麻豆av| 99视频精品全部免费 在线| av.在线天堂| 欧美精品人与动牲交sv欧美| a级毛片在线看网站| 欧美亚洲 丝袜 人妻 在线| 日韩一本色道免费dvd| 久久久午夜欧美精品| 国产黄频视频在线观看| 一区二区三区精品91| freevideosex欧美| 国产精品偷伦视频观看了| 视频区图区小说| 成人国产麻豆网| 中国三级夫妇交换| 七月丁香在线播放| 国产成人精品婷婷| 久久综合国产亚洲精品| 国产高清有码在线观看视频| 久久6这里有精品| 亚洲va在线va天堂va国产| 国产精品熟女久久久久浪| 国产片特级美女逼逼视频| 超碰97精品在线观看| 香蕉精品网在线| 国产91av在线免费观看| 中文在线观看免费www的网站| 亚洲三级黄色毛片| 伊人久久精品亚洲午夜| 在线天堂最新版资源| 日产精品乱码卡一卡2卡三| 亚洲美女视频黄频| av又黄又爽大尺度在线免费看| 久久人人爽人人爽人人片va| 精品少妇久久久久久888优播| 简卡轻食公司| 中文欧美无线码| 少妇人妻一区二区三区视频| 视频区图区小说| av又黄又爽大尺度在线免费看| 欧美另类一区| 一区二区三区四区激情视频| 啦啦啦视频在线资源免费观看| 久久国产亚洲av麻豆专区| 99久国产av精品国产电影| 曰老女人黄片| 狂野欧美白嫩少妇大欣赏| 久热久热在线精品观看| 秋霞在线观看毛片| 纯流量卡能插随身wifi吗| 嘟嘟电影网在线观看| 丁香六月天网| 人妻 亚洲 视频| 亚洲国产日韩一区二区| 黑丝袜美女国产一区| 欧美精品人与动牲交sv欧美| 国产精品无大码| 亚洲欧洲国产日韩| 精品久久久久久久久av| 成人18禁高潮啪啪吃奶动态图 | 婷婷色麻豆天堂久久| 自拍偷自拍亚洲精品老妇| 免费观看a级毛片全部| 亚洲av.av天堂| 在线观看国产h片| 大片免费播放器 马上看| 国产精品熟女久久久久浪| 狂野欧美激情性bbbbbb| 国产91av在线免费观看| av黄色大香蕉| 成年女人在线观看亚洲视频| 午夜福利,免费看| 岛国毛片在线播放| 成人亚洲欧美一区二区av| 日韩av不卡免费在线播放| 最新的欧美精品一区二区| 看非洲黑人一级黄片| 99热网站在线观看| 日韩强制内射视频| 免费观看性生交大片5| 欧美 亚洲 国产 日韩一| 麻豆成人午夜福利视频| 亚洲欧美日韩另类电影网站| 国产深夜福利视频在线观看| tube8黄色片| 在线观看一区二区三区激情| 内射极品少妇av片p| 最近中文字幕2019免费版| 亚洲精品自拍成人| 九九久久精品国产亚洲av麻豆| 伦理电影免费视频| 嫩草影院入口| 日韩欧美精品免费久久| 99热全是精品| 亚洲av综合色区一区| 国产精品麻豆人妻色哟哟久久| 国产精品99久久久久久久久| 久久国产精品大桥未久av | 精品久久久精品久久久| 久久精品久久久久久久性| 老熟女久久久| 18+在线观看网站| 免费大片黄手机在线观看| 91久久精品电影网| 26uuu在线亚洲综合色| 亚洲国产av新网站| 日本与韩国留学比较| 亚洲欧美精品专区久久| 国产精品福利在线免费观看| 女人精品久久久久毛片| 亚洲国产精品一区三区| 能在线免费看毛片的网站| 国产老妇伦熟女老妇高清| 韩国av在线不卡| 精品亚洲成国产av| 一个人看视频在线观看www免费| 国产免费福利视频在线观看| 少妇被粗大猛烈的视频| 国产欧美另类精品又又久久亚洲欧美| 黄色怎么调成土黄色| 国产片特级美女逼逼视频| 精品99又大又爽又粗少妇毛片| 精品一品国产午夜福利视频| 老司机影院毛片| 亚洲精品,欧美精品| 国产亚洲91精品色在线| 啦啦啦中文免费视频观看日本| 99久国产av精品国产电影| 中文字幕亚洲精品专区| 少妇被粗大的猛进出69影院 | 18禁裸乳无遮挡动漫免费视频| 亚洲av电影在线观看一区二区三区| 免费观看的影片在线观看| 晚上一个人看的免费电影| 赤兔流量卡办理| 亚洲欧美精品自产自拍| 日本午夜av视频| 国产精品免费大片| 国产在视频线精品| 男人狂女人下面高潮的视频| 久久精品熟女亚洲av麻豆精品| 欧美激情极品国产一区二区三区 | 国产精品女同一区二区软件| 亚洲综合精品二区| 欧美丝袜亚洲另类| 国产淫语在线视频| 黄色日韩在线| 亚洲熟女精品中文字幕| 男女国产视频网站| videossex国产| 热re99久久精品国产66热6| 亚洲丝袜综合中文字幕| 日本免费在线观看一区| 国内揄拍国产精品人妻在线| 人人澡人人妻人| 下体分泌物呈黄色| 自线自在国产av| 色5月婷婷丁香| 伊人亚洲综合成人网| 日韩大片免费观看网站| 免费久久久久久久精品成人欧美视频 | 一个人免费看片子| 国产又色又爽无遮挡免| h视频一区二区三区| 一个人看视频在线观看www免费| 一级,二级,三级黄色视频| 久久久久久久久大av| 在线观看美女被高潮喷水网站| 男人爽女人下面视频在线观看| 久久精品国产鲁丝片午夜精品| 亚洲精华国产精华液的使用体验| 久久99一区二区三区| 女人精品久久久久毛片| 天美传媒精品一区二区| 国产中年淑女户外野战色| 中文字幕免费在线视频6| 99国产精品免费福利视频| 亚洲精品aⅴ在线观看| 18禁裸乳无遮挡动漫免费视频| 高清不卡的av网站| 久久99热6这里只有精品| 美女中出高潮动态图| 亚洲成人av在线免费| 婷婷色综合大香蕉| 黑人猛操日本美女一级片| videossex国产| 女的被弄到高潮叫床怎么办| 在线播放无遮挡| 交换朋友夫妻互换小说| 久久久久视频综合| 亚洲婷婷狠狠爱综合网| 中文字幕制服av| 午夜福利,免费看| 国产极品天堂在线| 色哟哟·www| 麻豆成人午夜福利视频| 国产高清三级在线| 国产av精品麻豆| 国产欧美日韩综合在线一区二区 | 中文乱码字字幕精品一区二区三区| 色视频www国产| 免费黄频网站在线观看国产| 少妇高潮的动态图| 国产精品偷伦视频观看了| 男女国产视频网站| 黄色视频在线播放观看不卡| 欧美97在线视频| 全区人妻精品视频| 极品人妻少妇av视频| 啦啦啦中文免费视频观看日本| 免费看av在线观看网站| 久久毛片免费看一区二区三区| 亚洲欧美一区二区三区黑人 | 欧美国产精品一级二级三级 | 亚洲欧美精品自产自拍| 国产黄色免费在线视频| 91久久精品电影网| 你懂的网址亚洲精品在线观看| 在线观看国产h片| 热re99久久精品国产66热6| 成人毛片a级毛片在线播放| 成人毛片60女人毛片免费| 亚洲精品乱码久久久久久按摩| 亚洲自偷自拍三级| 日本欧美国产在线视频| 99久久人妻综合| 波野结衣二区三区在线| 又大又黄又爽视频免费| 在线观看免费高清a一片| 在线亚洲精品国产二区图片欧美 | 国产精品一区二区三区四区免费观看| av在线观看视频网站免费| 伦理电影免费视频| 熟女人妻精品中文字幕| 欧美97在线视频| 大话2 男鬼变身卡| 十分钟在线观看高清视频www | 91精品国产九色| 欧美激情国产日韩精品一区| 国产极品天堂在线| 成人18禁高潮啪啪吃奶动态图 | 天美传媒精品一区二区| 国产精品国产三级国产专区5o| 另类亚洲欧美激情| 亚洲精品亚洲一区二区| 欧美精品国产亚洲| 国产片特级美女逼逼视频| 精品亚洲成国产av| 色网站视频免费| 成人国产av品久久久| 精品国产一区二区三区久久久樱花| 国产 一区精品| 青青草视频在线视频观看| 在线免费观看不下载黄p国产| 国产在线男女| 日韩欧美 国产精品| 午夜激情久久久久久久| 性色avwww在线观看| 99热这里只有是精品50| 肉色欧美久久久久久久蜜桃| 日本wwww免费看| 午夜91福利影院| 一级爰片在线观看| 午夜福利在线观看免费完整高清在| 国产成人免费观看mmmm| 午夜福利视频精品| 国产精品伦人一区二区| 麻豆成人av视频| 亚洲精品,欧美精品| 国产黄片视频在线免费观看| 韩国av在线不卡| 大又大粗又爽又黄少妇毛片口| 国产免费一级a男人的天堂| 少妇的逼好多水| 一级毛片我不卡| 中文字幕久久专区| 边亲边吃奶的免费视频| 一级a做视频免费观看| 国产 精品1| 国产成人精品一,二区| 国产女主播在线喷水免费视频网站| 中文精品一卡2卡3卡4更新| 中文乱码字字幕精品一区二区三区| h视频一区二区三区| 亚洲精品亚洲一区二区| 色视频在线一区二区三区| 一级av片app| 国产高清不卡午夜福利| 九九久久精品国产亚洲av麻豆| 精品一区二区三区视频在线| 日韩一区二区视频免费看| 一区二区三区精品91| 久久青草综合色| 精品99又大又爽又粗少妇毛片| 成人亚洲欧美一区二区av| 日韩三级伦理在线观看| 午夜福利视频精品| 亚洲精品乱码久久久v下载方式| 久久人妻熟女aⅴ| 欧美精品国产亚洲| 麻豆乱淫一区二区| 国产在线免费精品| 日日撸夜夜添| 一本一本综合久久| 久久久欧美国产精品| 免费观看a级毛片全部| 一本久久精品| 久久狼人影院| 久久久午夜欧美精品| 插阴视频在线观看视频| 日韩欧美一区视频在线观看 | 国产淫片久久久久久久久| 亚洲国产日韩一区二区| 又粗又硬又长又爽又黄的视频| a级毛色黄片| 久久久久久久久久人人人人人人| 五月玫瑰六月丁香| 免费在线观看成人毛片| 99精国产麻豆久久婷婷| 国产成人91sexporn| 最近手机中文字幕大全| 蜜桃在线观看..| 99热国产这里只有精品6| 成年人午夜在线观看视频| 内射极品少妇av片p| 国产欧美另类精品又又久久亚洲欧美| 狂野欧美激情性bbbbbb| 不卡视频在线观看欧美| 99久久中文字幕三级久久日本| 女性被躁到高潮视频| 日韩,欧美,国产一区二区三区| av又黄又爽大尺度在线免费看| 另类精品久久| 波野结衣二区三区在线| 大片电影免费在线观看免费| 天堂俺去俺来也www色官网| 夫妻性生交免费视频一级片| 激情五月婷婷亚洲| 熟妇人妻不卡中文字幕| 亚洲成色77777| 2018国产大陆天天弄谢| 欧美丝袜亚洲另类| 国产高清不卡午夜福利| 久久ye,这里只有精品| 丰满饥渴人妻一区二区三| 国产免费又黄又爽又色| 如日韩欧美国产精品一区二区三区 | 国产一区有黄有色的免费视频| 精品久久久久久久久亚洲| 久久国内精品自在自线图片| 国产精品99久久久久久久久| 丰满人妻一区二区三区视频av| 免费黄色在线免费观看| 黄色一级大片看看| 亚洲欧洲日产国产| 99热这里只有精品一区| av网站免费在线观看视频| 卡戴珊不雅视频在线播放| 黄色欧美视频在线观看| 欧美日韩在线观看h| 午夜91福利影院| 欧美人与善性xxx| 啦啦啦在线观看免费高清www| 婷婷色av中文字幕| 欧美三级亚洲精品| 黄色怎么调成土黄色| 国产日韩欧美在线精品| 国产精品偷伦视频观看了| 插阴视频在线观看视频| 免费看不卡的av| 日日撸夜夜添| 妹子高潮喷水视频| 欧美xxⅹ黑人| 国产av国产精品国产| 99久久精品一区二区三区| 26uuu在线亚洲综合色| 水蜜桃什么品种好| 肉色欧美久久久久久久蜜桃| 观看美女的网站| 久久精品熟女亚洲av麻豆精品| 成人国产麻豆网| 青春草国产在线视频| 肉色欧美久久久久久久蜜桃| 中国国产av一级| 又爽又黄a免费视频| 免费人妻精品一区二区三区视频| av国产久精品久网站免费入址| 成人黄色视频免费在线看| 黄色怎么调成土黄色| 日韩视频在线欧美| 一区在线观看完整版| 人妻少妇偷人精品九色| 欧美激情极品国产一区二区三区 | 欧美变态另类bdsm刘玥| 极品人妻少妇av视频| 妹子高潮喷水视频| 男女免费视频国产| 免费观看无遮挡的男女| 人妻 亚洲 视频| 少妇人妻久久综合中文| 久久精品国产亚洲网站| 精品国产一区二区三区久久久樱花| 人人澡人人妻人| 日本欧美视频一区| 九九久久精品国产亚洲av麻豆| 一级av片app| 亚洲欧美一区二区三区黑人 | 国产美女午夜福利| 亚洲av二区三区四区| 噜噜噜噜噜久久久久久91| 日韩视频在线欧美| 亚洲一级一片aⅴ在线观看| 免费黄色在线免费观看| 国产伦理片在线播放av一区| 午夜激情久久久久久久| 久久久久久人妻| 欧美日韩视频精品一区| 午夜老司机福利剧场| 欧美变态另类bdsm刘玥| 国产精品秋霞免费鲁丝片| 秋霞在线观看毛片| 欧美老熟妇乱子伦牲交| 亚洲国产精品专区欧美| 纯流量卡能插随身wifi吗| 亚洲欧美精品专区久久| 一级黄片播放器| 久久久亚洲精品成人影院| 亚洲四区av| 免费av不卡在线播放| 国产精品久久久久久av不卡| 91精品一卡2卡3卡4卡| 亚洲一级一片aⅴ在线观看| 中文字幕免费在线视频6| 少妇人妻一区二区三区视频| 熟女人妻精品中文字幕| 国产熟女欧美一区二区| 好男人视频免费观看在线| 看免费成人av毛片| 久久ye,这里只有精品| 免费少妇av软件| 日韩大片免费观看网站| 久久女婷五月综合色啪小说| 亚洲av中文av极速乱| 精品少妇内射三级| 日韩欧美一区视频在线观看 | 亚洲精品中文字幕在线视频 | 一级a做视频免费观看| 亚洲第一av免费看| 免费不卡的大黄色大毛片视频在线观看| 亚洲精品国产av蜜桃| 日本与韩国留学比较| 日韩 亚洲 欧美在线| 日韩人妻高清精品专区| 男人和女人高潮做爰伦理| 亚洲电影在线观看av| 99九九在线精品视频 | 综合色丁香网| 人妻一区二区av| 我要看日韩黄色一级片| 久久久欧美国产精品| 精品久久国产蜜桃| 晚上一个人看的免费电影| 一区二区三区精品91| 国产 精品1| 三级国产精品欧美在线观看| 卡戴珊不雅视频在线播放| 麻豆成人av视频| 中国国产av一级| 18+在线观看网站| 日韩av免费高清视频| 免费少妇av软件| 成人国产av品久久久| 在线天堂最新版资源| 免费观看无遮挡的男女| 久久久久久久久久成人| a级毛片免费高清观看在线播放| 秋霞伦理黄片| 又粗又硬又长又爽又黄的视频| 国产精品麻豆人妻色哟哟久久| 亚洲av.av天堂| 日日啪夜夜爽| 日韩中字成人| 香蕉精品网在线| 街头女战士在线观看网站| 久久精品国产亚洲网站| 男女边摸边吃奶| 少妇熟女欧美另类| a级毛片在线看网站| 亚洲国产精品成人久久小说| 国产免费福利视频在线观看| 最后的刺客免费高清国语| 女性生殖器流出的白浆| 综合色丁香网| 人妻一区二区av| 免费大片18禁| 另类亚洲欧美激情| 国产精品三级大全| 男女无遮挡免费网站观看| 久久久久视频综合| 久久久久久久久久久丰满| 免费播放大片免费观看视频在线观看| 亚洲国产欧美日韩在线播放 | 日韩免费高清中文字幕av| 欧美最新免费一区二区三区| 高清不卡的av网站| 欧美激情极品国产一区二区三区 | 欧美区成人在线视频| 卡戴珊不雅视频在线播放| 人人妻人人澡人人爽人人夜夜| 国产精品久久久久久av不卡| 三级经典国产精品| 热re99久久国产66热| 少妇高潮的动态图| 少妇猛男粗大的猛烈进出视频| 亚洲成色77777| 一二三四中文在线观看免费高清| 久久久久久久久久久免费av| 亚洲成人手机| 精品国产一区二区三区久久久樱花| 大片免费播放器 马上看| 精品少妇黑人巨大在线播放| 观看美女的网站| h日本视频在线播放| 在线观看免费高清a一片| 精品国产一区二区三区久久久樱花| 一本大道久久a久久精品| 高清黄色对白视频在线免费看 | 高清午夜精品一区二区三区| 99久国产av精品国产电影| 天天操日日干夜夜撸| 韩国av在线不卡| 欧美亚洲 丝袜 人妻 在线| 建设人人有责人人尽责人人享有的| 国产精品成人在线| 热re99久久精品国产66热6| 哪个播放器可以免费观看大片| 亚洲国产成人一精品久久久| 国产精品成人在线| 18禁动态无遮挡网站| 久久精品国产自在天天线| 十八禁高潮呻吟视频 | av黄色大香蕉| 黑人高潮一二区| 秋霞伦理黄片| 亚洲av电影在线观看一区二区三区| 成年女人在线观看亚洲视频| 中文精品一卡2卡3卡4更新| 国产av一区二区精品久久| 在线观看国产h片| 亚洲在久久综合| 精品少妇久久久久久888优播| 国产在线视频一区二区| 伦理电影免费视频| 欧美激情国产日韩精品一区| 国产黄片视频在线免费观看| 女性生殖器流出的白浆| 亚洲国产最新在线播放| 成年美女黄网站色视频大全免费 | 国产精品99久久99久久久不卡 | a级片在线免费高清观看视频| 男女国产视频网站| 国产成人freesex在线| 黄色毛片三级朝国网站 | 99九九在线精品视频 | 免费观看的影片在线观看| 日本黄大片高清| 三级国产精品片| 婷婷色综合大香蕉| 中文在线观看免费www的网站| 大片电影免费在线观看免费| 午夜免费观看性视频| 伦理电影免费视频| 日韩,欧美,国产一区二区三区| 好男人视频免费观看在线| 午夜福利在线观看免费完整高清在| 水蜜桃什么品种好| 亚洲av国产av综合av卡| 热re99久久精品国产66热6| 久久久久精品久久久久真实原创| 久久久a久久爽久久v久久| 亚洲内射少妇av| 国产成人一区二区在线| 啦啦啦视频在线资源免费观看| 熟女电影av网| 精品人妻熟女av久视频|