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

    含隨機風(fēng)電出力及電網(wǎng)安全的電網(wǎng)調(diào)度策略研究

    2016-03-02 03:22:49羅亦梅
    四川電力技術(shù) 2016年1期
    關(guān)鍵詞:置信水平風(fēng)力發(fā)電

    羅亦梅,張 里,吳 琛,郝 瓏

    ?

    含隨機風(fēng)電出力及電網(wǎng)安全的電網(wǎng)調(diào)度策略研究

    羅亦梅1,張里2,吳琛3,郝瓏4

    0引言

    中國“十二五”規(guī)劃實施以來,作為一種綠色、清潔的可再生能源,風(fēng)電在中國北部和沿海地區(qū)迅速發(fā)展。大型風(fēng)電場的并網(wǎng)運行是目前中國利用風(fēng)能的主要模式,也是智能電網(wǎng)發(fā)展的重要組成部分[1]。大規(guī)模的風(fēng)電并網(wǎng)不僅緩解電力供應(yīng)緊張的形式,同時也減少了能源利用過程中大量溫室氣體的排放[2-3]。

    風(fēng)電與傳統(tǒng)的火電機組和水電機組相比,具有強烈的波動性和隨機性,其并網(wǎng)運行給電力系統(tǒng)調(diào)度帶來新的問題和挑戰(zhàn)[4-5]。針對風(fēng)電功率不確定性的經(jīng)濟調(diào)度問題,目前國內(nèi)外的一些學(xué)者進行了大量的研究。文獻[6]通過應(yīng)用機會約束規(guī)劃模擬風(fēng)電功率的隨機波動特性,建立了考慮機組組合的含風(fēng)電場電力系統(tǒng)經(jīng)濟調(diào)度模型;文獻[7]提出節(jié)能發(fā)電調(diào)度發(fā)電計劃的編制算法,首先根據(jù)優(yōu)先順序法確定火電機組開機方法,然后采用等微增率進行功率分配;文獻[8]提出在經(jīng)濟調(diào)度的目標函數(shù)中計及高估和低估風(fēng)電的經(jīng)濟代價,較合理地刻畫了風(fēng)電功率的概率分布特性,具有一定的工程實用性;文獻[9]為應(yīng)對風(fēng)電功率預(yù)測誤差給系統(tǒng)調(diào)度帶來的影響,在ED模型中同時考慮系統(tǒng)正、負旋轉(zhuǎn)備用容量約束;文獻[10]應(yīng)用模糊集理論建立了含風(fēng)電場經(jīng)濟調(diào)度的模糊模型;文獻[11]構(gòu)建含風(fēng)電場環(huán)境經(jīng)濟調(diào)度的隨機優(yōu)化模型,將燃煤機組的發(fā)電成本最小和污染氣體排放量最小同時最為目標函數(shù)進行優(yōu)化。

    上述發(fā)電調(diào)度模型中并未考慮隨機風(fēng)電出力,以及輸電元件傳輸容量的約束。在傳統(tǒng)經(jīng)濟調(diào)度模型中考慮了網(wǎng)絡(luò)安全約束,并將風(fēng)速的概率模型和風(fēng)機出力相結(jié)合,通過分段函數(shù)表達風(fēng)電場出力的概率分布函數(shù),利用風(fēng)電場出力的分布函數(shù)對模型進行確定性轉(zhuǎn)化,使得仿真結(jié)果更接近實際情況。最后利用量子行為粒子群算法對模型進行求解,同時對調(diào)度模型進行敏感性分析。

    1風(fēng)電出力的概率模型

    通過長期實測數(shù)據(jù)統(tǒng)計和擬合分析,普遍認為二參數(shù)的威布爾分布能夠較好描述風(fēng)速v的隨機分布規(guī)律[12]。分布函數(shù)為

    (1)

    式中:k為形狀參數(shù);c為威布爾分布的尺度參數(shù)。根據(jù)其分布函數(shù)可得威布爾分布的概率密度函數(shù)為

    (2)

    由于影響風(fēng)能輸入和功率輸出之間關(guān)系的因素很多,所以無法在建模中一一分析。研究人員使用通用、簡化的模型[13]來描述風(fēng)速和輸出功率的關(guān)系。

    (3)

    式中:W為風(fēng)電實際出力;V為實際風(fēng)速;wr為風(fēng)電場的額定容量;vin、vr、vout分別為風(fēng)電機的切入風(fēng)速、額定風(fēng)速和切除風(fēng)速。

    在經(jīng)濟調(diào)度中將采用以上模型,根據(jù)隨機變量函數(shù)的分布[14]可知,vin≤V

    (4)

    式中,h=(vr/vin)-1。在計算分布函數(shù)時需要采用分段線性化。

    p(W=0)=p(V

    (5)

    式中,p為概率。

    P(W=wr)=P(vr≤V

    (6)

    為求取連續(xù)部分的概率分布,對式(4)積分為

    (7)

    此外,

    P(W>wr)=0

    (8)

    綜合式(5)~式(8),得到W的概率分布函數(shù)為

    FW(w)=P(W≤w)=

    (9)

    2電力系統(tǒng)安全經(jīng)濟調(diào)度

    2.1 目標函數(shù)

    發(fā)電成本為

    (10)

    式中:N為系統(tǒng)內(nèi)燃煤機組數(shù);Pi為燃煤機組i的有功出力;FC為燃煤機組的總成本。

    2.2 約束條件

    1)功率平衡約束

    (11)

    式中:PD為系統(tǒng)負荷需求;Ω(PW)為風(fēng)電場出力PW的函數(shù)。由于PW為隨機變量,因此可將式(11)以概率的形式描述為

    (12)

    式中,η為滿足負荷需求的置信水平。

    2)機組出力約束

    燃煤機組的出力約束為

    Pimin≤Pi≤Pimax

    (13)

    風(fēng)電場的出力約束為

    0≤PW≤PWmax

    (14)

    3)網(wǎng)絡(luò)安全約束

    -PL≤BdiagLB-1[Pg-PD]≤PL

    (15)

    (16)

    式中:B為系統(tǒng)的導(dǎo)納矩陣;xi為支路i的電抗;Nl為系統(tǒng)支路條數(shù);L為系統(tǒng)支路節(jié)點的連接矩陣;PL為支路功率約束向量。

    2.3 模型的確定性轉(zhuǎn)化

    求解機會約束規(guī)劃問題的一種常用的方法是將其轉(zhuǎn)化為確定性優(yōu)化模型。式(12)可轉(zhuǎn)化為

    (17)

    整理后:

    (18)

    3求解算法

    3.1 量子行為粒子群算法

    含隨機風(fēng)電出力及電網(wǎng)安全的電網(wǎng)經(jīng)濟調(diào)度模型是一個多變量、非線性、混合整數(shù)規(guī)劃問題。在Matlab平臺上利用量子行為粒子群算法[16-17]求解計及網(wǎng)絡(luò)安全約束時的電力系統(tǒng)經(jīng)濟調(diào)度模型,該算法具有全局收斂性、控制參數(shù)更少,收斂速度更快、尋優(yōu)能力強等優(yōu)點。

    在粒子群中引入一個平均最好位置C(t)來計算粒子下一步迭代的變量,定義為所有粒子個體最好位置的平均,即

    (19)

    式中:M為粒子的個數(shù);Pi為粒子i的局部最好位置。于是粒子的迭代方程變?yōu)?/p>

    x(t+1)=p±λ×|Cj(t)-x(t)|×ln[1/u]

    (20)

    式中:u=rand(0,1);λ稱為收縮-擴張系數(shù),調(diào)節(jié)它可以控制算法的收斂速度。這里λ取值為

    λ=(1-0.5)×(Maxiter-t)/Maxiter+0.5

    (21)

    式中:Maxiter為最大迭代次數(shù);t為當前迭代次數(shù)。式(19)~式(21)稱為具有量子行為粒子群算法。

    3.2 算法流程

    根據(jù)數(shù)學(xué)模型,采用基于隨機模擬的量子行為粒子群算法流程如下:

    1) 讀取系統(tǒng)數(shù)據(jù),如輸入負荷、火電機組參數(shù)、線路參數(shù)和風(fēng)電出力的預(yù)測值等。再次是輸入量子行為粒子群算法參數(shù),如最大迭代次數(shù)、粒子群大小。

    2) 種群初始化。隨機生成各時段各火電機組有功功率和備用容量,組成一個個體,根據(jù)式(11)~式(16)檢驗個體的可行性;若個體不可行,則重新生成個體,直至初始群體生成完畢。

    3) 首先根據(jù)式(19)計算粒子群的平均最好位置,計算粒子當前位置的適應(yīng)值。

    4) 粒子位置的更新。根據(jù)式(20)對每個粒子的位置進行更新,同時驗證是否越限,如果越限則需要重新更新粒子。然后計算每個粒子的適應(yīng)值,如果優(yōu)于當前粒子個體極值,則更新個體極值。如果所有粒子的個體極值有優(yōu)于當前的全局極值,則更新全局最優(yōu)。

    5) 判斷是否達到迭代次數(shù)。如果未達到則返回步驟3)、步驟4);否則輸出最好的粒子作為最優(yōu)解。

    4算例分析

    4.1 算例描述

    利用文獻[15]中含有5臺燃煤機組和1個風(fēng)電場系統(tǒng)進行仿真計算。系統(tǒng)接線圖見圖1,燃煤機組的發(fā)電成本系數(shù)和機組出力限制見表1,線路參數(shù)和負荷分布參數(shù)見表2及表3。系統(tǒng)負荷為1 100 MW。網(wǎng)絡(luò)約束在這里只考慮線路的傳輸容量約束。在編程求解模型時,粒子群個數(shù)設(shè)置為40個,迭代運行次數(shù)為500。

    圖1 仿真系統(tǒng)接線

    機組a/(USD·MW-2)b/(USD·MW-1)c/USDPmin/MWPmax/MWG10.00506.0100.0050200G20.00258.0300.00100400G30.002010.0500.00100600G40.007014.560.001050G50.038026.880.001240

    表2 線路參數(shù)

    表3 節(jié)點負荷分布

    表4 風(fēng)電場參數(shù)

    表5 不同置信水平η下的機組出力

    風(fēng)電場參數(shù)見表4,威布爾分布的參數(shù)k=2,c=12 m/s。

    4.2 算例驗證

    4.2.1在是否計及網(wǎng)絡(luò)安全約束下,不同置信水平對機組出力的影響

    取形狀參數(shù)k=2,尺度參數(shù)c=12 m/s。分別對置信水平η為0.8、0.7、0.6和0.5下進行計算,表5給出了是否計及網(wǎng)絡(luò)安全約束下的機組出力。

    對計算結(jié)果進行分析可知,隨著滿足負荷需求置信水平η的下降,系統(tǒng)接入的風(fēng)電場出力增加,燃煤機組出力下降,同時燃煤機組發(fā)電成本降低。在不考慮網(wǎng)絡(luò)約束情況下,最小發(fā)電成本由11 067 $/h降至9 982 $/h;而考慮網(wǎng)絡(luò)安全約束情況下,最小發(fā)電成本由11 067 $/h降至10 205 $/h。

    同時針對相同的置信水平,考慮網(wǎng)絡(luò)安全約束比不考慮網(wǎng)絡(luò)安全約束所需燃煤機組成本高。這說明隨著置信水平的降低,系統(tǒng)接入的風(fēng)電場出力增加,并同時伴隨著部分線路潮流越線的情況,考慮網(wǎng)絡(luò)安全約束后成本增加的本質(zhì)是由發(fā)電成本高的機組代替發(fā)電成本低的機組出力,從而滿足電網(wǎng)處于安全運行狀態(tài)。

    4.2.2在計及網(wǎng)絡(luò)安全約束下,不同負荷對機組出力的影響

    取形狀參數(shù)k=2,尺度參數(shù)c=12 m/s。置信水平η為0.6,取負荷分別為1 000 MW、1 100 MW和1 200 MW進行優(yōu)化。表6給出了不同負荷下的調(diào)度結(jié)果。

    表6 不同負荷下的機組出力

    對計算結(jié)果進行分析,在相同的置信水平下,風(fēng)電場出力為不變量。隨著負荷的增加,未達到最大出力的機組會增加出力,同時所需成本也會相應(yīng)增加。

    4.2.3在計及網(wǎng)絡(luò)安全約束下,不同尺度參數(shù)對機組出力的影響

    取形狀參數(shù)k=2,置信水平η為0.6,分別取尺度參數(shù)為12 m/s、15 m/s、18 m/s和21 m/s進行優(yōu)化。表7給出了發(fā)電成本最少的調(diào)度結(jié)果。

    表7 不同尺度參數(shù)下的機組出力

    對計算結(jié)果進行分析可知,尺度參數(shù)從12 m/s增加到21 m/s,風(fēng)電場出力從90.3 MW增加到145.6 MW,從而降低了燃煤機組的成本。結(jié)果顯示,尺度參數(shù)反映了風(fēng)電場的平均風(fēng)速,因此隨著尺度參數(shù)的增加風(fēng)電場的出力增加,但是增幅卻在減少,同時伴隨著發(fā)電成本在逐漸降低。

    4.2.4在計及網(wǎng)絡(luò)安全約束下,不同形狀參數(shù)對機組出力的影響

    取置信水平η=0.6,尺度參數(shù)c=12 m/s。在計及網(wǎng)絡(luò)安全約束下,分別計算形狀參數(shù)k=1.8、2.0、2.2和2.4下的燃煤機組和風(fēng)電場調(diào)度結(jié)果。

    表8 不同形狀參數(shù)的機組出力

    比較計算結(jié)果可知,隨著形狀參數(shù)增加,系統(tǒng)風(fēng)電場出力從82.4 MW增加到101.5 MW,同時燃煤機組總出力減小,最小發(fā)電成本由10 520 $/h下降至10 370 $/h。

    圖2 不同形狀參數(shù)k下的風(fēng)速概率密度函數(shù)

    從圖2可知,隨著k值增加,風(fēng)速概率密度函數(shù)曲線逐漸變陡,且密度函數(shù)的峰值逐漸增加,同時出現(xiàn)概率峰值對應(yīng)的風(fēng)速值也逐漸增大,但增加幅度變小。這表明,隨著k值增加,風(fēng)速分布更接近于風(fēng)速峰值附近,概率峰值對應(yīng)的風(fēng)速值也增加,風(fēng)電場的出力增加。

    5結(jié)論

    1) 建立了含隨機風(fēng)電出力和電網(wǎng)安全的電網(wǎng)經(jīng)濟調(diào)度模型,并利用風(fēng)電出力的分布函數(shù)將風(fēng)電不確定性轉(zhuǎn)化為確定性模型;

    2) 在考慮網(wǎng)絡(luò)安全約束時,為了使電網(wǎng)安全運行,部分機組出力結(jié)果的變化導(dǎo)致系統(tǒng)成本增加;

    3) 分析了在考慮網(wǎng)絡(luò)安全約束下,對調(diào)度模型進行敏感性分析;

    4) 仿真結(jié)果表明了所提調(diào)度方法的合理性和有效性,為含隨機風(fēng)電出力和電網(wǎng)安全的電網(wǎng)經(jīng)濟調(diào)度提供參考。

    參考文獻

    [1] 王淑娟,郭元. 我國大規(guī)模發(fā)展風(fēng)電的機制設(shè)計[J].

    太陽能學(xué)報,2001,22 (1): 59-63.

    [2]Denny E. Wind Generation, Power System Operation,

    and Emissions Reduction [J]. IEEE Transaction on Power Systems,2006,21(1): 341-347.

    [3]王睿,張粒子,張麗娟,等. 含風(fēng)電場系統(tǒng)風(fēng)電與抽水蓄能匹配容量研究[J]. 太陽能學(xué)報,2012,33(6): 1039-1043.

    [4]J. Wang, A. Botterud, R. Bessa. Wind Power Forecasting Uncertainty and Unit Commitment [J]. Applied Energy,2011,88(11): 4014-4023.

    [5]K. Methaprayoon, W. J. Lee. An Integration of ANN Wind Power Estimation into UC Considering the Forecasting Uncertainty [C]. Industrial and Commercial Power Systems Technical Conference,2005.

    [6]孫元章,吳俊,李國杰,等. 基于風(fēng)速測量和隨機規(guī)劃的含風(fēng)電場電力系統(tǒng)動態(tài)經(jīng)濟調(diào)度[J].中國電機工程學(xué)報,2009,92(4): 41-47.

    [7]施建華,譚素梅.節(jié)能發(fā)電調(diào)度發(fā)電計劃編制算法[J].電力系統(tǒng)自動化,2008,32(24): 48-51.

    [8]董曉天,嚴正,馮冬涵,等. 計及風(fēng)電出力懲罰成本的電力系統(tǒng)經(jīng)濟調(diào)度[J]. 電網(wǎng)技術(shù),2012,36(8): 76-80.

    [9]周瑋,彭昱,孫輝,等. 含風(fēng)電場的電力系統(tǒng)動態(tài)經(jīng)濟調(diào)度[J]. 中國電機工程學(xué)報,2009,29(25): 13-18.

    [10]陳海焱,陳金富,段獻忠. 含風(fēng)電場電力系統(tǒng)經(jīng)濟調(diào)度的模糊建模及優(yōu)化算法[J]. 電力系統(tǒng)自動化,2006,30(2): 22-26.

    [11]邱威,張建華,劉念. 含大型風(fēng)電場的環(huán)境經(jīng)濟調(diào)度模型和解法[J]. 中國電機工程學(xué)報,2011,31(19): 8-16.

    [12]王成山,鄭海峰,謝瑩華,等. 計及分布式發(fā)電的配電系統(tǒng)隨機潮流計算[J]. 電力系統(tǒng)自動化,2005,29(24): 39-44.

    [13]雷亞洲,王偉勝,印永華,等. 風(fēng)電對電力系統(tǒng)運行的價值分析[J], 電網(wǎng)技術(shù),2002,26(5):10-14.

    [14]Stiebler M.風(fēng)力發(fā)電系統(tǒng)[M]. 北京:機械工業(yè)出版社,2010: 10-17.

    [15]Wood Allen J, Wollenberg Bruce F. Power Generation, Operation and Control[M]. New York: John Wiley & Sons Inc, 1996.

    [16]S.N. Omkar, G.Narayana Naik. Quantum Behaved P-

    article Swarm Optimization for Multi-objective Design Optimization of Composite Structures [J]. Expert System with Application, 2009, 36(8): 11312-11322.

    [17]Yejie Cai, Jun Sun. Optimizing the Codon Usage of Sy-

    nthetic Gene with QPSO Algorithm [J]. Journal of Theoretical Biology, 2008, 254 (1): 123-127.

    張里(1984),博士研究生,從事風(fēng)力發(fā)電方面研究。

    (1.成都理工大學(xué)核技術(shù)與自動化工程學(xué)院, 四川 成都610059;

    2. 國網(wǎng)四川省電力公司技能培訓(xùn)中心,四川 成都610072;

    3. 國網(wǎng)重慶市電力公司檢修分公司,重慶400039; 4. 國網(wǎng)保定供電公司,河北 保定071000)

    摘要:隨著風(fēng)電基地在中國北部和沿海地區(qū)大力建設(shè),風(fēng)電在電網(wǎng)中的比例不斷提高,這對電力系統(tǒng)的經(jīng)濟調(diào)度運行提出新的考驗。特別是風(fēng)電的間歇性和波動性使得經(jīng)濟調(diào)度具有更多不確定性。將風(fēng)速的概率模型和風(fēng)機出力相結(jié)合,利用風(fēng)電場出力的概率分布函數(shù)將風(fēng)電場出力的隨機模型轉(zhuǎn)化為確定性模型,構(gòu)建以機組的發(fā)電成本最小為目標函數(shù),同時計及電網(wǎng)網(wǎng)絡(luò)安全約束。在優(yōu)化算法上,利用量子行為粒子群算法對模型進行求解,最后對調(diào)度模型進行敏感性分析。仿真結(jié)果表明,所提模型和方法是可行、有效的,具有一定的實用價值。

    關(guān)鍵詞:風(fēng)力發(fā)電;經(jīng)濟調(diào)度;網(wǎng)絡(luò)安全約束;置信水平 ;量子行為粒子群算法

    Abstract:Wind generation has attracted great attention with the increase of the capacity of large wind farms connected to power grid, thus it brings new challenges to economic dispatch of traditional power system. Especially, more uncertainties in economic dispatch are produced by the random variation of wind power output. The wind velocity probability model is combined with wind generator output, and the probability distribution function of wind power output is expressed by piecewise function. According to the cumulative distribution function of wind power, the stochastic optimization problem is transformed to a deterministic one. Minimizing the fuel cost of thermal generators is taken as objective functions, at the same time, the network security constraints are considered. Finally, the effect of different load, shape parameter and scale factor of the scheduling results are considered too. The effectiveness of the proposed method is verified by the results of calculation examples.

    Key words:wind power; economic dispatch; network security constraint; confidence level; quantum behaved particle swarm optimization

    (收稿日期:2015-11-02)

    作者簡介:

    基金項目:國家自然科學(xué)基金項目(50977059);國家電網(wǎng)公司科技項目——大規(guī)模風(fēng)電并網(wǎng)調(diào)度運行支撐關(guān)鍵技術(shù)研究與應(yīng)用

    中圖分類號:TM73

    文獻標志碼:A

    文章編號:1003-6954(2016)01-0006-06

    猜你喜歡
    置信水平風(fēng)力發(fā)電
    產(chǎn)品控制與市場風(fēng)險之間的相互作用研究
    單因子方差分析法在卷煙均勻性檢驗中的研究與應(yīng)用
    用VaR方法分析中國A股市場的風(fēng)險
    風(fēng)力發(fā)電機組的電氣控制
    東汽FD77雙饋式變速恒頻風(fēng)力發(fā)電機的控制策略
    電池儲能技術(shù)在風(fēng)力發(fā)電系統(tǒng)中的應(yīng)用
    中國市場(2016年32期)2016-12-06 11:06:23
    風(fēng)力發(fā)電發(fā)展現(xiàn)狀及前景淺析
    價值工程(2016年31期)2016-12-03 00:07:43
    垂直軸風(fēng)力發(fā)電機葉片建模與研究農(nóng)用
    引信可靠性考核的系統(tǒng)性錯誤及統(tǒng)計檢驗方法
    考慮風(fēng)電置信水平的機組組合優(yōu)化方法
    久久这里只有精品19| av国产免费在线观看| 欧美乱码精品一区二区三区| 国产aⅴ精品一区二区三区波| 首页视频小说图片口味搜索| 久9热在线精品视频| 日本精品一区二区三区蜜桃| 熟妇人妻久久中文字幕3abv| 两个人看的免费小视频| 亚洲欧美日韩高清专用| 国产三级黄色录像| 国产私拍福利视频在线观看| 亚洲 国产 在线| 久久国产精品影院| 国产精品亚洲美女久久久| 亚洲在线自拍视频| 色综合欧美亚洲国产小说| 亚洲欧美激情综合另类| 人成视频在线观看免费观看| 国产一区二区在线av高清观看| 天堂动漫精品| 日韩 欧美 亚洲 中文字幕| 免费在线观看影片大全网站| 亚洲成av人片在线播放无| 色尼玛亚洲综合影院| 欧美黄色淫秽网站| 国产麻豆成人av免费视频| 亚洲七黄色美女视频| 亚洲av成人av| 国产一区在线观看成人免费| 女人高潮潮喷娇喘18禁视频| 欧美日韩黄片免| 国产亚洲精品综合一区在线观看 | 成人欧美大片| 国产高清videossex| 国产69精品久久久久777片 | 国产精品精品国产色婷婷| 一区福利在线观看| 国产精品av视频在线免费观看| 亚洲乱码一区二区免费版| 亚洲成人免费电影在线观看| 国产精品亚洲一级av第二区| 一二三四在线观看免费中文在| 国产精品乱码一区二三区的特点| 免费在线观看视频国产中文字幕亚洲| 精品一区二区三区av网在线观看| 国产精品久久电影中文字幕| 成人三级做爰电影| 久久久精品国产亚洲av高清涩受| 免费搜索国产男女视频| 欧美黑人巨大hd| 夜夜爽天天搞| 国产主播在线观看一区二区| 午夜激情av网站| 国产精品1区2区在线观看.| 久久久久国产精品人妻aⅴ院| 免费在线观看日本一区| 成人一区二区视频在线观看| 搡老熟女国产l中国老女人| 日本一本二区三区精品| 级片在线观看| 婷婷丁香在线五月| 三级毛片av免费| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲精品中文字幕一二三四区| 国产精品av久久久久免费| 久久久久久久午夜电影| 丰满人妻一区二区三区视频av | 舔av片在线| 国产av一区二区精品久久| 国产精品av视频在线免费观看| 国产精品一及| 色综合亚洲欧美另类图片| 亚洲,欧美精品.| 99热只有精品国产| 精品国产美女av久久久久小说| 午夜免费观看网址| 老汉色av国产亚洲站长工具| 亚洲精品粉嫩美女一区| 亚洲av日韩精品久久久久久密| 国产精品一及| 免费电影在线观看免费观看| 国产真人三级小视频在线观看| 女人高潮潮喷娇喘18禁视频| 亚洲真实伦在线观看| 成人av一区二区三区在线看| 亚洲一区二区三区色噜噜| 亚洲精品中文字幕在线视频| 免费无遮挡裸体视频| 欧美性长视频在线观看| 婷婷亚洲欧美| 国产精品1区2区在线观看.| 午夜激情av网站| 动漫黄色视频在线观看| 少妇粗大呻吟视频| 嫩草影视91久久| 亚洲 国产 在线| 琪琪午夜伦伦电影理论片6080| 国产99白浆流出| 身体一侧抽搐| 全区人妻精品视频| 欧美日韩黄片免| 婷婷精品国产亚洲av在线| 国产主播在线观看一区二区| 日韩欧美三级三区| 五月玫瑰六月丁香| 极品教师在线免费播放| 性色av乱码一区二区三区2| 欧美成人一区二区免费高清观看 | 国产一区二区激情短视频| 床上黄色一级片| 亚洲国产欧洲综合997久久,| 丝袜美腿诱惑在线| 男人的好看免费观看在线视频 | 国产成人精品久久二区二区免费| 看免费av毛片| 亚洲中文日韩欧美视频| 香蕉久久夜色| 欧美乱色亚洲激情| 国产成+人综合+亚洲专区| 在线十欧美十亚洲十日本专区| 免费无遮挡裸体视频| 国产一区二区三区视频了| av福利片在线观看| 91九色精品人成在线观看| 香蕉久久夜色| 国产亚洲精品av在线| 搞女人的毛片| 亚洲欧美日韩东京热| 超碰成人久久| 国产精品98久久久久久宅男小说| 国产精品av视频在线免费观看| 午夜福利18| 此物有八面人人有两片| 国产高清视频在线播放一区| 欧美三级亚洲精品| 麻豆久久精品国产亚洲av| av天堂在线播放| 亚洲精品在线美女| 亚洲第一电影网av| 午夜激情福利司机影院| 热99re8久久精品国产| 久久久久国产一级毛片高清牌| 亚洲一区高清亚洲精品| 午夜福利成人在线免费观看| 精品日产1卡2卡| 国产人伦9x9x在线观看| 性欧美人与动物交配| 成年女人毛片免费观看观看9| 成人三级做爰电影| 美女扒开内裤让男人捅视频| 色av中文字幕| 九色国产91popny在线| 精品久久久久久久久久久久久| 亚洲男人的天堂狠狠| 欧美另类亚洲清纯唯美| 国产高清videossex| 人妻丰满熟妇av一区二区三区| 啦啦啦韩国在线观看视频| 色综合婷婷激情| 亚洲国产精品合色在线| 午夜激情av网站| 日本熟妇午夜| 中文亚洲av片在线观看爽| 日韩免费av在线播放| 久久人人精品亚洲av| 国产三级在线视频| 国产不卡一卡二| 国产av不卡久久| 国产一区二区三区在线臀色熟女| 亚洲精华国产精华精| 后天国语完整版免费观看| 别揉我奶头~嗯~啊~动态视频| 免费看日本二区| 国产精品 欧美亚洲| 成人18禁高潮啪啪吃奶动态图| 国产激情偷乱视频一区二区| 欧美性长视频在线观看| 久久婷婷成人综合色麻豆| 国产成人系列免费观看| 麻豆国产97在线/欧美 | 国产精品久久视频播放| 午夜精品在线福利| 日本一区二区免费在线视频| 精品午夜福利视频在线观看一区| 99热这里只有是精品50| 国产伦一二天堂av在线观看| 久久久久国内视频| 中文字幕熟女人妻在线| 麻豆av在线久日| 国产日本99.免费观看| svipshipincom国产片| 91老司机精品| 香蕉久久夜色| 老汉色∧v一级毛片| 国产男靠女视频免费网站| 18禁国产床啪视频网站| av在线天堂中文字幕| 久久久久国内视频| 欧美3d第一页| 人人妻人人看人人澡| 日本成人三级电影网站| 亚洲真实伦在线观看| 长腿黑丝高跟| 欧美性猛交╳xxx乱大交人| 青草久久国产| 色av中文字幕| 看黄色毛片网站| 精品不卡国产一区二区三区| 每晚都被弄得嗷嗷叫到高潮| 99久久久亚洲精品蜜臀av| 国产三级在线视频| 波多野结衣高清作品| 91国产中文字幕| 日本三级黄在线观看| 国产亚洲精品一区二区www| 午夜视频精品福利| 欧美黑人精品巨大| 露出奶头的视频| 久久性视频一级片| 成年免费大片在线观看| 精品久久久久久久久久免费视频| 91在线观看av| 欧美一区二区精品小视频在线| 草草在线视频免费看| 亚洲无线在线观看| 欧美成狂野欧美在线观看| 熟女少妇亚洲综合色aaa.| netflix在线观看网站| 三级男女做爰猛烈吃奶摸视频| 最好的美女福利视频网| 久久久久久国产a免费观看| 桃红色精品国产亚洲av| 国产区一区二久久| 精品久久久久久久人妻蜜臀av| 精品少妇一区二区三区视频日本电影| 午夜福利在线在线| 长腿黑丝高跟| 桃红色精品国产亚洲av| 精品久久蜜臀av无| 亚洲在线自拍视频| 亚洲人成网站在线播放欧美日韩| 国产高清视频在线播放一区| 久久99热这里只有精品18| 欧美一区二区国产精品久久精品 | www国产在线视频色| 在线观看舔阴道视频| 国产精品国产高清国产av| 久久人人精品亚洲av| 大型av网站在线播放| 美女扒开内裤让男人捅视频| 99久久99久久久精品蜜桃| 欧美三级亚洲精品| 日韩欧美三级三区| 在线观看免费日韩欧美大片| 欧美另类亚洲清纯唯美| netflix在线观看网站| 国产一区二区三区视频了| 久久99热这里只有精品18| 亚洲av熟女| 国产97色在线日韩免费| 成人精品一区二区免费| 久久久久久九九精品二区国产 | 欧美另类亚洲清纯唯美| 天天一区二区日本电影三级| 午夜免费成人在线视频| 欧美黑人欧美精品刺激| 久久亚洲真实| 欧美成人免费av一区二区三区| 国产午夜精品论理片| 国产亚洲精品第一综合不卡| 变态另类丝袜制服| 国产av又大| 亚洲av成人精品一区久久| 亚洲成人免费电影在线观看| 99re在线观看精品视频| 亚洲七黄色美女视频| 国产成人欧美在线观看| 一本综合久久免费| 亚洲真实伦在线观看| www.999成人在线观看| 久久午夜亚洲精品久久| 最新美女视频免费是黄的| 曰老女人黄片| 亚洲精品一区av在线观看| 欧美日韩中文字幕国产精品一区二区三区| 久久久久国内视频| 亚洲av成人精品一区久久| 午夜两性在线视频| 黄色视频不卡| 青草久久国产| 中文字幕最新亚洲高清| 国产午夜精品论理片| 精品人妻1区二区| 欧美性猛交黑人性爽| 国产成人一区二区三区免费视频网站| 国产激情欧美一区二区| 欧美日韩黄片免| 亚洲18禁久久av| 日韩高清综合在线| 两个人看的免费小视频| 欧美性长视频在线观看| 美女免费视频网站| 国产成人精品久久二区二区91| 性欧美人与动物交配| 91在线观看av| 成人国产一区最新在线观看| 久久香蕉国产精品| 这个男人来自地球电影免费观看| 人妻久久中文字幕网| 国产伦在线观看视频一区| av天堂在线播放| 国产精品九九99| 一级作爱视频免费观看| 成人一区二区视频在线观看| 精品国产超薄肉色丝袜足j| 18禁国产床啪视频网站| 丝袜人妻中文字幕| 国产片内射在线| 小说图片视频综合网站| 国产精品久久久久久精品电影| 中亚洲国语对白在线视频| 亚洲成a人片在线一区二区| 日日夜夜操网爽| 亚洲 国产 在线| 床上黄色一级片| 亚洲五月婷婷丁香| 最新美女视频免费是黄的| 亚洲va日本ⅴa欧美va伊人久久| 黄色 视频免费看| 久久人人精品亚洲av| 制服丝袜大香蕉在线| 午夜老司机福利片| 欧美大码av| 欧美高清成人免费视频www| 精品国产超薄肉色丝袜足j| 国产三级中文精品| 亚洲国产欧美人成| 很黄的视频免费| 99久久99久久久精品蜜桃| 国产免费av片在线观看野外av| 一级作爱视频免费观看| 无限看片的www在线观看| 亚洲男人天堂网一区| 国产区一区二久久| 小说图片视频综合网站| 99精品欧美一区二区三区四区| 精品无人区乱码1区二区| av福利片在线观看| 久久久久久久午夜电影| 欧美国产日韩亚洲一区| ponron亚洲| 99国产综合亚洲精品| √禁漫天堂资源中文www| 亚洲精品在线美女| 亚洲一码二码三码区别大吗| 国内精品久久久久精免费| 黄色女人牲交| 中出人妻视频一区二区| 欧美+亚洲+日韩+国产| 国产一区二区三区在线臀色熟女| 亚洲第一欧美日韩一区二区三区| 久久精品91蜜桃| 天堂√8在线中文| 欧美一级毛片孕妇| 精品国产亚洲在线| 最近视频中文字幕2019在线8| 91麻豆av在线| 大型黄色视频在线免费观看| 欧美在线一区亚洲| 50天的宝宝边吃奶边哭怎么回事| 此物有八面人人有两片| 久久亚洲真实| 欧美日韩国产亚洲二区| 亚洲国产欧洲综合997久久,| 香蕉国产在线看| 亚洲av片天天在线观看| 村上凉子中文字幕在线| 日韩av在线大香蕉| 啦啦啦观看免费观看视频高清| 亚洲真实伦在线观看| 久久人妻福利社区极品人妻图片| 久久人妻av系列| 国产成人系列免费观看| 欧美人与性动交α欧美精品济南到| 一二三四在线观看免费中文在| 欧美性猛交╳xxx乱大交人| 精品国产乱子伦一区二区三区| 日韩欧美国产一区二区入口| 久久久久九九精品影院| 天天一区二区日本电影三级| 琪琪午夜伦伦电影理论片6080| 久久午夜综合久久蜜桃| 中文字幕人成人乱码亚洲影| 91国产中文字幕| 人人妻人人看人人澡| 人人妻人人澡欧美一区二区| 国产黄色小视频在线观看| 两个人看的免费小视频| 神马国产精品三级电影在线观看 | 99在线视频只有这里精品首页| 午夜日韩欧美国产| 一级a爱片免费观看的视频| 免费搜索国产男女视频| 男女那种视频在线观看| 每晚都被弄得嗷嗷叫到高潮| 亚洲性夜色夜夜综合| 蜜桃久久精品国产亚洲av| 亚洲 欧美 日韩 在线 免费| 老汉色av国产亚洲站长工具| 18禁黄网站禁片午夜丰满| netflix在线观看网站| 国产在线观看jvid| 色播亚洲综合网| 欧美大码av| 日韩欧美国产一区二区入口| 中国美女看黄片| 91九色精品人成在线观看| 亚洲,欧美精品.| 性欧美人与动物交配| 999久久久精品免费观看国产| 极品教师在线免费播放| 最近最新中文字幕大全电影3| 国产亚洲精品第一综合不卡| 亚洲成av人片免费观看| 老司机深夜福利视频在线观看| 国产成年人精品一区二区| 国产v大片淫在线免费观看| 亚洲精品中文字幕在线视频| 亚洲九九香蕉| 色av中文字幕| 每晚都被弄得嗷嗷叫到高潮| 久久精品国产亚洲av高清一级| 高潮久久久久久久久久久不卡| 无遮挡黄片免费观看| 狂野欧美白嫩少妇大欣赏| 悠悠久久av| www.www免费av| 曰老女人黄片| 男女之事视频高清在线观看| 久久久精品大字幕| 国产主播在线观看一区二区| 免费在线观看成人毛片| 777久久人妻少妇嫩草av网站| 久久久国产精品麻豆| 久久午夜综合久久蜜桃| 日韩欧美 国产精品| 国内毛片毛片毛片毛片毛片| www.熟女人妻精品国产| 久久久久久人人人人人| 久久精品国产亚洲av高清一级| 制服人妻中文乱码| 女人高潮潮喷娇喘18禁视频| 欧美最黄视频在线播放免费| 国产私拍福利视频在线观看| 久久香蕉激情| 激情在线观看视频在线高清| 国产午夜精品论理片| 精品欧美一区二区三区在线| 国产精品乱码一区二三区的特点| 欧美黑人精品巨大| 亚洲一卡2卡3卡4卡5卡精品中文| 少妇熟女aⅴ在线视频| 国产成年人精品一区二区| 国产亚洲欧美98| 亚洲国产高清在线一区二区三| 床上黄色一级片| 久久久久国产一级毛片高清牌| 亚洲av片天天在线观看| 人人妻人人澡欧美一区二区| 日韩大码丰满熟妇| 免费在线观看日本一区| 国产真人三级小视频在线观看| 久久香蕉激情| av福利片在线观看| 国产精品爽爽va在线观看网站| 国产av一区在线观看免费| 久久久久亚洲av毛片大全| a级毛片在线看网站| 中文字幕高清在线视频| 一进一出好大好爽视频| 精品无人区乱码1区二区| 一级黄色大片毛片| 狠狠狠狠99中文字幕| 国产1区2区3区精品| 欧美午夜高清在线| 在线看三级毛片| 成人特级黄色片久久久久久久| 久久中文字幕人妻熟女| 亚洲av成人精品一区久久| 中文字幕人成人乱码亚洲影| 欧美日韩福利视频一区二区| 亚洲精品一卡2卡三卡4卡5卡| 舔av片在线| 免费在线观看日本一区| 日本 av在线| 亚洲成a人片在线一区二区| 成人特级黄色片久久久久久久| 国产伦人伦偷精品视频| 欧美最黄视频在线播放免费| 国产高清激情床上av| 中文字幕精品亚洲无线码一区| 免费在线观看影片大全网站| 色精品久久人妻99蜜桃| 国产伦在线观看视频一区| 亚洲精品美女久久av网站| 欧美另类亚洲清纯唯美| 男人舔奶头视频| 青草久久国产| 婷婷亚洲欧美| 天天添夜夜摸| 日韩大尺度精品在线看网址| 午夜免费观看网址| 久久久久国产一级毛片高清牌| 999久久久国产精品视频| 久久久久国内视频| 99国产精品一区二区三区| 欧美日韩中文字幕国产精品一区二区三区| 国产aⅴ精品一区二区三区波| 国产免费av片在线观看野外av| 国产精品永久免费网站| 麻豆成人av在线观看| 一级片免费观看大全| 熟妇人妻久久中文字幕3abv| 久久精品综合一区二区三区| 黄色a级毛片大全视频| 国产精品av视频在线免费观看| 成人国语在线视频| 97碰自拍视频| 久久香蕉激情| 日韩中文字幕欧美一区二区| 日本a在线网址| 日韩欧美在线乱码| 亚洲精品美女久久av网站| 91大片在线观看| 婷婷六月久久综合丁香| 精品日产1卡2卡| 久久久国产欧美日韩av| 三级毛片av免费| 成人av一区二区三区在线看| 高清毛片免费观看视频网站| 色老头精品视频在线观看| 制服丝袜大香蕉在线| 亚洲国产欧美网| 19禁男女啪啪无遮挡网站| 99热这里只有是精品50| 丰满人妻一区二区三区视频av | 午夜精品一区二区三区免费看| 久久99热这里只有精品18| 午夜福利成人在线免费观看| 午夜福利18| 老司机午夜十八禁免费视频| 成人av一区二区三区在线看| 午夜精品久久久久久毛片777| 精品国产亚洲在线| 又大又爽又粗| 国产精品98久久久久久宅男小说| 亚洲七黄色美女视频| 搞女人的毛片| 99国产综合亚洲精品| 亚洲午夜理论影院| 亚洲无线在线观看| 一区二区三区国产精品乱码| 琪琪午夜伦伦电影理论片6080| 免费在线观看日本一区| 国产黄a三级三级三级人| 欧美久久黑人一区二区| 丰满的人妻完整版| 国产精品 欧美亚洲| 久久精品aⅴ一区二区三区四区| 人妻久久中文字幕网| 久久国产乱子伦精品免费另类| 99久久精品国产亚洲精品| 欧美又色又爽又黄视频| 亚洲人成电影免费在线| 精品电影一区二区在线| 非洲黑人性xxxx精品又粗又长| 大型黄色视频在线免费观看| 久99久视频精品免费| 国内少妇人妻偷人精品xxx网站 | 久久精品91无色码中文字幕| 国产乱人伦免费视频| 一级片免费观看大全| 亚洲色图av天堂| 亚洲aⅴ乱码一区二区在线播放 | 少妇裸体淫交视频免费看高清 | 午夜影院日韩av| 欧美3d第一页| 成人三级黄色视频| 久久精品国产99精品国产亚洲性色| 人人妻人人澡欧美一区二区| e午夜精品久久久久久久| 国产男靠女视频免费网站| 最新在线观看一区二区三区| 亚洲成人中文字幕在线播放| 国产成人欧美在线观看| 亚洲成人中文字幕在线播放| 美女扒开内裤让男人捅视频| 男女做爰动态图高潮gif福利片| 日本三级黄在线观看| 亚洲激情在线av| 九色国产91popny在线| 日本精品一区二区三区蜜桃| 看免费av毛片| 午夜激情福利司机影院| 伊人久久大香线蕉亚洲五| 夜夜夜夜夜久久久久| 国产一区二区在线观看日韩 | 国语自产精品视频在线第100页| 欧美成人一区二区免费高清观看 | 激情在线观看视频在线高清| 亚洲五月天丁香| 国产视频一区二区在线看| 麻豆国产av国片精品| 国语自产精品视频在线第100页|