• <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
    精品国产一区二区三区久久久樱花| 亚洲精品中文字幕在线视频| 女性被躁到高潮视频| 高清欧美精品videossex| 免费在线观看视频国产中文字幕亚洲 | 久久中文字幕一级| 男女之事视频高清在线观看 | 亚洲成国产人片在线观看| 亚洲av国产av综合av卡| 一边亲一边摸免费视频| 国产麻豆69| 中文字幕色久视频| 只有这里有精品99| 久久久久精品人妻al黑| 精品人妻熟女毛片av久久网站| 国产精品国产av在线观看| 十八禁高潮呻吟视频| 国产精品一二三区在线看| 麻豆乱淫一区二区| 777米奇影视久久| 精品久久久久久电影网| 国产一区二区三区av在线| 中国美女看黄片| 午夜福利在线免费观看网站| 亚洲人成网站在线观看播放| 午夜福利视频在线观看免费| 一边摸一边做爽爽视频免费| 一本久久精品| 丝袜美足系列| 高潮久久久久久久久久久不卡| 国产高清国产精品国产三级| 国产精品欧美亚洲77777| 桃花免费在线播放| 激情视频va一区二区三区| 精品高清国产在线一区| 咕卡用的链子| 少妇粗大呻吟视频| 女人久久www免费人成看片| 交换朋友夫妻互换小说| 桃花免费在线播放| 热re99久久国产66热| 国精品久久久久久国模美| 亚洲成人手机| 极品少妇高潮喷水抽搐| 久久精品人人爽人人爽视色| 欧美日韩视频精品一区| 欧美精品av麻豆av| 一区二区av电影网| 国产精品久久久久久精品电影小说| 脱女人内裤的视频| 亚洲伊人久久精品综合| bbb黄色大片| 又大又爽又粗| 国产精品久久久久久人妻精品电影 | 欧美亚洲日本最大视频资源| 激情视频va一区二区三区| 国产精品一二三区在线看| 国产亚洲一区二区精品| 久久国产精品男人的天堂亚洲| www.av在线官网国产| 国产亚洲精品久久久久5区| 首页视频小说图片口味搜索 | 国产91精品成人一区二区三区 | 在线亚洲精品国产二区图片欧美| 熟女av电影| 1024视频免费在线观看| 婷婷丁香在线五月| 午夜福利影视在线免费观看| 在线观看免费视频网站a站| 亚洲成色77777| 久9热在线精品视频| 日本猛色少妇xxxxx猛交久久| 99久久99久久久精品蜜桃| 欧美精品av麻豆av| 一二三四在线观看免费中文在| 亚洲av成人不卡在线观看播放网 | 三上悠亚av全集在线观看| 日韩视频在线欧美| 精品少妇黑人巨大在线播放| 中文字幕人妻丝袜一区二区| 19禁男女啪啪无遮挡网站| 亚洲人成网站在线观看播放| 在线观看免费视频网站a站| 欧美在线黄色| tube8黄色片| 下体分泌物呈黄色| 9191精品国产免费久久| 国产日韩一区二区三区精品不卡| 亚洲精品国产色婷婷电影| 自拍欧美九色日韩亚洲蝌蚪91| 日韩制服骚丝袜av| 精品久久久久久电影网| 国产精品成人在线| 免费高清在线观看视频在线观看| 成人手机av| 1024视频免费在线观看| 日韩人妻精品一区2区三区| 亚洲精品日韩在线中文字幕| 国产精品二区激情视频| www.自偷自拍.com| 亚洲精品成人av观看孕妇| 亚洲伊人久久精品综合| 黄色片一级片一级黄色片| 亚洲国产av新网站| 99久久精品国产亚洲精品| 国产一区二区在线观看av| 青草久久国产| 99香蕉大伊视频| 精品一区二区三区四区五区乱码 | 亚洲综合色网址| 777米奇影视久久| 欧美大码av| av国产精品久久久久影院| 国语对白做爰xxxⅹ性视频网站| 国产一区二区 视频在线| 久久九九热精品免费| 国产精品偷伦视频观看了| 人人妻人人澡人人看| 一级片'在线观看视频| av又黄又爽大尺度在线免费看| 在线观看免费高清a一片| 我要看黄色一级片免费的| 国产亚洲av高清不卡| av线在线观看网站| 免费少妇av软件| 黄色视频不卡| 国产成人精品久久久久久| 国产精品av久久久久免费| 婷婷色麻豆天堂久久| 国产一区二区在线观看av| 国产精品久久久久久精品古装| 亚洲国产精品成人久久小说| 另类精品久久| 无遮挡黄片免费观看| 亚洲欧美色中文字幕在线| 亚洲成av片中文字幕在线观看| www.av在线官网国产| 精品国产一区二区三区久久久樱花| 亚洲九九香蕉| 国产福利在线免费观看视频| 亚洲av国产av综合av卡| 你懂的网址亚洲精品在线观看| 亚洲人成77777在线视频| 中文字幕人妻丝袜制服| 国产福利在线免费观看视频| 嫁个100分男人电影在线观看 | 老汉色av国产亚洲站长工具| 女性生殖器流出的白浆| 亚洲激情五月婷婷啪啪| 日韩一区二区三区影片| 91精品国产国语对白视频| 久久热在线av| 亚洲,一卡二卡三卡| 大片免费播放器 马上看| 又黄又粗又硬又大视频| 日日爽夜夜爽网站| 好男人视频免费观看在线| 夜夜骑夜夜射夜夜干| 亚洲国产精品成人久久小说| 亚洲三区欧美一区| 国产精品一二三区在线看| 十八禁高潮呻吟视频| 一区二区日韩欧美中文字幕| a级片在线免费高清观看视频| 九草在线视频观看| 在线av久久热| 91九色精品人成在线观看| 亚洲图色成人| 欧美激情高清一区二区三区| 高清av免费在线| 免费日韩欧美在线观看| 人妻 亚洲 视频| 母亲3免费完整高清在线观看| 国产精品欧美亚洲77777| 亚洲欧美一区二区三区国产| 精品少妇黑人巨大在线播放| 日韩一卡2卡3卡4卡2021年| 亚洲国产欧美网| 国产真人三级小视频在线观看| 国产黄频视频在线观看| 热re99久久国产66热| 亚洲精品久久成人aⅴ小说| 久9热在线精品视频| 免费在线观看视频国产中文字幕亚洲 | 午夜视频精品福利| 大话2 男鬼变身卡| 亚洲精品av麻豆狂野| 欧美在线黄色| 精品国产一区二区久久| 久久国产精品影院| 丝袜美腿诱惑在线| 欧美久久黑人一区二区| 国产成人欧美| 亚洲中文av在线| 少妇人妻 视频| 精品少妇内射三级| 国产精品麻豆人妻色哟哟久久| 国产黄色免费在线视频| 国产99久久九九免费精品| 飞空精品影院首页| 免费高清在线观看日韩| 成人国语在线视频| 一本大道久久a久久精品| 免费在线观看完整版高清| 制服诱惑二区| 久久国产亚洲av麻豆专区| av天堂在线播放| 人人妻人人爽人人添夜夜欢视频| 欧美日韩黄片免| 欧美 日韩 精品 国产| 成人亚洲精品一区在线观看| 极品人妻少妇av视频| 亚洲自偷自拍图片 自拍| 国产日韩欧美在线精品| 婷婷色麻豆天堂久久| 黄片播放在线免费| 免费在线观看影片大全网站 | 首页视频小说图片口味搜索 | 狠狠婷婷综合久久久久久88av| 亚洲免费av在线视频| 女性被躁到高潮视频| 成人午夜精彩视频在线观看| 少妇的丰满在线观看| 久久久久网色| 宅男免费午夜| 波多野结衣一区麻豆| 脱女人内裤的视频| 丰满饥渴人妻一区二区三| 中文字幕另类日韩欧美亚洲嫩草| 亚洲综合色网址| 一级黄色大片毛片| 在线天堂中文资源库| 少妇人妻 视频| 最近手机中文字幕大全| 宅男免费午夜| 久久久精品国产亚洲av高清涩受| 中文字幕高清在线视频| 午夜免费观看性视频| 日韩大码丰满熟妇| 亚洲精品国产色婷婷电影| 黄网站色视频无遮挡免费观看| 大香蕉久久网| 免费久久久久久久精品成人欧美视频| 国产无遮挡羞羞视频在线观看| 久久中文字幕一级| 国产成人欧美在线观看 | av欧美777| 电影成人av| 亚洲,一卡二卡三卡| av天堂在线播放| 男女下面插进去视频免费观看| 多毛熟女@视频| 亚洲av综合色区一区| 波野结衣二区三区在线| 欧美另类一区| 亚洲精品国产av蜜桃| 亚洲黑人精品在线| 搡老乐熟女国产| 日韩中文字幕欧美一区二区 | 日本av免费视频播放| 老司机午夜十八禁免费视频| xxx大片免费视频| 亚洲国产精品一区三区| 一区二区三区四区激情视频| 日韩 亚洲 欧美在线| 涩涩av久久男人的天堂| av视频免费观看在线观看| av有码第一页| 一级毛片黄色毛片免费观看视频| 国产色视频综合| 国产成人欧美在线观看 | 黑丝袜美女国产一区| 成人手机av| 激情五月婷婷亚洲| 午夜免费鲁丝| 欧美精品高潮呻吟av久久| 久久天躁狠狠躁夜夜2o2o | 免费高清在线观看日韩| 亚洲综合色网址| 脱女人内裤的视频| 日韩 欧美 亚洲 中文字幕| 丝袜人妻中文字幕| 最近手机中文字幕大全| 国产不卡av网站在线观看| 婷婷色麻豆天堂久久| 欧美精品啪啪一区二区三区 | 日韩免费高清中文字幕av| 美女午夜性视频免费| 国产精品偷伦视频观看了| 我的亚洲天堂| 人妻一区二区av| 国产成人91sexporn| 手机成人av网站| 老汉色av国产亚洲站长工具| 亚洲国产欧美在线一区| 免费看av在线观看网站| 国产精品av久久久久免费| 少妇的丰满在线观看| a 毛片基地| 亚洲av电影在线进入| 18禁裸乳无遮挡动漫免费视频| 又紧又爽又黄一区二区| 久久精品人人爽人人爽视色| 国产人伦9x9x在线观看| 欧美av亚洲av综合av国产av| 免费人妻精品一区二区三区视频| 热re99久久国产66热| 一二三四社区在线视频社区8| 国产成人91sexporn| 国产视频一区二区在线看| 日韩熟女老妇一区二区性免费视频| 美女扒开内裤让男人捅视频| 18禁观看日本| 交换朋友夫妻互换小说| 精品国产国语对白av| 久久久欧美国产精品| 菩萨蛮人人尽说江南好唐韦庄| 久久久久国产精品人妻一区二区| 国产精品.久久久| av网站在线播放免费| 免费在线观看影片大全网站 | av网站在线播放免费| 亚洲精品中文字幕在线视频| 高清黄色对白视频在线免费看| 精品少妇久久久久久888优播| 在线观看人妻少妇| 最黄视频免费看| 丝袜脚勾引网站| 国产在线免费精品| √禁漫天堂资源中文www| 天天躁夜夜躁狠狠躁躁| 日韩av在线免费看完整版不卡| 午夜日韩欧美国产| 亚洲视频免费观看视频| 亚洲欧美成人综合另类久久久| netflix在线观看网站| 国产真人三级小视频在线观看| 看免费成人av毛片| 国产激情久久老熟女| 国产成人欧美| 交换朋友夫妻互换小说| 97人妻天天添夜夜摸| 成人亚洲精品一区在线观看| 欧美少妇被猛烈插入视频| 精品视频人人做人人爽| 精品久久久久久电影网| 老汉色∧v一级毛片| 久久久久久久久久久久大奶| 日日爽夜夜爽网站| 亚洲人成网站在线观看播放| 十分钟在线观看高清视频www| 日韩中文字幕欧美一区二区 | 国产男女超爽视频在线观看| 国产精品 欧美亚洲| 日本五十路高清| 亚洲国产av影院在线观看| 国产无遮挡羞羞视频在线观看| 高清黄色对白视频在线免费看| 久久国产亚洲av麻豆专区| 午夜福利一区二区在线看| 国产成人一区二区三区免费视频网站 | 久久久亚洲精品成人影院| 色播在线永久视频| 性高湖久久久久久久久免费观看| www.熟女人妻精品国产| 97在线人人人人妻| 香蕉国产在线看| 人妻人人澡人人爽人人| av在线老鸭窝| a 毛片基地| 黑人巨大精品欧美一区二区蜜桃| 午夜免费男女啪啪视频观看| 午夜福利在线免费观看网站| 成人手机av| av一本久久久久| 精品久久久精品久久久| 久久久久久久大尺度免费视频| 女警被强在线播放| 精品国产一区二区久久| 亚洲欧美日韩高清在线视频 | 一区二区三区乱码不卡18| 涩涩av久久男人的天堂| 欧美激情高清一区二区三区| av国产精品久久久久影院| 啦啦啦 在线观看视频| 久久精品久久精品一区二区三区| 九草在线视频观看| 亚洲欧美色中文字幕在线| 电影成人av| 国精品久久久久久国模美| 欧美日韩亚洲高清精品| 久久人妻熟女aⅴ| 国产精品 欧美亚洲| 亚洲 国产 在线| 久久人妻福利社区极品人妻图片 | 中文字幕亚洲精品专区| 国产有黄有色有爽视频| 精品一区在线观看国产| 韩国精品一区二区三区| 久久性视频一级片| 亚洲av综合色区一区| 日本午夜av视频| 18禁观看日本| 国产精品亚洲av一区麻豆| 免费少妇av软件| 欧美中文综合在线视频| a级片在线免费高清观看视频| 久久综合国产亚洲精品| 19禁男女啪啪无遮挡网站| 男男h啪啪无遮挡| 成在线人永久免费视频| 欧美 亚洲 国产 日韩一| 国产成人精品久久久久久| 欧美日韩av久久| a级毛片在线看网站| 亚洲欧美激情在线| 操出白浆在线播放| 免费黄频网站在线观看国产| 亚洲欧洲日产国产| a级片在线免费高清观看视频| 亚洲欧美日韩高清在线视频 | 成人免费观看视频高清| 免费看十八禁软件| 婷婷丁香在线五月| 亚洲av在线观看美女高潮| 久久九九热精品免费| 99国产精品一区二区三区| 精品国产一区二区三区久久久樱花| 国产免费视频播放在线视频| 国产一卡二卡三卡精品| 啦啦啦在线观看免费高清www| 极品少妇高潮喷水抽搐| 国产成人一区二区三区免费视频网站 | 91九色精品人成在线观看| 亚洲中文av在线| 久久精品久久精品一区二区三区| 91麻豆精品激情在线观看国产 | e午夜精品久久久久久久| 婷婷色av中文字幕| 亚洲国产欧美一区二区综合| 99国产精品99久久久久| 捣出白浆h1v1| 精品一区二区三区av网在线观看 | 男女国产视频网站| 久久国产精品男人的天堂亚洲| 久久这里只有精品19| 一本综合久久免费| 人人妻人人澡人人爽人人夜夜| 国产成人a∨麻豆精品| 久久精品国产亚洲av涩爱| 久久国产精品人妻蜜桃| 久久精品aⅴ一区二区三区四区| 两人在一起打扑克的视频| 建设人人有责人人尽责人人享有的| 肉色欧美久久久久久久蜜桃| 国产在视频线精品| 国产成人一区二区在线| 精品国产一区二区久久| 精品少妇一区二区三区视频日本电影| 亚洲av欧美aⅴ国产| 国产精品熟女久久久久浪| 九草在线视频观看| 色视频在线一区二区三区| a级片在线免费高清观看视频| 国产成人av教育| 亚洲三区欧美一区| 欧美日韩成人在线一区二区| 777久久人妻少妇嫩草av网站| 91精品三级在线观看| 亚洲欧洲精品一区二区精品久久久| 一区二区日韩欧美中文字幕| 狂野欧美激情性xxxx| 精品亚洲乱码少妇综合久久| 国产精品成人在线| 九色亚洲精品在线播放| 黄频高清免费视频| 丰满迷人的少妇在线观看| 美女午夜性视频免费| bbb黄色大片| 欧美精品啪啪一区二区三区 | 建设人人有责人人尽责人人享有的| 久久久国产精品麻豆| 制服人妻中文乱码| 久久精品aⅴ一区二区三区四区| 精品一区二区三区四区五区乱码 | 亚洲成人国产一区在线观看 | 免费少妇av软件| 精品一区二区三区四区五区乱码 | 国产在线视频一区二区| 久久热在线av| 精品人妻在线不人妻| 亚洲免费av在线视频| 国产爽快片一区二区三区| 欧美日韩视频精品一区| 一边摸一边做爽爽视频免费| 99国产精品一区二区蜜桃av | 18禁观看日本| 精品熟女少妇八av免费久了| 欧美成人午夜精品| 欧美日韩视频精品一区| 丝袜喷水一区| 婷婷丁香在线五月| 高清视频免费观看一区二区| 女警被强在线播放| 美女午夜性视频免费| 国产成人av激情在线播放| 老司机影院毛片| 五月天丁香电影| 欧美日韩视频精品一区| 在线观看人妻少妇| 两个人看的免费小视频| 成人亚洲欧美一区二区av| 亚洲专区国产一区二区| 精品国产乱码久久久久久小说| 国产在线观看jvid| 久久热在线av| 久久国产精品大桥未久av| 久久久久国产精品人妻一区二区| 日本一区二区免费在线视频| 国产淫语在线视频| 国产99久久九九免费精品| 成人亚洲精品一区在线观看| 老司机影院成人| 国产熟女午夜一区二区三区| 久久久国产欧美日韩av| 天堂8中文在线网| 美女脱内裤让男人舔精品视频| 国产视频一区二区在线看| 国产成人av教育| 午夜免费观看性视频| 久久久精品94久久精品| 国产不卡av网站在线观看| 亚洲欧洲精品一区二区精品久久久| 一区二区三区激情视频| 亚洲av成人不卡在线观看播放网 | 免费高清在线观看日韩| 久久精品久久精品一区二区三区| 久久久久精品人妻al黑| 黄色a级毛片大全视频| 午夜福利,免费看| 青春草亚洲视频在线观看| 国产午夜精品一二区理论片| 久久久国产一区二区| 国产真人三级小视频在线观看| 成年女人毛片免费观看观看9 | 国产高清视频在线播放一区 | 亚洲人成电影免费在线| 少妇粗大呻吟视频| bbb黄色大片| 日本欧美视频一区| 精品一区二区三区av网在线观看 | 国产亚洲精品第一综合不卡| svipshipincom国产片| 国产精品.久久久| 国产黄色视频一区二区在线观看| 亚洲av成人精品一二三区| 另类亚洲欧美激情| 91精品伊人久久大香线蕉| 亚洲一区二区三区欧美精品| 午夜视频精品福利| 日本猛色少妇xxxxx猛交久久| 国产在线免费精品| 久热爱精品视频在线9| 久久这里只有精品19| 国产不卡av网站在线观看| 大片免费播放器 马上看| 十八禁人妻一区二区| 韩国精品一区二区三区| 成人18禁高潮啪啪吃奶动态图| 欧美日韩综合久久久久久| 亚洲成人免费电影在线观看 | 人妻一区二区av| 欧美人与性动交α欧美软件| 国产精品一国产av| 成年动漫av网址| 亚洲,一卡二卡三卡| 日韩一本色道免费dvd| 亚洲国产精品国产精品| 青青草视频在线视频观看| 免费观看av网站的网址| 亚洲精品国产区一区二| 免费久久久久久久精品成人欧美视频| 女人久久www免费人成看片| 九色亚洲精品在线播放| 极品人妻少妇av视频| 欧美黑人欧美精品刺激| www日本在线高清视频| 国产成人av教育| 国产欧美日韩一区二区三 | 丰满饥渴人妻一区二区三| 国产成人欧美在线观看 | a级毛片黄视频| 久久精品久久久久久噜噜老黄| 国产精品国产av在线观看| 无限看片的www在线观看| 波多野结衣一区麻豆| 久热这里只有精品99| av国产久精品久网站免费入址| 国产精品一国产av| 久久精品aⅴ一区二区三区四区| 色播在线永久视频| 久久九九热精品免费| 日韩 欧美 亚洲 中文字幕| 午夜日韩欧美国产| 制服诱惑二区| 精品一品国产午夜福利视频| 欧美乱码精品一区二区三区| 大香蕉久久成人网| 久久国产精品男人的天堂亚洲| 热re99久久国产66热| 免费人妻精品一区二区三区视频| 男人添女人高潮全过程视频| 精品少妇内射三级| 日本av免费视频播放| 两个人免费观看高清视频|