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

    基于改進馬爾科夫鏈的風電功率時間序列模型

    2017-07-10 10:26:20趙宇肖白顧兵王逍祎張鈺王成龍
    電力建設 2017年7期
    關鍵詞:馬爾科夫電功率時刻

    趙宇,肖白,顧兵,王逍祎,張鈺,王成龍

    (東北電力大學,吉林省吉林市 132000)

    基于改進馬爾科夫鏈的風電功率時間序列模型

    趙宇,肖白,顧兵,王逍祎,張鈺,王成龍

    (東北電力大學,吉林省吉林市 132000)

    模擬風電功率時間序列在風電并網(wǎng)系統(tǒng)的規(guī)劃和評估研究中具有重要意義,針對原始馬爾科夫鏈在風電功率建模上無法保留其自相關性的不足,構建了一種基于改進馬爾科夫鏈的風電功率時間序列模型。首先分析了風電功率的季節(jié)特性、日特性和波動特性;然后將風電功率數(shù)據(jù)按照不同月份及時段進行了細致劃分,生成相應的狀態(tài)轉移概率矩陣;最后,對風電功率波動量的概率分布進行擬合,并疊加波動量,建立了基于改進馬爾科夫鏈的風電功率時間序列模型。實例分析表明,本文所建新模型生成的風電功率序列能夠保留歷史序列自相關性,同時在一般統(tǒng)計參數(shù)、概率密度分布和自相關性三方面的準確性也優(yōu)于已有模型。

    馬爾科夫鏈;風電功率;時間序列模擬;波動特性

    0 引 言

    能源緊缺與環(huán)境危機問題促使了以風力發(fā)電為代表的新能源發(fā)電技術飛快發(fā)展[1-3]。由于風力發(fā)電的隨機性與不確定性,使得風電并網(wǎng)后電力系統(tǒng)穩(wěn)定、安全、可靠運行面臨巨大挑戰(zhàn),因此研究風電功率的特性,并生成風電功率模擬序列,在電力系統(tǒng)規(guī)劃及安全性評估等領域具有重要意義。模擬生成風電功率序列是指以歷史觀測數(shù)據(jù)為基礎,生成多個在統(tǒng)計特征上與歷史序列相吻合的新序列。在目前研究中生成風電功率序列的方法可分為兩類:風速法和風功率法[4-5]。

    風速法是以風速時間序列作為輸入量,利用風速-功率轉換模型生成風電功率序列。受尾流效應、地形、能量轉換曲線等因素影響,風速法生成的模擬風電功率序列較歷史觀測序列存在較大誤差。風功率法直接利用觀測數(shù)據(jù)生成風電功率時間序列,避免了經(jīng)過風速-功率轉換而進一步造成的偏差。文獻[6-7]在生成風電功率序列時分別采用了自回歸滑動模型(autoregressive integrated moving average,ARIMA)和馬爾科夫鏈蒙特卡洛法(Markov chain Monte Carlo method,MCMC),這2種方法在概率密度函數(shù)和自相關性上都表現(xiàn)出較好的性能,且后者優(yōu)于前者。文獻[8]提出了3種離散馬爾科夫模型,分別為基本馬爾科夫模型、貝葉斯馬爾科夫模型和birth-and-death馬爾科夫模型,相比之下貝葉斯馬爾科夫模型生成的風電功率序列具有更好的自相關性。但這些方法均沒有考慮風電功率自身特性,所以導致生成的序列與原始序列相比均存在較大誤差。

    本文提出一種改進的馬爾科夫鏈模型,根據(jù)風電功率的原始數(shù)據(jù),綜合考慮風電功率的季節(jié)特性、日特性及波動特性直接生成風電功率時間序列,既避免了構建風速-功率轉換模型帶來的偏差又保留序列自身的特點[9-17]。最后基于某實際風電場數(shù)據(jù)進行仿真。其結果表明,改進馬爾科夫鏈模型生成的風電功率序列與原始馬爾科夫鏈模型生成的序列相比具有更高的準確性。

    1 馬爾科夫鏈模型

    馬爾科夫鏈是指離散的馬爾科夫過程,即隨機模擬過程當中下一時刻的狀態(tài)與過去狀態(tài)無關。

    1.1 離散狀態(tài)

    風力發(fā)電的每一個輸出功率值都可看做是一個特定的狀態(tài),根據(jù)精度要求的不同可以設定不同的離散狀態(tài)數(shù),進而將風電功率劃分為不同的出力范圍。假設風電輸出功率區(qū)間為(amin,amax),設定S個狀態(tài)后,每個狀態(tài)對應的風電功率范圍大小為(amin,amax)/S。

    1.2 基本原理

    假設在一階馬爾科夫鏈X={X1,X2,…,Xt,…,XN}中,Xt為模擬時間序列中的元素。其中,N為模擬時間序列的數(shù)據(jù)個數(shù),風電功率的值用S個狀態(tài)構成的狀態(tài)空間表示為E={E1,E2,…,ES},每一時刻的功率值只可對應其中1個狀態(tài),因此每一個狀態(tài)都有S種轉移方向(包括轉向自身),在t時刻風電功率值處于狀態(tài)Ei的前提下,在t+1時刻轉向狀態(tài)Ej的條件概率為

    pij=Pr(Xt+1=Ej|Xt=Ei)

    (1)

    式中:Xt+1、Xt分別為t+1時刻和t時刻風電功率所處的狀態(tài);pij(i,j=1,2,…,S)為從t時刻的狀態(tài)Ei轉移到t+1時刻的狀態(tài)Ej的概率;Pr為概率函數(shù)。

    所有狀態(tài)轉移概率的集合形成的矩陣稱為狀態(tài)轉移概率矩陣P。P為S階方陣,且各行元素和為1,其形式為

    (2)

    狀態(tài)轉移概率矩陣是馬爾科夫鏈模型模擬生成風電功率時間序列的基礎。

    2 考慮季節(jié)特性、日特性及波動特性的改進一階馬爾科夫模型

    2.1 季節(jié)特性

    風電功率的季節(jié)特性主要表現(xiàn)為1年中不同月份輸出功率的大小存在差異。為了在生成的風電功率時間序列中體現(xiàn)出季節(jié)特性,本文將1年的數(shù)據(jù)分為12個片段,用λ表示,λ∈{1,2,…,12},分別對應1年中的12個月份,對每一個片段單獨建立狀態(tài)轉移概率矩陣。

    2.2 日特性

    風電功率的日特性主要表現(xiàn)為由于日照的原因,1天當中不同時段風電功率的大小存在差異,且不同月份的日特性一般也不同。在電力系統(tǒng)規(guī)劃中,考慮風電功率的日特性在研究負荷和風電之間的影響,以及風電與光伏互補發(fā)電中具有重要作用。本文將1天的數(shù)據(jù)分為4個部分,并分別用θi表示。θi∈{θ1,θ2,θ3,θ4}。θ1對應00:00—06:00,θ2對應06:00—12:00,θ3對應12:00—18:00,θ4對應18:00—24:00。

    2.3 波動特性

    由于風速的間歇性及不確定性,導致了風電輸出功率存在一定波動,其波動的特點對電力系統(tǒng)的穩(wěn)定運行、規(guī)劃及調(diào)度都會帶來不利影響,因此在研究電力系統(tǒng)安全評估及長期規(guī)劃中都必須考慮這種波動特性。風電功率的波動性是指后一時刻與當前時刻的風電功率值存在差異,其波動量有不同的定義方式,本文采用一階差分量(即相鄰2個單位時間內(nèi)的功率差值)來刻畫:

    ΔP=Pi-Pi-1,i=2,3,…,N

    (3)

    式中:ΔP為風電功率的波動量;Pi為第i個單位時間內(nèi)風電功率的平均值;N為模擬時間序列的數(shù)據(jù)個數(shù)。

    根據(jù)文獻[12]的研究成果,選取t location-scale分布函數(shù)對風電功率波動量的概率密度函數(shù)進行擬合:

    (4)

    式中:μ為位置參數(shù);σ為尺度參數(shù);v為形狀參數(shù)。

    2.4 改進一階馬爾科夫鏈模型

    改進的馬爾科夫鏈模型可分為2個部分:其一是考慮風電功率季節(jié)特性、日特性形成的多個狀態(tài)轉移概率矩陣;其二是生成滿足概率分布函數(shù)的波動量的集合。模型建立步驟如下詳述。

    (1)設定狀態(tài)數(shù)S,并根據(jù)歷史風電功率序列Ω,將風電功率的取值區(qū)間均分為S個部分;

    (2)將風電功率序列Ω根據(jù)相應的月份λ和時段θ進行分類,每一類對應的序列記為Ωλ,θ;

    (3)求取序列Ωλ,θ中狀態(tài)i到j的轉移概率:

    (5)

    式中:pij,λ,θ為序列Ωλ,θ中狀態(tài)i到j的轉移概率;nij,λ,θ為序列Ωλ,θ中從狀態(tài)i到狀態(tài)j的轉移頻數(shù)(只經(jīng)過1步轉移)。

    對應的Ωλ,θ序列的狀態(tài)轉移概率矩陣Pλ,θ為

    (6)

    (4)根據(jù)式(3)計算出歷史風電功率序列的波動量,進而用t location-scale分布函數(shù)對波動量的概率密度進行擬合,得到參數(shù)μ、σ、v,然后生成符合該參數(shù)的波動值的集合,從集合中隨機生成1個波動值作為下一時刻相對于當前時刻的波動量。

    3 基于改進馬爾科夫鏈的風電功率序列模擬

    應用改進馬爾科夫鏈模型模擬生成風電功率時間序列的具體步驟如下詳述。

    (1)假設當前時刻為t,所處狀態(tài)為Et,風電功率大小為αt,所處月份為λt,所處時段為θt,生成1個服從均勻分布的隨機數(shù)ε∈(0,1);

    (2)根據(jù)已知的λt、θt,選取相應的狀態(tài)轉移矩陣Pλt,θt,進而得出累計狀態(tài)轉移概率矩陣Qλt,θt:

    (7)

    Qλt,θt中元素的取值如下:

    (8)

    式中qij,t為t時刻狀態(tài)i到狀態(tài)j的累計轉移概率。

    (3)假設風電功率下一時刻的狀態(tài)為Et+1,若 0<ε

    (4) 確定了下一時刻模擬風電功率所處狀態(tài)Et+1后,還需進一步確定該狀態(tài)的具體風電功率值αt+1,在符合t location scale分布函數(shù)的參數(shù)的波動值集合中抽取波動量γt,令η=αt+γt,并判斷η是否在狀態(tài)Et+1的取值范圍內(nèi),若是,則αt+1=η,否則重新抽取波動量再進行判斷;

    (5) 判斷t是否為結束時刻,若是,則結束,若不是,則令t=t+1,返回步驟(2)。

    具體流程如圖1所示。

    圖1 改進模型序列生成流程圖Fig.1 Flow chart of improved model sequences generation

    4 算例分析

    本節(jié)利用西北地區(qū)某風電場的觀測數(shù)據(jù),計算并對比觀測風電功率序列與原始馬爾科夫鏈模型和改進馬爾科夫鏈模型模擬生成的風電功率序列的統(tǒng)計特性,其中一般統(tǒng)計參數(shù)包括平均值、標準差、最大值。概率分布特性用概率密度函數(shù)(probability density function,PDF)來衡量,序列的自相關性用自相關函數(shù)(autocorrelation function,ACF)來衡量。

    4.1 模型建立

    根據(jù)觀測風電功率的區(qū)間,將其平均分為40個狀態(tài),其中第i個狀態(tài)對應的功率區(qū)間為[(i-1),i]× 2 MW。按照2.4節(jié)給出的方法建立改進一階馬爾科夫鏈模型,3月06:00—12:00對應的狀態(tài)轉移矩陣P3,θ2為

    (9)

    由于觀測數(shù)據(jù)的時間間隔為5 min,因此生成的模擬序列也是每隔5 min的數(shù)據(jù),生成的風電功率部分序列如圖2所示。為了能夠呈現(xiàn)清晰、直觀的對比,只給出前750 min的風電功率序列。

    4.2 一般統(tǒng)計參數(shù)比較

    將改進馬爾科夫鏈模型和原始馬爾科夫鏈模型生成的風電功率序列按照不同的月份分為12類,分別計算每類的平均值、標準差和最大值,以及各統(tǒng)計參數(shù)的絕對誤差和相對誤差,并與原始風電功率序列的參數(shù)值進行對比,結果見表1—3。

    西北地區(qū)屬大陸性氣候,春季平均風速最大,冬夏次之,秋季最小。從表1中可看出,原始馬爾科夫鏈生成的風電功率序列中,每個月的功率平均值并沒有太大變動,且與原始序列的平均值相差較大,而改進馬爾科夫鏈生成的模擬序列不僅在平均值上更接近于原始序列,且每個月份風電功率平均值的變化情況也符合該地區(qū)風電場輸出功率的季節(jié)特性。由此可知,改進馬爾科夫鏈模型在考慮季節(jié)特性、日特性和波動特性之后擁有更高的精度。表2和表3與表1的分析結果類似。

    圖2 部分歷史序列、改進模型序列及原始模型序列Fig.2 Part of historical sequence, improved model sequence and primitive model sequence

    表2 不同月份風電功率標準差比較Table 2 Comparison of standard deviation of wind power in different months

    表3 不同月份風電功率最大值比較Table 3 Comparison of maximum wind power in different months

    4.3 PDF與ACF的比較

    改進馬爾科夫鏈模型和原始馬爾科夫鏈模型生成的風電功率序列與歷史觀測序列的PDF對比如圖3所示。

    圖3 不同風電功率序列PDF對比Fig.3 PDF comparison of different wind power sequences

    從圖3的對比中可得出,改進馬爾科夫鏈模型生成的風電功率序列較好地保留了歷史觀測序列的概率分布特性,而原始馬爾科夫鏈模型生成的風電功率序列在概率分布上與原始序列相比相差很大。

    原始風電功率序列、原始模型序列與改進模型序列,在不同滯后時間下的ACF如圖4—6所示。

    圖4 原始風電功率序列ACFFig.4 ACF of original wind power sequence

    圖5 原始馬爾科夫鏈模型生成的風電功率序列ACFFig.5 ACF of wind power sequence by original Markov chain model

    通過比較可以看出,改進馬爾科夫鏈生成的風電功率序列基本保留了歷史序列的自相關性,原始馬爾科夫鏈生成的風電功率序列的ACF與歷史序列相比相差較大,這是由于風電功率在1天中不同時段的特性不同。改進馬爾科夫鏈引入了日特性,故生成的序列保留了良好的自相關性。

    圖6 改進馬爾科夫鏈模型生成的風電功率序列ACFFig.6 ACF of wind power sequence by improved Markov chain model

    5 結 論

    本文基于馬爾科夫鏈模型直接利用歷史數(shù)據(jù)生成風電功率序列,在原始一階馬爾科夫鏈的基礎上考慮了風電功率的季節(jié)特性、日特性及波動特性,生成不同的狀態(tài)轉移矩陣,并疊加波動量,得出改進一階馬爾科夫鏈模型。經(jīng)過實例驗證,改進模型在較好地保留原始序列概率分布特性和自相關特性的同時,還能提高風電功率序列的模擬精度,為研究風電并網(wǎng)后的電力系統(tǒng)可靠性分析工作奠定了基礎。

    [1]CHAI C I, LEE W J, FUANGFOO P, et al. System impact study for the interconnection of wind generation and utility system[J]. IEEE Transactions on Industry Applications, 2005, 41(1): 163-168.

    [2]張麗英,葉廷路,辛耀中,等.大規(guī)模風電接入電網(wǎng)的相關問題及措施[J].中國電機工程學報,2010,30(25):1-9. ZHANG Liying, YE Tinglu, XIN Yaozhong, et al. Problems and measures of power grid accommodating large scale wind power[J]. Proceedings of the CSEE, 2010, 30(25): 1-9.

    [3]李軍軍,吳政球,譚勛瓊,等.風力發(fā)電及其技術發(fā)展綜述[J].電力建設,2011,32(8):64-72. LI Junjun, WU Zhengqiu, TAN Xunqiong, et al. Review of wind power generation and relative technology development[J]. Electric Power Construction, 2011, 32(8):64-72.

    [4]孟巖峰,胡書舉,鄧雅,等.風電功率預測誤差分析及預測誤差評價方法[J].電力建設,2013,34(7):6-9. MENG Yanfeng, HU Shuju, DENG Ya, et al. Analysis and evaluation method of wind power predicted-error[J]. Electric Power Construction, 2013, 34(7):6-9.

    [5]汪寧渤,馬彥宏,王健東.大規(guī)模風電集中并網(wǎng)對電力系統(tǒng)安全穩(wěn)定的影響[J].電力建設,2011,32(11):77-80. WANG Ningbo, MA Yanhong, WANG Jiandong. Analysis of power system security and stability caused by large-scale wind power grid integration[J]. Electric Power Construction, 2011, 32(11): 77-80.

    [6]CHEN P, PEDERSEN T, BAK-JENSEN B, et al. ARIMA-based time series model of stochastic wind power generation [J]. IEEE Transactions on Power System, 2010, 25(2): 667-676.

    [7]PAPAEFTHYMIOU G, KLOCKL B. MCMC for wind power simulation[J]. IEEE Transactions on Energy Conversion, 2007, 23(1): 234-240.

    [8]CHEN P, BERTHELSEN K K, BAK-JENSEN B, et al. Markov model of wind power time series using bayesian inference of ransition matrix[C]//35th Annual Conference of IEEE on Industrial Electronics. Porto: IEEE, 2009: 627-632.

    [9]蔣平,霍雨翀,張龍,等.基于改進一階馬爾可夫鏈的風速時間序列模型[J].電力系統(tǒng)自動化,2014,38(19):22-27. JANG Ping, HUO Yuchong, ZHANG long, et al. Wind speed time series model based on improved first-order markov chain[J]. Automation of Electric Power Systems, 2014, 38(19): 22-27.

    [10]丁明,鮑玉瑩,畢銳.應用改進馬爾科夫鏈的光伏出力時間序列模擬[J].電網(wǎng)技術,2016,40(2):459-464. DING Ming, BAO Yuying, BI Rui. Simulation of PV output time series used improved markov chain[J]. Power System Technology, 2016, 40(2): 459-464.

    [11]于鵬,黎靜華,文勁宇,等.含風電功率時域特性的風電功率序列建模方法[J],中國電機工程學報,2014,34(22):3715-3723. YU Peng, LI Jinghua, WEN Jinyu, et al. A wind power time series modeling method based on its time domain characteristics[J]. Proceedings of the CSEE, 2014, 34(22): 3715-3723.

    [12]林衛(wèi)星,文勁宇,艾小猛,等.風電功率波動特性的概率分布研究[J].中國電機工程學報,2012,32(1):38-46. LIN Weixing, WEN Jinyu, AI Xiaomeng, et al. Probability density function of wind power variations[J]. Proceedings of the CSEE, 2012, 32(1): 38-46.

    [13]丁明,徐寧舟.基于馬爾可夫鏈的光伏發(fā)電系統(tǒng)輸出功率短期預測方法[J].電網(wǎng)技術,2011,35(1):152-157. DING Ming, XU Ningzhou. Short-term prediction method of output power of photovoltaic power generation system based on markov chain[J]. Power System Technology, 2011, 35(1): 152-157.

    [14]LI Y Z, LUAN R, NIU J C. Forecast of power generation for grid-connected photovoltaic system based on grey model and markov chain[C]//3rd IEEE Conference on Industrial Electronics and Applications. Singapore: IEEE, 2008: 1729-1733.

    [15]張遠實,潘延林.風電功率模擬生成模型研究[J].電網(wǎng)與清潔能源,2014,30(12):126-138. ZHANG Yuanshi, PAN Yanlin. Studies on wind power curves simulaton via a stochastic process[J]. Power System and Clean Energy, 2014, 30(12): 126-138.

    [16]張宏宇,印永華,申洪,等.基于概率測度變換的風速時間序列建模方法[J].電力系統(tǒng)自動化,2013,37(2):7-10. ZHANG Hongyu, YIN Yonghua, SHEN Hong, et al. Modeling of wind speed time series based on probability measure transformation[J]. Automation of Electric Power Systems, 2013, 37(2): 7-10.

    [17]羅鋼,石東源,陳金富,等.風光發(fā)電功率時間序列模擬的MCMC方法[J].電網(wǎng)技術,2014,38(2):311-327. LUO Gang, SHI Dongyuan, CHEN Jinfu, et al. A markov chain monte carlo method for simulation of wind and solar power time series[J]. Power System Technology, 2014, 38(2): 311-327.

    (編輯 郭文瑞)

    Wind Power Time Series Model Based on Improved Markov Chain

    ZHAO Yu,XIAO Bai,GU Bing,WANG Xiaoyi,ZHANG Yu,WANG Chenglong

    (Northeast Electric Power University, Jilin 132000, Jilin Province, China)

    The simulation of wind power time series is of great significance in the planning and evaluation of wind power grid-connected systems. In order to solve the problem that the original Markov chain can not keep its autocorrelation in wind power modeling, this paper constructs a wind power time series model based on improved Markov chain. This paper firstly analyzes the seasonal characteristics, daily characteristics and fluctuation characteristics of wind power; and then subdivides the wind power data according to different months and time periods to generate the corresponding state transition probability matrix. Finally, this paper fits the probability distribution of wind power fluctuation and increase the amount of fluctuation to establish the wind power time series model based on improved Markov chain. The case analysis shows that the wind power series generated by the proposed model is superior to the existing model in the aspects of general statistical parameters, probability density distribution and autocorrelation, while preserving the historical sequence autocorrelation.

    Markov chain; wind power; time series simulation; fluctuation characteristics

    TM 614

    A

    1000-7229(2017)07-0018-07

    10.3969/j.issn.1000-7229.2017.07.003

    2017-02-11

    趙宇(1992),女,碩士研究生,主要從事含風力發(fā)電的電力系統(tǒng)規(guī)劃方面的研究工作;

    肖白(1973),男,博士,教授,主要從事電力系統(tǒng)規(guī)劃、空間負荷預測、城市電網(wǎng)風險評估和電力系統(tǒng)繼電保護等方面的研究工作;

    顧兵(1979),女,碩士,副教授,主要從事全壽命周期理論在電力系統(tǒng)中的應用等方面的研究工作;

    王逍祎(1993),男,碩士研究生,主要從事風電并網(wǎng)技術方面的研究工作;

    張鈺(1996),女,碩士研究生,主要從事含風力發(fā)電的電力系統(tǒng)規(guī)劃方面的研究工作;

    王成龍(1994),男,碩士研究生,主要從事含風力發(fā)電的電力系統(tǒng)規(guī)劃方面的研究工作。

    猜你喜歡
    馬爾科夫電功率時刻
    基于PCC-CNN-GRU的短期風電功率預測
    冬“傲”時刻
    捕獵時刻
    基于疊加馬爾科夫鏈的邊坡位移預測研究
    輕松上手電功率
    你會計算電功率嗎
    基于改進的灰色-馬爾科夫模型在風機沉降中的應用
    解讀電功率
    馬爾科夫鏈在教學評價中的應用
    街拍的歡樂時刻到來了
    亚洲av中文字字幕乱码综合| 久久久久免费精品人妻一区二区| 免费看十八禁软件| 五月伊人婷婷丁香| 亚洲成人免费电影在线观看| 欧美黄色淫秽网站| 成人午夜高清在线视频| avwww免费| 在线a可以看的网站| 精品一区二区三区av网在线观看| 国模一区二区三区四区视频| 18禁国产床啪视频网站| 老汉色∧v一级毛片| netflix在线观看网站| 亚洲精品乱码久久久v下载方式 | 十八禁人妻一区二区| 美女黄网站色视频| 此物有八面人人有两片| 国产精品乱码一区二三区的特点| www国产在线视频色| 欧美xxxx黑人xx丫x性爽| 久久久久久久久久黄片| 一进一出抽搐gif免费好疼| 免费无遮挡裸体视频| 看片在线看免费视频| 国产单亲对白刺激| 性色av乱码一区二区三区2| 免费av毛片视频| 日韩,欧美,国产一区二区三区| 美女被艹到高潮喷水动态| 男女下面进入的视频免费午夜| 黄色配什么色好看| 精品久久久噜噜| 成人av在线播放网站| 一个人观看的视频www高清免费观看| 看非洲黑人一级黄片| 久久久久久久久久成人| kizo精华| 国产免费又黄又爽又色| 国产淫片久久久久久久久| 国产精品日韩av在线免费观看| 亚洲第一区二区三区不卡| 国产成人精品一,二区| 三级毛片av免费| 搡老乐熟女国产| 哪个播放器可以免费观看大片| 免费大片黄手机在线观看| 亚洲国产成人一精品久久久| 久久6这里有精品| 日本色播在线视频| 人人妻人人澡人人爽人人夜夜 | 亚洲欧洲国产日韩| 夫妻性生交免费视频一级片| 精品久久久久久成人av| 成人无遮挡网站| 九九久久精品国产亚洲av麻豆| 国产午夜福利久久久久久| 中文在线观看免费www的网站| 国语对白做爰xxxⅹ性视频网站| 亚洲av福利一区| 22中文网久久字幕| 国产伦一二天堂av在线观看| 亚洲美女搞黄在线观看| 卡戴珊不雅视频在线播放| 菩萨蛮人人尽说江南好唐韦庄| 亚洲av日韩在线播放| 国产在线男女| 午夜精品国产一区二区电影 | 美女cb高潮喷水在线观看| 亚洲国产精品专区欧美| 欧美激情国产日韩精品一区| 欧美成人午夜免费资源| 在现免费观看毛片| 亚洲丝袜综合中文字幕| 久99久视频精品免费| 亚洲av一区综合| 久久久久久久久大av| 又粗又硬又长又爽又黄的视频| 久久久精品94久久精品| 啦啦啦中文免费视频观看日本| 丰满人妻一区二区三区视频av| 高清午夜精品一区二区三区| 日韩一区二区视频免费看| 晚上一个人看的免费电影| 一区二区三区乱码不卡18| 精品久久国产蜜桃| 深夜a级毛片| 成人欧美大片| 中文字幕免费在线视频6| 国语对白做爰xxxⅹ性视频网站| av一本久久久久| 五月伊人婷婷丁香| 国产精品一区二区性色av| 久久久国产一区二区| 99re6热这里在线精品视频| 国产亚洲一区二区精品| 有码 亚洲区| 亚洲成人精品中文字幕电影| 白带黄色成豆腐渣| 成人综合一区亚洲| 18禁动态无遮挡网站| 国产淫片久久久久久久久| 在线观看美女被高潮喷水网站| 一区二区三区高清视频在线| 亚洲av中文字字幕乱码综合| 亚洲精品乱码久久久v下载方式| 男人和女人高潮做爰伦理| 夫妻午夜视频| 国产午夜精品久久久久久一区二区三区| 国产精品一区二区三区四区免费观看| 国产 一区 欧美 日韩| 少妇熟女欧美另类| 成人av在线播放网站| 午夜爱爱视频在线播放| 18禁动态无遮挡网站| 韩国av在线不卡| 大片免费播放器 马上看| av在线天堂中文字幕| 久久久久久久久久黄片| 少妇丰满av| 精品熟女少妇av免费看| 中文字幕av在线有码专区| 美女内射精品一级片tv| 久久久国产一区二区| 日日撸夜夜添| 日韩av在线免费看完整版不卡| 免费av不卡在线播放| 秋霞伦理黄片| av在线老鸭窝| 欧美极品一区二区三区四区| 美女cb高潮喷水在线观看| 欧美zozozo另类| 中国国产av一级| 97超视频在线观看视频| 哪个播放器可以免费观看大片| 成人毛片a级毛片在线播放| 搞女人的毛片| 亚洲三级黄色毛片| 蜜桃久久精品国产亚洲av| 国产成人aa在线观看| 亚洲综合精品二区| 啦啦啦中文免费视频观看日本| 热99在线观看视频| av在线观看视频网站免费| 成人毛片a级毛片在线播放| 99视频精品全部免费 在线| 特级一级黄色大片| 亚洲av成人精品一区久久| 青春草国产在线视频| 美女xxoo啪啪120秒动态图| 欧美日韩一区二区视频在线观看视频在线 | ponron亚洲| www.av在线官网国产| 国产又色又爽无遮挡免| 国产v大片淫在线免费观看| 精品久久久久久久末码| 街头女战士在线观看网站| 日本熟妇午夜| 久久精品熟女亚洲av麻豆精品 | 久热久热在线精品观看| 亚洲高清免费不卡视频| 一级片'在线观看视频| 91午夜精品亚洲一区二区三区| 国产av码专区亚洲av| 成人亚洲精品一区在线观看 | 99久久精品国产国产毛片| 狂野欧美激情性xxxx在线观看| 婷婷六月久久综合丁香| 亚洲第一区二区三区不卡| 亚洲欧美日韩东京热| 秋霞伦理黄片| 久久久久性生活片| 校园人妻丝袜中文字幕| 国产成人福利小说| 五月玫瑰六月丁香| 男人和女人高潮做爰伦理| 观看美女的网站| 久久精品综合一区二区三区| 亚洲精品色激情综合| 国内揄拍国产精品人妻在线| 国产精品福利在线免费观看| freevideosex欧美| 国产精品综合久久久久久久免费| 三级男女做爰猛烈吃奶摸视频| 一个人看的www免费观看视频| 亚洲欧美日韩东京热| 精品少妇黑人巨大在线播放| 欧美日本视频| 亚洲av中文字字幕乱码综合| 黄色欧美视频在线观看| 十八禁国产超污无遮挡网站| 免费黄频网站在线观看国产| 精品亚洲乱码少妇综合久久| 精品亚洲乱码少妇综合久久| 精品人妻熟女av久视频| 一区二区三区乱码不卡18| 好男人视频免费观看在线| 在线观看免费高清a一片| 国产欧美另类精品又又久久亚洲欧美| 天美传媒精品一区二区| 搡老妇女老女人老熟妇| 国产精品女同一区二区软件| 国产高清国产精品国产三级 | 亚洲熟妇中文字幕五十中出| 一级毛片 在线播放| 又大又黄又爽视频免费| 亚洲成人av在线免费| 岛国毛片在线播放| 人人妻人人澡人人爽人人夜夜 | 看黄色毛片网站| 中文字幕av在线有码专区| 免费人成在线观看视频色| 免费黄频网站在线观看国产| 国产色爽女视频免费观看| 国产色爽女视频免费观看| 老师上课跳d突然被开到最大视频| 女的被弄到高潮叫床怎么办| 听说在线观看完整版免费高清| 97精品久久久久久久久久精品| 好男人视频免费观看在线| 国产精品精品国产色婷婷| 免费观看的影片在线观看| 国产av国产精品国产| 一级毛片我不卡| 日韩 亚洲 欧美在线| 欧美成人精品欧美一级黄| 亚洲av成人精品一区久久| 国产成人一区二区在线| 国产av码专区亚洲av| 午夜日本视频在线| 亚洲综合色惰| 欧美日韩视频高清一区二区三区二| 午夜亚洲福利在线播放| 一本久久精品| 欧美3d第一页| 久久这里有精品视频免费| 国产高潮美女av| av一本久久久久| 极品教师在线视频| 天堂网av新在线| 内射极品少妇av片p| 国产一级毛片七仙女欲春2| 久久久午夜欧美精品| 精品国产三级普通话版| 爱豆传媒免费全集在线观看| 99九九线精品视频在线观看视频| 精品久久久精品久久久| 三级国产精品片| 中文字幕亚洲精品专区| 亚洲av不卡在线观看| 国产高清三级在线| 国产伦在线观看视频一区| 永久网站在线| 精华霜和精华液先用哪个| 麻豆国产97在线/欧美| 国产精品久久视频播放| 亚洲成人久久爱视频| 老女人水多毛片| 伊人久久精品亚洲午夜| 麻豆久久精品国产亚洲av| 看非洲黑人一级黄片| 欧美另类一区| 欧美xxxx黑人xx丫x性爽| 99热这里只有是精品在线观看| 国产大屁股一区二区在线视频| 亚洲精华国产精华液的使用体验| 你懂的网址亚洲精品在线观看| 亚洲国产精品sss在线观看| 美女黄网站色视频| 最近中文字幕高清免费大全6| 亚洲欧美中文字幕日韩二区| 极品教师在线视频| 欧美潮喷喷水| 亚洲色图av天堂| 成人国产麻豆网| 最近视频中文字幕2019在线8| 能在线免费看毛片的网站| 久久人人爽人人片av| 中文字幕人妻熟人妻熟丝袜美| 欧美人与善性xxx| 我的老师免费观看完整版| 在线观看一区二区三区| 噜噜噜噜噜久久久久久91| 国产亚洲精品久久久com| 午夜福利成人在线免费观看| 欧美一级a爱片免费观看看| 人人妻人人澡人人爽人人夜夜 | 26uuu在线亚洲综合色| 一级毛片电影观看| 18禁在线无遮挡免费观看视频| 99热全是精品| videos熟女内射| 黄色配什么色好看| 神马国产精品三级电影在线观看| 国产欧美另类精品又又久久亚洲欧美| 91在线精品国自产拍蜜月| 韩国av在线不卡| 夜夜爽夜夜爽视频| 成人亚洲欧美一区二区av| 亚洲成人精品中文字幕电影| 日韩av免费高清视频| 精品久久久久久久末码| 国产综合懂色| 人妻一区二区av| 亚洲精品久久久久久婷婷小说| 一级片'在线观看视频| 欧美成人一区二区免费高清观看| 精品人妻熟女av久视频| 老师上课跳d突然被开到最大视频| 欧美+日韩+精品| 日韩人妻高清精品专区| 亚洲精品色激情综合| 亚洲国产精品成人综合色| 男人狂女人下面高潮的视频| 免费观看精品视频网站| 插逼视频在线观看| 亚洲国产av新网站| 欧美性感艳星| 国产人妻一区二区三区在| 成人综合一区亚洲| 你懂的网址亚洲精品在线观看| 欧美日韩精品成人综合77777| 一级二级三级毛片免费看| 久99久视频精品免费| 亚洲av中文字字幕乱码综合| 成人亚洲欧美一区二区av| 3wmmmm亚洲av在线观看| 听说在线观看完整版免费高清| 久久久久精品性色| 久久精品国产鲁丝片午夜精品| 国内精品一区二区在线观看| 亚洲精品,欧美精品| 国产国拍精品亚洲av在线观看| 久久午夜福利片| 国产单亲对白刺激| 色5月婷婷丁香| 午夜日本视频在线| 日韩欧美三级三区| 好男人视频免费观看在线| 国产精品久久久久久久久免| 久久久久免费精品人妻一区二区| 男女边摸边吃奶| 丰满人妻一区二区三区视频av| 免费观看无遮挡的男女| 亚洲精品日本国产第一区| 亚洲国产成人一精品久久久| 国产一级毛片在线| 乱码一卡2卡4卡精品| 国产一区二区三区综合在线观看 | 一区二区三区乱码不卡18| 免费看日本二区| 精品少妇黑人巨大在线播放| 日韩一本色道免费dvd| 2018国产大陆天天弄谢| 老司机影院成人| 国产大屁股一区二区在线视频| 一个人看的www免费观看视频| 久久久亚洲精品成人影院| 97超视频在线观看视频| 亚洲欧美日韩卡通动漫| 久久人人爽人人片av| 成人欧美大片| 国产免费福利视频在线观看| 久久精品久久久久久噜噜老黄| 男女边吃奶边做爰视频| 久久国内精品自在自线图片| 国产在线一区二区三区精| 亚洲,欧美,日韩| 日韩av不卡免费在线播放| 国产 亚洲一区二区三区 | 亚洲欧洲日产国产| 好男人视频免费观看在线| 亚洲精品国产av蜜桃| 黄色欧美视频在线观看| 亚洲av日韩在线播放| 国产成人精品婷婷| 麻豆成人午夜福利视频| 国产探花在线观看一区二区| 两个人的视频大全免费| 国产av不卡久久| 人妻少妇偷人精品九色| 中文在线观看免费www的网站| 国产熟女欧美一区二区| 九九久久精品国产亚洲av麻豆| 观看美女的网站| 亚洲av电影不卡..在线观看| 国产毛片a区久久久久| av免费观看日本| 少妇猛男粗大的猛烈进出视频 | 中文字幕制服av| 波野结衣二区三区在线| 热99在线观看视频| 青青草视频在线视频观看| 秋霞在线观看毛片| 久久久久久久久中文| 久久久久久九九精品二区国产| 国产精品国产三级国产av玫瑰| 国产精品久久视频播放| 日日啪夜夜爽| 亚洲成人av在线免费| 亚洲自偷自拍三级| 亚洲国产色片| 精品国内亚洲2022精品成人| 午夜福利网站1000一区二区三区| 一级毛片电影观看| 一级av片app| 天天一区二区日本电影三级| 久久这里只有精品中国| 成人无遮挡网站| 免费黄频网站在线观看国产| 人妻夜夜爽99麻豆av| 欧美性猛交╳xxx乱大交人| 91精品国产九色| 国产黄片美女视频| 七月丁香在线播放| av女优亚洲男人天堂| 国产激情偷乱视频一区二区| 男的添女的下面高潮视频| 国国产精品蜜臀av免费| 三级国产精品片| 色综合色国产| 亚洲综合色惰| 男人爽女人下面视频在线观看| 亚洲无线观看免费| 亚洲综合色惰| 日韩 亚洲 欧美在线| 老司机影院毛片| 男女国产视频网站| 久热久热在线精品观看| 一个人看的www免费观看视频| 日韩制服骚丝袜av| 18禁在线播放成人免费| 欧美极品一区二区三区四区| 国产老妇女一区| 久久久精品94久久精品| 婷婷色av中文字幕| 精品酒店卫生间| 国产男人的电影天堂91| 秋霞在线观看毛片| 亚洲经典国产精华液单| 又粗又硬又长又爽又黄的视频| 狂野欧美白嫩少妇大欣赏| 午夜视频国产福利| 国产 一区精品| 毛片一级片免费看久久久久| 久久久久久久久中文| 免费观看a级毛片全部| 久久99热这里只有精品18| 精品久久久久久电影网| 国产老妇伦熟女老妇高清| 美女被艹到高潮喷水动态| 亚洲欧美日韩无卡精品| 性色avwww在线观看| 免费看不卡的av| 91aial.com中文字幕在线观看| 精品久久久精品久久久| 日韩在线高清观看一区二区三区| 精品国产三级普通话版| 精品一区二区三区视频在线| 特级一级黄色大片| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲国产精品sss在线观看| 久久午夜福利片| kizo精华| 高清日韩中文字幕在线| 久久人人爽人人片av| av福利片在线观看| 国产亚洲最大av| 欧美成人午夜免费资源| ponron亚洲| 国产免费视频播放在线视频 | 久久99热这里只有精品18| 三级毛片av免费| 三级男女做爰猛烈吃奶摸视频| 日韩av免费高清视频| 午夜福利成人在线免费观看| 永久免费av网站大全| 免费观看无遮挡的男女| 少妇的逼水好多| 亚洲欧美成人综合另类久久久| 久99久视频精品免费| 免费无遮挡裸体视频| 久久久精品94久久精品| 亚洲精品乱码久久久v下载方式| 中文字幕av成人在线电影| 成人午夜高清在线视频| 99久久精品一区二区三区| 搡女人真爽免费视频火全软件| 韩国av在线不卡| 在线天堂最新版资源| 青青草视频在线视频观看| 夜夜爽夜夜爽视频| 少妇丰满av| 久久精品夜色国产| 美女脱内裤让男人舔精品视频| 欧美性感艳星| 精品国内亚洲2022精品成人| 国内揄拍国产精品人妻在线| 51国产日韩欧美| 精品不卡国产一区二区三区| 99九九线精品视频在线观看视频| 亚洲图色成人| 中文天堂在线官网| 国产黄色小视频在线观看| 建设人人有责人人尽责人人享有的 | 亚洲精品乱码久久久久久按摩| eeuss影院久久| 69人妻影院| 亚洲三级黄色毛片| 亚洲av国产av综合av卡| 69av精品久久久久久| 亚洲在久久综合| 免费不卡的大黄色大毛片视频在线观看 | 午夜久久久久精精品| 黄片无遮挡物在线观看| 大香蕉97超碰在线| 少妇高潮的动态图| 九色成人免费人妻av| 久久精品国产亚洲av天美| 亚洲av国产av综合av卡| 嫩草影院入口| 精品久久久久久久人妻蜜臀av| www.色视频.com| 晚上一个人看的免费电影| 国产永久视频网站| 国语对白做爰xxxⅹ性视频网站| 亚洲无线观看免费| 国产欧美另类精品又又久久亚洲欧美| 91精品国产九色| 亚洲av国产av综合av卡| 日本欧美国产在线视频| 久久久久国产网址| 2018国产大陆天天弄谢| 国产乱人偷精品视频| 日韩欧美精品v在线| 亚洲精品第二区| 午夜福利网站1000一区二区三区| 国产免费视频播放在线视频 | 中文精品一卡2卡3卡4更新| 别揉我奶头 嗯啊视频| 久久国产乱子免费精品| 午夜福利视频1000在线观看| 伦理电影大哥的女人| 乱人视频在线观看| 国产伦精品一区二区三区四那| 色吧在线观看| 嘟嘟电影网在线观看| 男女啪啪激烈高潮av片| 国内精品宾馆在线| 日本-黄色视频高清免费观看| 最近视频中文字幕2019在线8| 国产综合懂色| 婷婷色综合大香蕉| 嘟嘟电影网在线观看| 亚洲激情五月婷婷啪啪| videossex国产| 免费av不卡在线播放| 免费观看性生交大片5| 少妇的逼水好多| 亚洲18禁久久av| 真实男女啪啪啪动态图| 美女国产视频在线观看| 国产精品蜜桃在线观看| av免费观看日本| 女人被狂操c到高潮| 免费高清在线观看视频在线观看| av女优亚洲男人天堂| 高清午夜精品一区二区三区| 一区二区三区乱码不卡18| 国产伦一二天堂av在线观看| 日韩强制内射视频| 可以在线观看毛片的网站| 免费不卡的大黄色大毛片视频在线观看 | 一级爰片在线观看| 国产伦精品一区二区三区视频9| 中文字幕免费在线视频6| 欧美日韩精品成人综合77777| 成人鲁丝片一二三区免费| 国产av在哪里看| 韩国高清视频一区二区三区| 亚洲精品自拍成人| 久久久久精品性色| 亚洲精品aⅴ在线观看| 99久久精品国产国产毛片| 国产午夜福利久久久久久| 国产亚洲精品久久久com| 肉色欧美久久久久久久蜜桃 | 久99久视频精品免费| 高清午夜精品一区二区三区| 毛片女人毛片| 天堂中文最新版在线下载 | av免费观看日本| 亚洲一区高清亚洲精品| 深夜a级毛片| 男人狂女人下面高潮的视频| 日本午夜av视频| 成年av动漫网址| 97人妻精品一区二区三区麻豆| 国产美女午夜福利| 天天躁夜夜躁狠狠久久av| 韩国av在线不卡| 久久99精品国语久久久| 麻豆av噜噜一区二区三区| 2021天堂中文幕一二区在线观| 99热全是精品| 深夜a级毛片| 蜜桃亚洲精品一区二区三区| 亚洲精品成人久久久久久| 成人亚洲精品av一区二区| 亚洲国产av新网站| 丝袜喷水一区| 亚洲av福利一区| 80岁老熟妇乱子伦牲交| 亚洲精品国产成人久久av| 国产亚洲精品av在线|