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

    考慮發(fā)電機(jī)閥點(diǎn)效應(yīng)的改進(jìn)拉格朗日松弛有功調(diào)度策略

    2021-12-29 07:26:46陳建華于希娟
    電力系統(tǒng)自動(dòng)化 2021年24期
    關(guān)鍵詞:拉格朗調(diào)度效應(yīng)

    陳建華,丁 冬,于希娟,丁 紅,忻 煜

    (1. 國(guó)網(wǎng)冀北電力有限公司,北京市 100053;2. 國(guó)網(wǎng)北京市電力公司,北京市 100031;3. 國(guó)網(wǎng)北京市電力公司電力科學(xué)研究院,北京市 100075)

    0 引言

    經(jīng)濟(jì)調(diào)度問題是電力系統(tǒng)運(yùn)行中最重要的優(yōu)化問題之一[1-2],其目的是在保證電力系統(tǒng)安全運(yùn)行的前提下,通過優(yōu)化調(diào)整機(jī)組出力實(shí)現(xiàn)系統(tǒng)運(yùn)行經(jīng)濟(jì)性最優(yōu)的目的[3-4]。目前,經(jīng)濟(jì)調(diào)度運(yùn)行中,一般認(rèn)為機(jī)組運(yùn)行成本與機(jī)組出力之間為簡(jiǎn)單的二次多項(xiàng)式關(guān)系。但實(shí)際上,對(duì)大型發(fā)電機(jī)組來說,其進(jìn)氣閥突然開啟時(shí)出現(xiàn)的拔絲現(xiàn)象會(huì)在機(jī)組耗量曲線上疊加一個(gè)脈動(dòng)效果,產(chǎn)生閥點(diǎn)效應(yīng),從而導(dǎo)致經(jīng)濟(jì)調(diào)度問題呈現(xiàn)非光滑和非凸特性,增加了模型的復(fù)雜度,給優(yōu)化求解帶來較大難度[5]。

    以往對(duì)于考慮閥點(diǎn)效應(yīng)的有功調(diào)度問題的研究主要以啟發(fā)式算法為主,包括粒子群優(yōu)化算法[6-7]、自組織遷移算法、人工蜂群算法、人工免疫算法、共生生物搜索算法、和聲搜索算法[8-14]等。部分研究將2 種或多種算法進(jìn)行混合求解[15-19],其中先由一種方法進(jìn)行初步搜索,以縮小搜索范圍,然后用另外一種方法進(jìn)行精細(xì)搜索,以達(dá)到最優(yōu)。啟發(fā)式算法的優(yōu)點(diǎn)在于對(duì)優(yōu)化問題的形狀要求較少,對(duì)于非凸非光滑優(yōu)化問題也能獲得優(yōu)化解。但啟發(fā)式算法每一次優(yōu)化結(jié)果均不相同,存在前后兩次結(jié)果差別較大的情況,且結(jié)果不能復(fù)現(xiàn),因此,目前尚難以在實(shí)際現(xiàn)場(chǎng)應(yīng)用。

    在數(shù)學(xué)優(yōu)化領(lǐng)域,對(duì)于此類非凸非光滑優(yōu)化問題,一般采用分段線性化方法進(jìn)行求解,但由此會(huì)產(chǎn)生2 個(gè)問題:一是優(yōu)化模型中引入了整數(shù)變量,從而變?yōu)榛旌险麛?shù)規(guī)劃問題,降低了求解效率;二是對(duì)于存在多個(gè)局部最優(yōu)點(diǎn)的優(yōu)化問題,可能會(huì)出現(xiàn)最優(yōu)解來回跳變,從而導(dǎo)致優(yōu)化過程不收斂問題。

    本文提出一種改進(jìn)的拉格朗日松弛有功調(diào)度算法,以考慮發(fā)電機(jī)閥點(diǎn)效應(yīng)影響。提出一種二次多項(xiàng)式分段擬合算法及基于近端梯度法的增廣求解策略,能夠在確保結(jié)果最優(yōu)性的基礎(chǔ)上加快收斂速度。

    1 考慮閥點(diǎn)效應(yīng)的有功經(jīng)濟(jì)調(diào)度問題模型

    有功經(jīng)濟(jì)調(diào)度問題模型的目標(biāo)函數(shù)和約束條件分析如下。

    1)目標(biāo)函數(shù)

    一般以火電機(jī)組煤耗最小作為有功經(jīng)濟(jì)調(diào)度問題的優(yōu)化目標(biāo)函數(shù),如式(1)所示。

    式中:f(pi)為火電機(jī)組i的煤耗成本,pi為機(jī)組i的有功出力;ai,bi,ci為機(jī)組i的煤耗成本系數(shù)。

    當(dāng)考慮閥點(diǎn)效應(yīng)時(shí),式(1)將轉(zhuǎn)變?yōu)槭剑?)。

    式中:n為系統(tǒng)中所有火電機(jī)組數(shù)量;ei和fi為常數(shù);為第i臺(tái)機(jī)組的有功出力下限。

    2)約束條件

    發(fā)電負(fù)荷平衡約束為:

    式中:D為系統(tǒng)負(fù)荷需求預(yù)測(cè)值。

    發(fā)電機(jī)出力限值約束為:

    輸電斷面容量約束為:

    可以看出,當(dāng)考慮常規(guī)機(jī)組閥點(diǎn)效應(yīng)時(shí),優(yōu)化目標(biāo)將由傳統(tǒng)的凸優(yōu)化問題轉(zhuǎn)變?yōu)榉峭狗枪饣瑑?yōu)化問題。由于目前未有針對(duì)這一問題的有效求解策略,使得經(jīng)濟(jì)調(diào)度問題優(yōu)化求解困難。

    2 改進(jìn)拉格朗日松弛法求解算法

    2.1 整體求解流程

    對(duì)第1 章中的有功經(jīng)濟(jì)調(diào)度模型,采用拉格朗日松弛法進(jìn)行求解,對(duì)應(yīng)的拉格朗日對(duì)偶問題如式(6)所示。

    式中:C為常數(shù)項(xiàng)。

    可以看出,式(7)可以看成是多個(gè)單機(jī)約束的優(yōu)化子問題之和。因此,原優(yōu)化問題可以表示為多個(gè)并行的單機(jī)優(yōu)化子問題,每個(gè)子問題中的變量及約束僅與一臺(tái)機(jī)組相關(guān),而與其他機(jī)組無關(guān)。子問題的表達(dá)式如式(8)所示。

    對(duì)于傳統(tǒng)有功調(diào)度問題而言,f(pi)為二次多項(xiàng)式函數(shù),因此,子問題為二次規(guī)劃問題,最優(yōu)解可以較容易獲得。

    在子問題獲得最優(yōu)解后,通過次梯度法更新主問題的拉格朗日乘子,并不斷對(duì)主問題、子問題進(jìn)行迭代,即可逐漸逼近最優(yōu)解。

    2.2 考慮閥點(diǎn)效應(yīng)的改進(jìn)拉格朗日松弛算法

    當(dāng)考慮閥點(diǎn)效應(yīng)時(shí),f(pi)轉(zhuǎn)變?yōu)槭剑?)的非凸非光滑函數(shù)形式,如圖1 所示,導(dǎo)致傳統(tǒng)的數(shù)值優(yōu)化算法不再適用。

    圖1 考慮閥點(diǎn)效應(yīng)的機(jī)組煤耗成本與有功出力關(guān)系Fig.1 Relationship between coal consumption cost and active power output of generator considering valve-point effect

    式中:αi,k、βi,k、γi,k為第i臺(tái)機(jī)組第k分段二次擬合多項(xiàng)式系數(shù)。

    以圖1 所示的發(fā)電機(jī)有功出力曲線為例,二次多項(xiàng)式擬合的誤差曲線如圖2 所示。其中,擬合誤差=|擬合值- 真實(shí)值|/機(jī)組額定出力下的煤耗成本。

    圖2 擬合誤差曲線Fig.2 Fitting error curve

    由圖2 可知,二次擬合最大誤差為0.25%左右,能夠滿足實(shí)際應(yīng)用需要。

    將擬合后的二次多項(xiàng)式(9)代入式(2)可以得到火電機(jī)組煤耗成本為:

    將式(14)代入式(8),可以得到此時(shí)拉格朗日子問題的表達(dá)式為:

    式(15)為分段函數(shù)形式,其中,對(duì)第k分段來說,其優(yōu)化目標(biāo)函數(shù)為二次函數(shù),開口方向取決于二次項(xiàng)系數(shù)大小,對(duì)應(yīng)的二次項(xiàng)系數(shù)為:

    因此,當(dāng)ai>4ei f2i/π2時(shí),最優(yōu)輸出功率為:

    當(dāng)ai<4ei f2i/π2時(shí),最優(yōu)輸出功率為:

    實(shí)際上,對(duì)絕大多數(shù)火電機(jī)組來說,ai取值一般在10-4~10-3之間,而ei的取值一般為102量級(jí)、fi的取值一般為10-2量級(jí),因此,ai<4ei f2i/π2總是成立,優(yōu)化過程中一般只需要考慮式(18)一種情況即可。式(18)表明,對(duì)于開口向下的二次函數(shù),其最小值將在某一邊界處取得。因此,優(yōu)化問題簡(jiǎn)化為只考慮每臺(tái)機(jī)組少數(shù)幾個(gè)閥點(diǎn)即可,此時(shí)需要松弛發(fā)電負(fù)荷平衡約束(v≥0),確保問題能夠收斂。

    2.3 基于近端梯度法的改進(jìn)子問題求解策略

    拉格朗日松弛法對(duì)原問題進(jìn)行對(duì)偶求解做法的一個(gè)主要缺點(diǎn)是要求原問題具有凸結(jié)構(gòu),否則對(duì)偶問題的最優(yōu)解與原問題之間不僅存在對(duì)偶間隙,而且導(dǎo)致對(duì)偶問題收斂性較差,容易出現(xiàn)振蕩不收斂情況。

    針對(duì)子問題的非連續(xù)可微特性,本文提出一種基于近端梯度法的改進(jìn)子問題求解方法。近端梯度法是一種特殊的梯度下降算法,主要用于求解目標(biāo)函數(shù)不可微的最優(yōu)化問題。如果目標(biāo)函數(shù)在某些點(diǎn)是不可微的,那么該點(diǎn)的梯度無法求解,傳統(tǒng)的梯度下降法也就無法使用。

    近端梯度法的核心思想是通過proximal 映射,將不可微函數(shù)轉(zhuǎn)化為易求解的proximal 映射函數(shù),從而實(shí)現(xiàn)近似求解。假設(shè)約束函數(shù)f(x)的定義域?yàn)閁,則自變量x的proximal 映射為:

    本文在上式基礎(chǔ)上,通過在目標(biāo)函數(shù)中引入懲罰因子?,在不影響原優(yōu)化問題的優(yōu)化結(jié)果的基礎(chǔ)上,加快問題收斂速度。

    式中:?為懲罰項(xiàng)系數(shù);y為輔助變量。為加快收斂速度,在第λ次迭代時(shí),本文令y取為第λ-1 次迭代時(shí)對(duì)應(yīng)的有功出力。

    可以看出,式(20)優(yōu)化問題的最優(yōu)解p*i與式(15)優(yōu)化問題最優(yōu)解相同,且在最優(yōu)解處,有y=p*i。因此,二次懲罰項(xiàng)的加入不影響原優(yōu)化問題的優(yōu)化結(jié)果。

    對(duì)于式(20)的優(yōu)化問題,其目標(biāo)函數(shù)的Hessian矩陣為:

    當(dāng)?取值過小時(shí),將導(dǎo)致優(yōu)化問題Hesssian 矩陣條件數(shù)過大,從而導(dǎo)致優(yōu)化問題病態(tài);而當(dāng)?取值過大時(shí),優(yōu)化問題凸性不足,將導(dǎo)致優(yōu)化問題難以收斂。因此,合適的懲罰項(xiàng)系數(shù)選取對(duì)于保證優(yōu)化問題的最優(yōu)性及高效性具有重要意義。一般來說,為保證問題的收斂性,一般要求Hessian 矩陣為正定矩陣,即:

    由式(22)、式(23)可知,?的取值區(qū)間建議為:

    以2.2 節(jié)中發(fā)電負(fù)荷平衡約束松弛后獲得的解作為優(yōu)化初值,用式(20)替換式(15),通過對(duì)每一分段分別求解并比較各分段區(qū)間優(yōu)化解的經(jīng)濟(jì)性即可得到最優(yōu)煤耗及對(duì)應(yīng)的發(fā)電機(jī)組出力。

    3 算例分析

    3.1 40 機(jī)測(cè)試系統(tǒng)仿真結(jié)果

    選取文獻(xiàn)[20]中的40 機(jī)測(cè)試系統(tǒng)作為研究對(duì)象,系統(tǒng)負(fù)荷需求為10 500 MW,懲罰項(xiàng)系數(shù)?取0.01。

    不考慮閥點(diǎn)效應(yīng)時(shí),有功調(diào)度問題轉(zhuǎn)化為傳統(tǒng)的二次規(guī)劃問題,可以采用傳統(tǒng)的拉格朗日松弛法快速求解,得到最小運(yùn)行成本為118 660.2 美元。

    考慮閥點(diǎn)效應(yīng)時(shí),采用本文方法及與其他方法優(yōu)化結(jié)果比較如表1 所示。

    表1 幾種方法優(yōu)化結(jié)果比較(40 機(jī)測(cè)試系統(tǒng))Table 1 Comparison of optimization results of several methods(40-machine test system)

    由表1 可見,當(dāng)考慮閥點(diǎn)效應(yīng)時(shí),系統(tǒng)運(yùn)行成本由118 660.2 美元至少增加到121 412.5 美元,升高了2.32%,可見發(fā)電機(jī)閥點(diǎn)效應(yīng)特性對(duì)系統(tǒng)調(diào)度運(yùn)行的經(jīng)濟(jì)性具有明顯影響。

    同時(shí),本文方法計(jì)算結(jié)果為121 459.6 美元,僅比全局最優(yōu)解(121 412.5 美元)高0.03%左右,但計(jì)算時(shí)間僅為1.34 s,計(jì)算效率明顯高于混合整數(shù)規(guī)劃法。同時(shí),本文方法計(jì)算結(jié)果也明顯優(yōu)于其他6 種計(jì)算方法,且計(jì)算所需時(shí)間最短、效率最高。本文方法的收斂特性曲線如圖3 所示,經(jīng)過200 次左右迭代,本文方法能夠快速逼近最優(yōu)點(diǎn)。

    圖3 所提算法收斂特性Fig.3 Convergence characteristic of proposed algorithm

    3.2 48 機(jī)測(cè)試系統(tǒng)仿真結(jié)果

    以文獻(xiàn)[24]中的48 機(jī)測(cè)試系統(tǒng)為研究對(duì)象,系統(tǒng)包括26 臺(tái)傳統(tǒng)火力發(fā)電機(jī)組、12 臺(tái)熱電聯(lián)產(chǎn)機(jī)組及10 臺(tái)供熱機(jī)組。系統(tǒng)負(fù)荷需求為4 700 MW,供熱需求為2 500 MWth,懲罰項(xiàng)系數(shù)?取值為0.9。

    采用本文方法與其他方法的優(yōu)化結(jié)果比較如表2 所示,每一臺(tái)機(jī)組的優(yōu)化結(jié)果比較如附錄A 表A1所示。

    表2 幾種方法優(yōu)化結(jié)果比較(48 機(jī)測(cè)試系統(tǒng))Table 2 Comparison of optimization results of several methods(48-machine test system)

    由表2 可見,受優(yōu)化模型目標(biāo)函數(shù)、約束條件均存在強(qiáng)非凸性的影響,傳統(tǒng)拉格朗日松弛法無法收斂,采用混合整數(shù)規(guī)劃法可以得到最優(yōu)的優(yōu)化結(jié)果,但算法運(yùn)行時(shí)間較長(zhǎng),為84.82 s,效率較低。本文方法煤耗成本僅為115 806.55 美元,僅高于全局最優(yōu)解0.12%。與其他幾種方法相比,本文方法煤耗成本也最低。同時(shí),本文方法運(yùn)行時(shí)間僅為6.29 s,遠(yuǎn)小于其他幾種方法,計(jì)算效率高。同時(shí),對(duì)于現(xiàn)場(chǎng)實(shí)際應(yīng)用來說,本文方法具有可重復(fù)性,且只需要運(yùn)行一次即可獲得最優(yōu)解,而啟發(fā)式算法需要反復(fù)運(yùn)行多次,且一般不能復(fù)現(xiàn)結(jié)果。

    本文方法的收斂特性曲線如附錄B 圖B1和圖B2所示??梢姡?jīng)過600 次左右迭代,本文方法能夠逐漸逼近最優(yōu)點(diǎn)。根據(jù)式(23)可知,要保證結(jié)果的收斂性,懲罰項(xiàng)系數(shù)?的取值范圍應(yīng)為[0,1.749]。對(duì)不同?取值下的優(yōu)化結(jié)果進(jìn)行分析,如表3 所示。

    表3 不同懲罰項(xiàng)系數(shù)時(shí)的優(yōu)化結(jié)果比較Table 3 Comparison of optimization results with different penalty term coefficients

    由表3 可見,當(dāng)?在允許區(qū)間范圍內(nèi)取值時(shí),優(yōu)化結(jié)果及計(jì)算效率差別不大;當(dāng)?取值超出允許區(qū)間時(shí),優(yōu)化過程不收斂。特別的,當(dāng)?取為∞時(shí),優(yōu)化模型轉(zhuǎn)變?yōu)閭鹘y(tǒng)的拉格朗日松弛算法,此時(shí),優(yōu)化過程不收斂。

    進(jìn)一步,考慮不同負(fù)荷需求水平及供熱需求水平,本文方法的收斂性如附錄B 圖B3 所示。由圖B3 可見,在不同負(fù)荷需求及供熱需求水平下,本文方法均能夠快速收斂,證明了本文方法的魯棒性。同時(shí)可以看出,負(fù)荷水平變化時(shí)的收斂曲線具有一定的振蕩性,與不同供熱需求水平下的收斂曲線具有較大差異,這是由于考慮常規(guī)機(jī)組閥點(diǎn)效應(yīng)時(shí)優(yōu)化目標(biāo)函數(shù)的非凸非光滑特性導(dǎo)致的。

    4 結(jié)語

    考慮閥點(diǎn)效應(yīng)的有功調(diào)度問題,以往的研究主要集中于啟發(fā)式算法。本文提出一種基于數(shù)學(xué)規(guī)劃方法的求解策略,通過對(duì)拉格朗日松弛法進(jìn)行改進(jìn),并結(jié)合近端梯度算法,消除了發(fā)電機(jī)閥點(diǎn)效應(yīng)導(dǎo)致的非凸非光滑特性的影響。同時(shí),提出一種二次多項(xiàng)式分段擬合及基于近端梯度法的增廣優(yōu)化求解策略。40 機(jī)及48 機(jī)兩個(gè)測(cè)試系統(tǒng)上的優(yōu)化結(jié)果表明:本文方法收斂性好、運(yùn)行效率較高,較為適合在現(xiàn)場(chǎng)實(shí)際應(yīng)用。

    同時(shí),需要指出的是,考慮閥點(diǎn)效應(yīng)的有功調(diào)度問題為非凸優(yōu)化問題,由于對(duì)偶間隙的存在,拉格朗日松弛法在求解此類問題時(shí),只能得到局部最優(yōu)解。下一步需要對(duì)此進(jìn)行改進(jìn),結(jié)合其他數(shù)學(xué)規(guī)劃算法的優(yōu)點(diǎn),尋求能夠快速獲得全局最優(yōu)解的高效算法。

    附錄見本刊網(wǎng)絡(luò)版(http://www.aeps-info.com/aeps/ch/index.aspx),掃英文摘要后二維碼可以閱讀網(wǎng)絡(luò)全文。

    猜你喜歡
    拉格朗調(diào)度效應(yīng)
    鈾對(duì)大型溞的急性毒性效應(yīng)
    懶馬效應(yīng)
    《調(diào)度集中系統(tǒng)(CTC)/列車調(diào)度指揮系統(tǒng)(TDCS)維護(hù)手冊(cè)》正式出版
    一種基于負(fù)載均衡的Kubernetes調(diào)度改進(jìn)算法
    虛擬機(jī)實(shí)時(shí)遷移調(diào)度算法
    Nearly Kaehler流形S3×S3上的切觸拉格朗日子流形
    應(yīng)變效應(yīng)及其應(yīng)用
    拉格朗日代數(shù)方程求解中的置換思想
    基于拉格朗日的IGS精密星歷和鐘差插值分析
    拉格朗日點(diǎn)
    太空探索(2014年3期)2014-07-10 14:59:39
    午夜影院日韩av| 久久人人精品亚洲av| 美女黄网站色视频| 99久久无色码亚洲精品果冻| 欧美一级毛片孕妇| 手机成人av网站| 午夜福利18| a级毛片a级免费在线| 在线观看免费日韩欧美大片| 男女那种视频在线观看| 精品一区二区三区视频在线观看免费| 精品国产美女av久久久久小说| 我要搜黄色片| 精品欧美一区二区三区在线| 国产精品国产高清国产av| 丁香六月欧美| 精品电影一区二区在线| 变态另类丝袜制服| 非洲黑人性xxxx精品又粗又长| 一级毛片女人18水好多| 精品电影一区二区在线| 国产三级黄色录像| 午夜福利在线观看吧| 精品午夜福利视频在线观看一区| 国产亚洲精品一区二区www| 18禁观看日本| 国产精品爽爽va在线观看网站| 两性夫妻黄色片| 国产精品 欧美亚洲| 99国产精品一区二区三区| www.自偷自拍.com| 午夜免费观看网址| 亚洲精华国产精华精| 91麻豆av在线| 欧美日韩精品网址| 国内揄拍国产精品人妻在线| 色精品久久人妻99蜜桃| 91麻豆av在线| 精品少妇一区二区三区视频日本电影| 国产成年人精品一区二区| 99久久无色码亚洲精品果冻| 最近最新中文字幕大全电影3| 一二三四在线观看免费中文在| 一区福利在线观看| 久久精品91无色码中文字幕| 久久国产精品人妻蜜桃| 人妻夜夜爽99麻豆av| 国产成年人精品一区二区| 91国产中文字幕| 制服人妻中文乱码| 一区二区三区高清视频在线| 亚洲精品粉嫩美女一区| aaaaa片日本免费| 国产亚洲精品综合一区在线观看 | 18禁黄网站禁片午夜丰满| 国产探花在线观看一区二区| 国产亚洲av高清不卡| 五月伊人婷婷丁香| 久久中文字幕一级| 亚洲男人的天堂狠狠| 又黄又爽又免费观看的视频| 男女下面进入的视频免费午夜| 精品免费久久久久久久清纯| 亚洲第一欧美日韩一区二区三区| 国产精品98久久久久久宅男小说| 久久久久久久久久黄片| 久久午夜亚洲精品久久| 国内久久婷婷六月综合欲色啪| 男女那种视频在线观看| 999久久久国产精品视频| 亚洲成人精品中文字幕电影| 精品久久久久久久毛片微露脸| 99热6这里只有精品| 色综合婷婷激情| 人人妻,人人澡人人爽秒播| 亚洲精品美女久久久久99蜜臀| 后天国语完整版免费观看| 人妻丰满熟妇av一区二区三区| 亚洲精品中文字幕在线视频| 91成年电影在线观看| 国产熟女午夜一区二区三区| 国产亚洲精品综合一区在线观看 | 国内久久婷婷六月综合欲色啪| 在线a可以看的网站| 1024香蕉在线观看| 99精品久久久久人妻精品| 日日爽夜夜爽网站| 日韩欧美精品v在线| 一个人免费在线观看电影 | 日本熟妇午夜| 日韩欧美精品v在线| av视频在线观看入口| 亚洲男人的天堂狠狠| 亚洲专区国产一区二区| 免费电影在线观看免费观看| 99久久国产精品久久久| 又黄又粗又硬又大视频| 午夜免费观看网址| 非洲黑人性xxxx精品又粗又长| 中文资源天堂在线| 男女做爰动态图高潮gif福利片| 日韩高清综合在线| 成人欧美大片| 久久久久亚洲av毛片大全| 最近视频中文字幕2019在线8| 成人国语在线视频| 国产视频一区二区在线看| 亚洲,欧美精品.| 真人一进一出gif抽搐免费| 成人一区二区视频在线观看| 午夜福利18| 两个人视频免费观看高清| 午夜成年电影在线免费观看| 精品不卡国产一区二区三区| 亚洲中文av在线| 97超级碰碰碰精品色视频在线观看| 欧美日韩国产亚洲二区| 国产蜜桃级精品一区二区三区| 成人手机av| 成人一区二区视频在线观看| 五月伊人婷婷丁香| 一边摸一边抽搐一进一小说| 精品久久久久久久人妻蜜臀av| 午夜福利在线在线| 51午夜福利影视在线观看| 国产精品久久电影中文字幕| 人妻夜夜爽99麻豆av| 亚洲av电影在线进入| 欧美丝袜亚洲另类 | 精品人妻1区二区| 我要搜黄色片| 黑人欧美特级aaaaaa片| 亚洲乱码一区二区免费版| 欧美成人性av电影在线观看| 色尼玛亚洲综合影院| 欧美三级亚洲精品| 制服丝袜大香蕉在线| 午夜福利成人在线免费观看| 欧美日韩黄片免| 久久精品综合一区二区三区| 亚洲avbb在线观看| 天天躁夜夜躁狠狠躁躁| 国产成人精品久久二区二区91| 亚洲第一电影网av| 又粗又爽又猛毛片免费看| 男女之事视频高清在线观看| 亚洲av成人精品一区久久| 18禁观看日本| 色精品久久人妻99蜜桃| 亚洲一区二区三区色噜噜| 欧美一区二区国产精品久久精品 | 日韩欧美 国产精品| 天天添夜夜摸| 最近最新中文字幕大全电影3| 小说图片视频综合网站| 国产一区二区三区在线臀色熟女| 正在播放国产对白刺激| 日韩精品青青久久久久久| 男女之事视频高清在线观看| 人妻丰满熟妇av一区二区三区| av中文乱码字幕在线| 岛国在线免费视频观看| 免费观看精品视频网站| 99re在线观看精品视频| 1024手机看黄色片| 亚洲国产精品合色在线| 久久久久久久午夜电影| 国产av一区在线观看免费| 三级国产精品欧美在线观看 | 99久久国产精品久久久| 老熟妇乱子伦视频在线观看| 亚洲国产欧美人成| 一卡2卡三卡四卡精品乱码亚洲| 国产精品久久久久久久电影 | 国产欧美日韩精品亚洲av| 在线看三级毛片| 免费无遮挡裸体视频| 亚洲av成人av| 夜夜夜夜夜久久久久| 777久久人妻少妇嫩草av网站| 欧美黑人欧美精品刺激| 欧美另类亚洲清纯唯美| 999久久久国产精品视频| 欧美成人性av电影在线观看| 国产1区2区3区精品| 巨乳人妻的诱惑在线观看| 午夜福利视频1000在线观看| 中文资源天堂在线| 亚洲熟妇中文字幕五十中出| 精品久久久久久成人av| 久久久久久九九精品二区国产 | 亚洲成人免费电影在线观看| 少妇裸体淫交视频免费看高清 | 美女免费视频网站| 99在线视频只有这里精品首页| 国产亚洲精品一区二区www| 午夜免费观看网址| 精品午夜福利视频在线观看一区| 欧美日韩黄片免| 国产亚洲精品久久久久5区| 国产人伦9x9x在线观看| 日本五十路高清| a级毛片在线看网站| 一级毛片女人18水好多| 午夜福利欧美成人| 午夜老司机福利片| 亚洲精品在线美女| 亚洲五月天丁香| 亚洲欧美精品综合久久99| 757午夜福利合集在线观看| 国产私拍福利视频在线观看| 久久性视频一级片| 国产成+人综合+亚洲专区| 久久久久久久久久黄片| 又爽又黄无遮挡网站| 日韩中文字幕欧美一区二区| 国产精品精品国产色婷婷| 高清在线国产一区| 日韩精品青青久久久久久| 色精品久久人妻99蜜桃| 免费在线观看视频国产中文字幕亚洲| 国产一区二区在线av高清观看| 99久久综合精品五月天人人| 国产精品久久久久久亚洲av鲁大| 香蕉国产在线看| 少妇熟女aⅴ在线视频| 久久精品影院6| 欧美日韩亚洲国产一区二区在线观看| 精品人妻1区二区| 久久久久九九精品影院| 国产精品亚洲美女久久久| 99热只有精品国产| 少妇的丰满在线观看| 久久久久久免费高清国产稀缺| 国产黄色小视频在线观看| 夜夜躁狠狠躁天天躁| av国产免费在线观看| 精品熟女少妇八av免费久了| 全区人妻精品视频| 香蕉丝袜av| 精品第一国产精品| 国产蜜桃级精品一区二区三区| 99久久精品国产亚洲精品| 日本黄色视频三级网站网址| 日日干狠狠操夜夜爽| 看免费av毛片| 国产私拍福利视频在线观看| 精品久久蜜臀av无| 午夜福利高清视频| 亚洲成av人片免费观看| 欧美绝顶高潮抽搐喷水| 90打野战视频偷拍视频| 亚洲在线自拍视频| 在线永久观看黄色视频| 亚洲av电影不卡..在线观看| 亚洲精品一区av在线观看| 毛片女人毛片| 麻豆成人av在线观看| 国产精品av久久久久免费| 99热这里只有是精品50| 亚洲美女视频黄频| 欧美中文日本在线观看视频| 日韩精品青青久久久久久| 十八禁人妻一区二区| 亚洲欧美日韩东京热| 国产三级黄色录像| 日韩大尺度精品在线看网址| 两个人免费观看高清视频| 在线看三级毛片| 不卡一级毛片| 久久精品成人免费网站| 亚洲av五月六月丁香网| 久久中文字幕人妻熟女| 啦啦啦观看免费观看视频高清| 久久亚洲精品不卡| 人妻久久中文字幕网| 久久久久亚洲av毛片大全| 久久午夜亚洲精品久久| 日本黄色视频三级网站网址| 日韩大码丰满熟妇| 三级男女做爰猛烈吃奶摸视频| 久久中文字幕一级| 欧美乱码精品一区二区三区| 欧美一区二区国产精品久久精品 | 色哟哟哟哟哟哟| 久久久久久久午夜电影| 男人舔女人的私密视频| 国内揄拍国产精品人妻在线| cao死你这个sao货| 日本撒尿小便嘘嘘汇集6| 亚洲一区二区三区不卡视频| 国产高清激情床上av| 午夜精品久久久久久毛片777| 亚洲精品av麻豆狂野| 午夜福利18| 精品久久久久久久末码| 精品高清国产在线一区| 久久精品夜夜夜夜夜久久蜜豆 | 午夜福利欧美成人| 亚洲黑人精品在线| 99精品在免费线老司机午夜| 天堂影院成人在线观看| 熟妇人妻久久中文字幕3abv| 久久精品国产亚洲av香蕉五月| 亚洲色图 男人天堂 中文字幕| 欧美乱码精品一区二区三区| 2021天堂中文幕一二区在线观| 亚洲精华国产精华精| e午夜精品久久久久久久| 三级男女做爰猛烈吃奶摸视频| 欧美黑人欧美精品刺激| 麻豆国产av国片精品| 亚洲专区字幕在线| 中文字幕人妻丝袜一区二区| 国产97色在线日韩免费| 国产视频一区二区在线看| www.999成人在线观看| 热99re8久久精品国产| 婷婷亚洲欧美| 亚洲精品国产精品久久久不卡| 99久久国产精品久久久| 亚洲va日本ⅴa欧美va伊人久久| 最好的美女福利视频网| 99国产精品一区二区蜜桃av| 丰满人妻熟妇乱又伦精品不卡| 日本熟妇午夜| 欧美日本亚洲视频在线播放| or卡值多少钱| 欧美日韩福利视频一区二区| 在线观看舔阴道视频| 黄色女人牲交| tocl精华| 免费在线观看视频国产中文字幕亚洲| www.999成人在线观看| 欧美极品一区二区三区四区| 亚洲午夜精品一区,二区,三区| 久久久国产精品麻豆| 午夜福利成人在线免费观看| 亚洲人成网站在线播放欧美日韩| videosex国产| svipshipincom国产片| 国产99白浆流出| 一夜夜www| 日本在线视频免费播放| 黄色 视频免费看| 国产午夜精品久久久久久| 一级作爱视频免费观看| 免费在线观看视频国产中文字幕亚洲| 97人妻精品一区二区三区麻豆| 欧美成狂野欧美在线观看| 久久午夜综合久久蜜桃| 少妇粗大呻吟视频| 国产伦人伦偷精品视频| 日本三级黄在线观看| 女同久久另类99精品国产91| 十八禁网站免费在线| 美女 人体艺术 gogo| 国产免费男女视频| 黄片小视频在线播放| 欧美日韩福利视频一区二区| 国产99久久九九免费精品| 欧美激情久久久久久爽电影| 美女午夜性视频免费| 国内毛片毛片毛片毛片毛片| 国产av又大| 国产伦在线观看视频一区| 国产免费av片在线观看野外av| 午夜免费激情av| 看片在线看免费视频| 午夜福利欧美成人| 操出白浆在线播放| 精品一区二区三区视频在线观看免费| 好男人在线观看高清免费视频| 成人特级黄色片久久久久久久| 嫁个100分男人电影在线观看| 精品不卡国产一区二区三区| 毛片女人毛片| 国产精品综合久久久久久久免费| 午夜激情福利司机影院| 国产三级在线视频| 老司机福利观看| 亚洲国产看品久久| √禁漫天堂资源中文www| 久久精品aⅴ一区二区三区四区| 日本三级黄在线观看| 丰满的人妻完整版| 国产精品一区二区免费欧美| 亚洲精品粉嫩美女一区| 91在线观看av| а√天堂www在线а√下载| 欧美性长视频在线观看| 免费人成视频x8x8入口观看| 一级黄色大片毛片| 欧美成狂野欧美在线观看| 人妻夜夜爽99麻豆av| av福利片在线| 亚洲精品一区av在线观看| 欧美黄色淫秽网站| 一区福利在线观看| xxx96com| 国产亚洲欧美在线一区二区| 好看av亚洲va欧美ⅴa在| 一个人免费在线观看电影 | 久久欧美精品欧美久久欧美| 欧美中文日本在线观看视频| av欧美777| 国产片内射在线| 午夜免费观看网址| 成人一区二区视频在线观看| 国产精品精品国产色婷婷| 搞女人的毛片| 亚洲五月天丁香| 亚洲人成网站在线播放欧美日韩| 麻豆国产av国片精品| 亚洲 欧美 日韩 在线 免费| 在线观看午夜福利视频| 少妇粗大呻吟视频| 最新在线观看一区二区三区| 一a级毛片在线观看| 在线播放国产精品三级| 五月玫瑰六月丁香| 中文字幕av在线有码专区| 叶爱在线成人免费视频播放| 亚洲午夜精品一区,二区,三区| 1024视频免费在线观看| 麻豆成人av在线观看| 日韩大码丰满熟妇| 久久久久久人人人人人| 国产亚洲欧美98| 日韩 欧美 亚洲 中文字幕| 国产精品永久免费网站| 男女床上黄色一级片免费看| 色综合婷婷激情| 日本免费a在线| 757午夜福利合集在线观看| 老鸭窝网址在线观看| 亚洲中文字幕日韩| 亚洲人成电影免费在线| 日韩欧美精品v在线| 亚洲男人的天堂狠狠| 宅男免费午夜| 99国产精品99久久久久| 在线观看美女被高潮喷水网站 | 国产精品自产拍在线观看55亚洲| 岛国在线观看网站| 在线观看午夜福利视频| 老熟妇仑乱视频hdxx| 精品无人区乱码1区二区| 亚洲国产精品sss在线观看| 亚洲国产精品久久男人天堂| 亚洲熟妇中文字幕五十中出| 在线观看舔阴道视频| 国产成人精品久久二区二区91| 精品不卡国产一区二区三区| 免费在线观看黄色视频的| 亚洲美女黄片视频| 蜜桃久久精品国产亚洲av| 少妇被粗大的猛进出69影院| 欧美最黄视频在线播放免费| 舔av片在线| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美一区二区三区在线观看| 国产激情欧美一区二区| 午夜亚洲福利在线播放| 女同久久另类99精品国产91| 99精品久久久久人妻精品| 久久人人精品亚洲av| 成年女人毛片免费观看观看9| 一区福利在线观看| 91老司机精品| 亚洲人成77777在线视频| 少妇熟女aⅴ在线视频| www.www免费av| 麻豆久久精品国产亚洲av| 99久久国产精品久久久| 午夜福利18| 国产激情久久老熟女| 免费在线观看完整版高清| 成人一区二区视频在线观看| 性色av乱码一区二区三区2| 国产乱人伦免费视频| 精品欧美国产一区二区三| 午夜免费成人在线视频| 日本在线视频免费播放| АⅤ资源中文在线天堂| 欧美日本视频| 国产高清有码在线观看视频 | 国产精品av久久久久免费| 精品福利观看| 好男人在线观看高清免费视频| 午夜福利在线在线| svipshipincom国产片| 国产精品一区二区三区四区久久| 欧美日韩国产亚洲二区| 中国美女看黄片| 妹子高潮喷水视频| x7x7x7水蜜桃| 久久99热这里只有精品18| 免费观看精品视频网站| 制服丝袜大香蕉在线| 久久香蕉精品热| 操出白浆在线播放| 一a级毛片在线观看| 黑人欧美特级aaaaaa片| 香蕉丝袜av| 亚洲专区国产一区二区| 美女 人体艺术 gogo| 美女午夜性视频免费| 夜夜躁狠狠躁天天躁| 国产精品精品国产色婷婷| av福利片在线观看| 国产精品美女特级片免费视频播放器 | 亚洲最大成人中文| 禁无遮挡网站| 久久久国产欧美日韩av| 级片在线观看| 中文字幕熟女人妻在线| 黄色成人免费大全| 久久久精品国产亚洲av高清涩受| 男女做爰动态图高潮gif福利片| 中文字幕久久专区| av视频在线观看入口| 999精品在线视频| 国产区一区二久久| 欧美黄色片欧美黄色片| 国产熟女xx| 久久精品成人免费网站| 嫩草影视91久久| 99热这里只有精品一区 | 精品免费久久久久久久清纯| 日韩中文字幕欧美一区二区| 成年免费大片在线观看| 国产私拍福利视频在线观看| 色综合婷婷激情| 午夜福利在线观看吧| 亚洲av成人不卡在线观看播放网| 欧美中文综合在线视频| 午夜福利免费观看在线| 精品少妇一区二区三区视频日本电影| 热99re8久久精品国产| 久久精品国产亚洲av香蕉五月| 成人欧美大片| 男女视频在线观看网站免费 | 亚洲avbb在线观看| АⅤ资源中文在线天堂| 亚洲欧美日韩东京热| 真人一进一出gif抽搐免费| 五月伊人婷婷丁香| 中文字幕久久专区| 色播亚洲综合网| 亚洲 国产 在线| 法律面前人人平等表现在哪些方面| 日韩三级视频一区二区三区| 色老头精品视频在线观看| 日韩有码中文字幕| 国产精品香港三级国产av潘金莲| 成人一区二区视频在线观看| 老鸭窝网址在线观看| 亚洲全国av大片| √禁漫天堂资源中文www| 琪琪午夜伦伦电影理论片6080| 一二三四在线观看免费中文在| 亚洲人与动物交配视频| 免费av毛片视频| 一区二区三区国产精品乱码| 这个男人来自地球电影免费观看| 精品人妻1区二区| 两个人看的免费小视频| 国产亚洲精品久久久久5区| 精品国产超薄肉色丝袜足j| 国产区一区二久久| 亚洲男人天堂网一区| 免费在线观看视频国产中文字幕亚洲| 久9热在线精品视频| 日韩 欧美 亚洲 中文字幕| 麻豆一二三区av精品| 久久天躁狠狠躁夜夜2o2o| 国产精品久久久久久亚洲av鲁大| 正在播放国产对白刺激| 国产亚洲精品久久久久5区| 精品国产超薄肉色丝袜足j| 人妻久久中文字幕网| 麻豆国产97在线/欧美 | 久久精品国产亚洲av高清一级| 麻豆久久精品国产亚洲av| 老司机午夜福利在线观看视频| 亚洲精品粉嫩美女一区| 婷婷亚洲欧美| 久久久久性生活片| 国产精品一及| 熟妇人妻久久中文字幕3abv| 毛片女人毛片| 一区福利在线观看| 麻豆一二三区av精品| 一本久久中文字幕| 午夜影院日韩av| 国产99久久九九免费精品| 亚洲黑人精品在线| 此物有八面人人有两片| 日韩欧美国产在线观看| а√天堂www在线а√下载| 老熟妇仑乱视频hdxx| 级片在线观看| 免费观看人在逋| 每晚都被弄得嗷嗷叫到高潮| 成熟少妇高潮喷水视频| videosex国产| 999久久久精品免费观看国产| 免费看日本二区| 亚洲av日韩精品久久久久久密| 亚洲欧美精品综合一区二区三区| 性欧美人与动物交配| 十八禁网站免费在线| 亚洲专区国产一区二区| 欧美3d第一页|