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

    考慮精細模型的電-氣綜合能源系統(tǒng)優(yōu)化運行方法

    2022-11-23 11:17:58董文杰李俊格田志強孫英英曾順奇艾芊
    電測與儀表 2022年11期
    關(guān)鍵詞:潮流壓縮機約束

    董文杰,李俊格,田志強,孫英英,曾順奇,艾芊

    (1.東方電子股份有限公司, 山東 煙臺 264000; 2.廣東電網(wǎng)有限公司公司廣州供電局, 廣州 510620;3. 上海交通大學, 上海 200240)

    0 引 言

    近年來,隨著環(huán)境污染與能源危機不斷加劇,能源的高效清潔利用成為亟待解決的問題[1]。綜合能源系統(tǒng)能夠通過協(xié)調(diào)多種能源形式,實現(xiàn)能源的高效清潔利用[2-5]。為此,綜合能源系統(tǒng)的優(yōu)化運行成為當前能源領(lǐng)域的研究熱點。燃氣機組作為基本的電、氣耦合元件,并具有排放低、爬坡速率高等優(yōu)點[6],近年來其裝機容量大幅增加。燃氣機組的廣泛應(yīng)用使電-氣綜合能源系統(tǒng)的耦合程度更加緊密[7],也為電-氣綜合能源系統(tǒng)的聯(lián)合優(yōu)化運行帶來了挑戰(zhàn)。

    與電力系統(tǒng)相比,天然氣系統(tǒng)的動態(tài)過程緩慢,描述管道物理狀態(tài)的偏微分方程難以求解,文獻[8]在一定簡化條件下將偏微分方程轉(zhuǎn)化為代數(shù)方程;然而代數(shù)方程仍具有很強的非凸非線性特點,使得優(yōu)化問題的求解難度大、收斂性無法保證。為解決天然氣系統(tǒng)的優(yōu)化求解問題,該文獻提出了處理管道方程的分段線性化方法。在此基礎(chǔ)上,文獻[9-11]求解了電-氣綜合能源系統(tǒng)的聯(lián)合優(yōu)化問題。在電-氣綜合能源系統(tǒng)的聯(lián)合優(yōu)化領(lǐng)域,以往研究大多采用二次等式形式的Weymouth方程表征天然氣管道的物理特性,文獻[12]將Weymouth方程約束松弛為二階錐約束,進而將非凸優(yōu)化問題轉(zhuǎn)化為凸優(yōu)化問題進行求解。在此基礎(chǔ)上,文獻[13-14]采用迭代的方法進行修正以克服松弛解不可行的問題,文獻[15-18]考慮了天然氣管道流向未知的情況對二階錐松弛方法進行改進。

    在不同運行條件下,不同管道方程的擬合精度存在差異[19],例如Panhandle公式適用于高壓氣網(wǎng),而Spitzglass公式適用于低壓氣網(wǎng)。因此,Weymouth方程無法適用于所有天然氣管道,否則會導致誤差較大。然而,二階錐松弛方法僅能適用于二次等式形式的Weymouth方程;分段線性化方法難以權(quán)衡分段精度與求解代價,可能導致所求運行策略不可行的問題,無法滿足電-氣綜合能源系統(tǒng)的安全穩(wěn)定運行要求。此外,上述文獻中采用的模型還存在兩點不足:(1)壓縮機模型均為線性近似模型,認為壓縮機功率與流入流量成正比,實際上壓縮機耗氣量還與首末兩端壓強相關(guān);(2)電力系統(tǒng)模型均采用線性化的直流潮流模型或者忽略線路損耗的交流潮流模型。文獻[9-18]采用了精確的數(shù)學優(yōu)化算法,盡管可以得到所求模型的最優(yōu)解,但是所求模型為簡化模型,可能導致調(diào)度策略不可行。

    為了解決模型復雜度為優(yōu)化求解帶來的問題,考慮從智能算法的角度尋找解決方案,以得到精確模型下的可行解。粒子群優(yōu)化(Particle Swarm Optimization,PSO) 算法作為一種新興的群體智能優(yōu)化技術(shù),適用于目前科學領(lǐng)域、工程領(lǐng)域和經(jīng)濟領(lǐng)域中很多復雜的、非凸非線性形式的優(yōu)化問題。目前,PSO算法在電氣工程領(lǐng)域有著廣泛的應(yīng)用[20-21],涉及負荷經(jīng)濟分配[22]、無功優(yōu)化[23]、最優(yōu)潮流計算[24]、機組組合[25-26]、經(jīng)濟調(diào)度[27]等方面,并衍生出適用于特殊場景的改進算法。

    隨著綜合能源系統(tǒng)逐漸成為研究熱點,諸多研究將粒子群算法在電力系統(tǒng)中的應(yīng)用方法遷移至綜合能源系統(tǒng)領(lǐng)域[28]。文獻[29]建立了考慮用能成本、棄風損失、負荷曲線峰谷差的多目標優(yōu)化模型,研究了電-氣-熱綜合能源系統(tǒng)的魯棒優(yōu)化問題。針對傳統(tǒng)粒子群算法易陷入局部最優(yōu)的問題,文獻[30]將小生境技術(shù)與粒子群算法結(jié)合并提出一種非線性遞減慣性權(quán)重策略。文獻[31]建立了規(guī)劃和運行兩階段優(yōu)化模式。

    在粒子群算法的應(yīng)用中,為了考慮能源系統(tǒng)中的約束,已有文獻采用的方法大致分為兩種,一種是根據(jù)約束條件構(gòu)造罰項并添加至評價函數(shù)中[25],另一種是根據(jù)約束條件對結(jié)果進行檢驗,若不滿足則采用越限代限的方法重新求解,直至滿足約束條件[28]。然而,上述涉及粒子群算法應(yīng)用的文獻僅考慮了能量平衡和供能上下限約束,未考慮能源系統(tǒng)網(wǎng)絡(luò)約束(如潮流方程約束、管道方程約束等),可能導致調(diào)度策略不可行,出現(xiàn)電壓、電流、壓強越限等問題。

    在電力系統(tǒng)中,根據(jù)機組出力等決策變量可以通過潮流計算得到如線路潮流、母線電壓等狀態(tài)變量,從而可以考慮潮流方程約束以及狀態(tài)變量的上下限約束。在天然氣系統(tǒng)中,根據(jù)氣井產(chǎn)氣量等決策變量也可以通過氣網(wǎng)能量流計算得到如管道流量、節(jié)點壓強等狀態(tài)變量,從而可以考慮氣網(wǎng)管道方程約束以及節(jié)點壓強上下限約束。對于電力系統(tǒng)與天然氣系統(tǒng)通過燃氣輪機、電壓縮機雙向耦合的場景,文獻[32-33]給出一種基于分別循環(huán)迭代的求解算法,通過迭代求解氣網(wǎng)與電網(wǎng)得到綜合能源系統(tǒng)的能量流分布。

    在已有研究成果的基礎(chǔ)上,文章旨在提出一種考慮電-氣綜合能源系統(tǒng)精細模型的粒子群優(yōu)化方法,將電-氣綜合能源系統(tǒng)能量流計算嵌入到粒子群優(yōu)化算法中,從而能夠考慮較為精細的電-氣綜合能源系統(tǒng)模型。文中首先根據(jù)粒子位置得到機組出力、氣井產(chǎn)氣量,通過分別循環(huán)迭代方法迭代求解電網(wǎng)、氣網(wǎng)潮流,迭代收斂后即可得到電氣綜合能源系統(tǒng)的運行狀態(tài),進而將運行上下限約束作為罰項添加到粒子的評價函數(shù)中,并賦以較高的權(quán)重,從而得到電-氣綜合能源系統(tǒng)的協(xié)同優(yōu)化運行策略。

    1 電-氣綜合能源系統(tǒng)協(xié)同優(yōu)化運行模型

    研究對象為電-氣綜合能源系統(tǒng),其包含兩個子系統(tǒng),即電力網(wǎng)絡(luò)和天然氣網(wǎng)絡(luò),電力網(wǎng)絡(luò)與天然氣網(wǎng)絡(luò)通過燃氣輪機、電壓縮機等裝置耦合。由于能源系統(tǒng)的優(yōu)化不涉及動態(tài)過程,采用精細的配電網(wǎng)交流潮流模型以及天然氣網(wǎng)的穩(wěn)態(tài)模型。電-氣綜合能源系統(tǒng)示意圖如圖1所示,文中考慮電力網(wǎng)絡(luò)與天然氣網(wǎng)絡(luò)雙向耦合的情況。

    圖1 電-氣綜合能源系統(tǒng)結(jié)構(gòu)

    1.1 電力系統(tǒng)約束條件

    (1)常規(guī)機組和燃氣機組功率約束:

    (1)

    (2)

    (3)

    (4)

    (2)常規(guī)機組和燃氣機組爬坡約束:

    (5)

    (6)

    (3)母線電壓約束:

    (7)

    (4)線路容量約束:

    (8)

    (5)線路電流約束:

    (9)

    該式表示線路電流與傳輸有功、無功功率、首端母線電壓的關(guān)系。式中,plt為線路流過有功功率;qlt為線路流過無功功率;vl+t為線路首端電壓;Ilt為支路電流的平方。

    (6)潮流方程約束:

    (10)

    (11)

    (12)

    1.2 天然氣系統(tǒng)約束條件

    (1)天然氣氣井產(chǎn)氣量約束

    (13)

    (2)天然氣網(wǎng)絡(luò)節(jié)點壓強約束

    (14)

    (3)天然氣管道方程約束

    (15)

    考慮到氣網(wǎng)的運行壓力通常較高,因此采用了適用于0.7 MPa以上壓力運行的高壓氣網(wǎng)Panhandle’A’公式。式中,fpt為管道流過的流量;τp+t/τp-t為管道首端、末端壓強;φp為管道常數(shù)。對于其他形式的管道方程約束,所提方法依然適用。

    (4)天然氣網(wǎng)絡(luò)的節(jié)點流量平衡約束:

    (16)

    式中(n)為與節(jié)點n相連的氣井集合; (n+)/(n-)為以節(jié)點n為首端、末端節(jié)點的管道集合;(n+)/(n-)為以節(jié)點n為首端、末端節(jié)點的壓縮機集合;g(n)為與節(jié)點n相連的天然氣負荷集合;g(n)為與節(jié)點n相連的燃氣輪機集合;Fdgt為天然氣負荷流量;fugt為燃氣機組耗氣量;fct為流過壓縮機的流量。

    采用了精細化的壓縮機模型,電壓縮機約束如下:

    (17)

    (18)

    式中τc-t為入口壓力;τc+t為出口壓力;α為多變指數(shù)并取1.27;ηc為壓縮機的效率;λc為電壓縮機的功率因數(shù)。考慮壓縮機處于定出口壓強工作狀態(tài),出口壓強作為系統(tǒng)決策變量。

    1.3 耦合約束條件

    燃氣輪機耗氣量約束:

    pugt=ηugKfugt,?ug,t

    (19)

    式中ηug為燃氣輪機的效率;K為天然氣熱值。

    1.4 電-氣綜合能源系統(tǒng)的目標函數(shù)

    電-氣綜合能源系統(tǒng)優(yōu)化模型的目標函數(shù)即為一個調(diào)度周期內(nèi)的運行成本,包括常規(guī)機組成本與氣井產(chǎn)氣成本。

    (20)

    2 粒子群算法及其應(yīng)用

    2.1 粒子群算法的應(yīng)用方法

    根據(jù)粒子群算法的求解流程可知,粒子群算法求解的是一個無約束優(yōu)化問題,盡管可以限制粒子的位置,但僅能表示簡單的決策變量上下限約束。例如,電網(wǎng)優(yōu)化調(diào)度中機組出力為決策變量,機組出力上下限約束可以通過限制粒子在空間中的位置實現(xiàn)。但是,電網(wǎng)優(yōu)化模型中還含有其他復雜約束,如潮流方程約束、機組爬坡約束、線路電流約束等,對于電-氣綜合能源系統(tǒng)優(yōu)化問題亦是如此。

    對于電-氣綜合能源系統(tǒng),其決策變量即為機組有功無功出力、氣井產(chǎn)氣量、壓縮機出口壓強,當決策變量通過粒子位置確定時,耦合系統(tǒng)的運行工況即可被唯一確定,從而將各類約束作為罰項添加至粒子的評價函數(shù)中。根據(jù)現(xiàn)有的潮流計算方法,上述求解思路存在如下問題:

    (1)從物理意義角度,機組有功無功出力、氣井產(chǎn)氣量全部由粒子位置確定,難以保證電網(wǎng)、氣網(wǎng)的能量流平衡,繼而無法獲得耦合系統(tǒng)能量流分布;

    (2)從計算方法角度,現(xiàn)有的電力系統(tǒng)潮流計算方法需要設(shè)定平衡節(jié)點,以平衡線路上的有功無功損失;氣網(wǎng)潮流計算方法也需要設(shè)定平衡氣井,保證全網(wǎng)的平衡。

    針對上述問題,將電力系統(tǒng)中燃氣機組所在節(jié)點作為平衡節(jié)點,其節(jié)點電壓、相角作為決策變量由粒子位置確定;常規(guī)機組節(jié)點作為PV節(jié)點,其有功功率、節(jié)點電壓作為決策變量由粒子位置確定。經(jīng)過潮流計算,可以得到燃氣機組的有功無功功率、常規(guī)機組的無功功率,以及全網(wǎng)的潮流分布。對于氣網(wǎng),設(shè)定某一氣井作為平衡氣井,其節(jié)點壓強作為決策變量,其余氣井的產(chǎn)氣量、壓縮機出口壓強作為決策變量。經(jīng)過潮流計算,可以得到平衡氣井的產(chǎn)氣量以及氣網(wǎng)的能量流分布。電網(wǎng)與氣網(wǎng)通過燃氣輪機、電壓縮機耦合,通過分別循環(huán)迭代算法迭代求解電網(wǎng)、氣網(wǎng)能量流即可得到耦合系統(tǒng)的潮流分布[33],求解方法見2.2節(jié)所示。

    因此,通過給定電網(wǎng)、氣網(wǎng)的決策變量,即可得到耦合系統(tǒng)的能量流分布,結(jié)果必然滿足潮流方程約束與氣網(wǎng)中的管道方程、節(jié)點流量平衡約束。由于平衡機組的出力、平衡氣井的產(chǎn)氣量無法通過粒子位置進行限制,只能通過添加罰項的方式進行處理。此外,爬坡約束、母線電壓約束、線路容量約束、節(jié)點壓強約束均需要通過添加罰項的方法進行處理。添加罰項后的目標函數(shù)如下:

    (21)

    式中第一項表示常規(guī)機組的發(fā)電成本和氣井的產(chǎn)氣成本;第二項表示正爬約束對應(yīng)的罰項;第三項表示負爬坡約束對應(yīng)的罰項;第四項表示線路功率約束對應(yīng)的罰項;第五項表示全網(wǎng)功率平衡約束對應(yīng)的罰項;第六項為氣井的產(chǎn)氣成本;第七項為節(jié)點壓強約束對應(yīng)的罰項,M1、M2、M3、M4、M5、M6為對應(yīng)的罰因子。

    2.2 電-氣綜合能源系統(tǒng)潮流計算方法

    對于電-氣耦合系統(tǒng),考慮燃氣輪機、電壓縮機為耦合元件的情況,在根據(jù)決策變量求解耦合系統(tǒng)能量流時,為提高求解速度采用一種分別循環(huán)迭代的求解方法[33]首先假設(shè)電壓縮機功率初值,然后由機組出力確定電網(wǎng)潮流分布,根據(jù)式(19)確定燃氣機組的耗氣量,求解氣網(wǎng)潮流,根據(jù)式(17)~式(18)確定電壓縮機功率作為電網(wǎng)負荷重新進行電力系統(tǒng)潮流計算,直至燃氣輪機耗氣量及電壓縮機功率不再變化為止,計算流程如圖2所示。其中,氣網(wǎng)潮流計算采用牛頓拉夫遜法,考慮到多臺燃氣輪機同時作為平衡機組,電網(wǎng)潮流計算采用基于注入電流的多平衡節(jié)點潮流算法。

    圖2 電-氣耦合系統(tǒng)能量流求解流程圖

    3 算例分析

    為了驗證優(yōu)化程序的正確性,選取IEEE 13節(jié)點電力系統(tǒng)和7節(jié)點天然氣系統(tǒng)通過燃氣輪機、電壓縮機耦合得到如圖3所示電-氣綜合能源系統(tǒng)作為測試系統(tǒng),測試系統(tǒng)及負荷參數(shù)來自文獻[34]。

    圖3 測試算例拓撲圖

    3.1 仿真結(jié)果分析

    對上述系統(tǒng),取粒子種群數(shù)為100,迭代次數(shù)為200,罰系數(shù)取為20 000,得到的粒子群算法迭代曲線如圖4所示,圖4(a)表示不含罰項的目標函數(shù)迭代曲線,圖4(b)表示含有罰項的目標函數(shù)迭代曲線,即為粒子群算法的評價函數(shù)。由圖可知,含罰項的目標函數(shù)值隨迭代次數(shù)增加而逐漸減小,并最終達到5.1324×106,不含罰項的目標函數(shù)值隨迭代次數(shù)的變化并無明顯規(guī)律,最終值5.12875×106與初值較為接近,但解的可行性得到改善。此外,迭代200次時罰項仍不為零,說明還存在部分約束越限,隨著迭代次數(shù)的增大罰項5.0438×106,此時罰項較小,得到可行性更好的解。

    圖4 粒子群算法的迭代曲線

    粒子群算法迭代1 000次后得到的機組最優(yōu)出力和氣井最優(yōu)產(chǎn)氣量曲線如圖5與圖6所示。

    圖5 機組最優(yōu)出力曲線

    圖6 氣井最優(yōu)產(chǎn)氣曲線

    罰項的大小與罰系數(shù)、迭代次數(shù)、種群數(shù)量有關(guān),增加迭代次數(shù)、種群數(shù)量也會導致計算代價的增大。在固定迭代次數(shù)為200、種群數(shù)量為100的前提下,設(shè)置不同的罰系數(shù)進行測試,測試結(jié)果如表1所示,罰項與罰系數(shù)的比值描述了約束條件的越限量。

    表1 不同罰系數(shù)下的對比測試

    由表1可知,當罰系數(shù)的數(shù)量級增大時,罰項與罰系數(shù)的比值先減小后出現(xiàn)增大的趨勢,當罰系數(shù)較小時,增加罰系數(shù)的數(shù)量級有利于降低約束條件的越限量,當罰系數(shù)過大時反而引起越限量的增大。在同一數(shù)量級下,由于粒子群算法的求解結(jié)果具有一定的隨機性,使得罰系數(shù)與越限量并無明顯的負相關(guān)性。

    為分析評價函數(shù)中各約束越限量所占比重,在罰系數(shù)取20 000、種群數(shù)取100、迭代次數(shù)取200的情況下,得到各約束越限量所占比重如圖7中餅狀圖所示,圖中未展示的約束均未產(chǎn)生越限??梢钥吹郊s束越限量主要集中于機組有功約束和爬坡約束中。受限于粒子群算法難以找到全局最優(yōu)解,約束越限量無法完全消除,只能通過增加罰系數(shù)盡可能減少約束越限量。

    圖7 各約束越限量所占比例

    3.2 模型對比測試分析

    諸多關(guān)于粒子群算法在能源系統(tǒng)優(yōu)化中應(yīng)用的研究工作采用不考慮網(wǎng)絡(luò)拓撲的簡單模型[35-36],也有研究工作考慮了線性化的電網(wǎng)模型,及直流潮流模型。這些模型雖然具有易于求解的優(yōu)點,但由于簡化可能導致求解結(jié)果過于樂觀。為了體現(xiàn)考慮復雜網(wǎng)絡(luò)模型的必要性,在相同的算例參數(shù)下,采用三種模型得到的優(yōu)化結(jié)果如表2所示。其中,模型1考慮了電網(wǎng)的交流潮流模型以及氣網(wǎng)穩(wěn)態(tài)模型,模型2考慮了電網(wǎng)和直流潮流模型及氣網(wǎng)穩(wěn)態(tài)模型,模型3忽略了電網(wǎng)與氣網(wǎng)的網(wǎng)絡(luò)拓撲約束,僅考慮其中的能量平衡約束。

    表2 三種模型優(yōu)化結(jié)果對比

    從表2中結(jié)果可知,模型1的目標函數(shù)值最大,其他模型的結(jié)果均過于樂觀。主要原因在于模型3中忽略了網(wǎng)絡(luò)拓撲約束,無法考慮線路、管道容量限制,模型2采用的直流潮流模型忽略了線路電阻及其產(chǎn)生的功率損耗。

    與現(xiàn)有電-氣綜合能源系統(tǒng)優(yōu)化相關(guān)文獻相比[37-38],文中采用了精細化的壓縮機模型,而非假設(shè)壓縮消耗功率與流入流量成正比,得到的壓縮機有功功率、流入流量曲線如圖8所示。由圖8可知,壓縮機有功功率與流入流量并非完全成正比,在某些時段功率變化趨勢甚至與流量變化趨勢相反。因此,考慮精細化的壓縮機建模方法具有一定的意義。

    4 結(jié)束語

    提出了一種基于粒子群算法的電-氣綜合能源系統(tǒng)的優(yōu)化方法,并考慮了精細化的電網(wǎng)交流潮流約束和氣網(wǎng)穩(wěn)態(tài)約束模型。為處理電網(wǎng)與氣網(wǎng)中的網(wǎng)絡(luò)拓撲約束,采用基于分別循環(huán)迭代的電-氣綜合能源系統(tǒng)潮流計算方法,根據(jù)耦合系統(tǒng)決策變量計算得到全網(wǎng)的能量流分布,繼而將上下限約束作為罰項添加到目標函數(shù)中。最后,通過算例測試驗證了所提求解方法的有效性,分析了罰系數(shù)對約束越限量的影響,并驗證了考慮精細化模型的必要性。

    猜你喜歡
    潮流壓縮機約束
    “碳中和”約束下的路徑選擇
    約束離散KP方程族的完全Virasoro對稱
    BOG壓縮機在小型LNG船舶上的應(yīng)用
    淺析特大型往復壓縮機組的安裝
    壓縮機組智能控制與節(jié)能增效
    潮流
    足球周刊(2016年14期)2016-11-02 11:47:59
    潮流
    足球周刊(2016年15期)2016-11-02 11:44:02
    潮流
    足球周刊(2016年10期)2016-10-08 18:50:29
    從2014到2015潮流就是“貪新厭舊”
    Coco薇(2015年1期)2015-08-13 21:35:10
    適當放手能讓孩子更好地自我約束
    人生十六七(2015年6期)2015-02-28 13:08:38
    亚洲精品色激情综合| 欧美最新免费一区二区三区| 麻豆久久精品国产亚洲av| 亚洲美女视频黄频| 91麻豆精品激情在线观看国产| 日本五十路高清| 毛片一级片免费看久久久久 | 在线观看一区二区三区| 日本 av在线| 香蕉av资源在线| 亚洲av.av天堂| 亚洲av免费在线观看| 亚洲av成人av| 国产成人福利小说| 女生性感内裤真人,穿戴方法视频| 中文字幕熟女人妻在线| 深爱激情五月婷婷| 欧美xxxx黑人xx丫x性爽| 色视频www国产| a级一级毛片免费在线观看| 日本-黄色视频高清免费观看| 精品人妻视频免费看| 琪琪午夜伦伦电影理论片6080| 嫁个100分男人电影在线观看| 校园春色视频在线观看| 色综合色国产| 国产主播在线观看一区二区| 色精品久久人妻99蜜桃| 又爽又黄a免费视频| 色综合色国产| 久久精品夜夜夜夜夜久久蜜豆| 欧美日韩乱码在线| 国产v大片淫在线免费观看| 看黄色毛片网站| 悠悠久久av| 在线免费观看的www视频| 99国产精品一区二区蜜桃av| 亚洲不卡免费看| 亚洲va日本ⅴa欧美va伊人久久| 联通29元200g的流量卡| 全区人妻精品视频| 深夜精品福利| 精品欧美国产一区二区三| 亚洲欧美精品综合久久99| 赤兔流量卡办理| 亚洲成人久久性| 亚洲av免费在线观看| 在线观看午夜福利视频| 九色国产91popny在线| 国产高清不卡午夜福利| 一进一出好大好爽视频| 国产亚洲av嫩草精品影院| 日本在线视频免费播放| 久久香蕉精品热| 久久99热这里只有精品18| 69av精品久久久久久| 免费看日本二区| 久久久久免费精品人妻一区二区| 国产精品无大码| 超碰av人人做人人爽久久| 亚洲欧美激情综合另类| 内射极品少妇av片p| 国产男靠女视频免费网站| 久久久久久国产a免费观看| 欧美一区二区亚洲| 国产男人的电影天堂91| 国产欧美日韩一区二区精品| 麻豆一二三区av精品| 22中文网久久字幕| 国产精品日韩av在线免费观看| 麻豆国产av国片精品| 亚洲av电影不卡..在线观看| av黄色大香蕉| 成年女人看的毛片在线观看| 亚洲av免费在线观看| videossex国产| 男插女下体视频免费在线播放| 日韩欧美在线二视频| 国产大屁股一区二区在线视频| 欧美激情国产日韩精品一区| 日本色播在线视频| 两个人视频免费观看高清| 老司机午夜福利在线观看视频| 国产伦在线观看视频一区| 国产av在哪里看| 乱人视频在线观看| 欧美成人免费av一区二区三区| 欧美人与善性xxx| 国产午夜精品久久久久久一区二区三区 | 久久久久久久亚洲中文字幕| 国内少妇人妻偷人精品xxx网站| 亚洲精品粉嫩美女一区| 在线免费观看不下载黄p国产 | 日本免费a在线| 日本黄色片子视频| 老熟妇乱子伦视频在线观看| 国产91精品成人一区二区三区| 久久久久国内视频| 国产探花在线观看一区二区| 亚洲 国产 在线| 亚洲av日韩精品久久久久久密| 午夜激情欧美在线| 少妇熟女aⅴ在线视频| 精华霜和精华液先用哪个| av视频在线观看入口| 欧美另类亚洲清纯唯美| 一区二区三区激情视频| 久久精品综合一区二区三区| 又爽又黄无遮挡网站| 日韩欧美在线乱码| 亚洲美女黄片视频| 午夜亚洲福利在线播放| 精品福利观看| 99久久精品热视频| 久久久色成人| 真人做人爱边吃奶动态| 97人妻精品一区二区三区麻豆| 69人妻影院| 内地一区二区视频在线| 女同久久另类99精品国产91| 国产免费一级a男人的天堂| 在线天堂最新版资源| 欧美一区二区国产精品久久精品| 国产高清有码在线观看视频| 国产亚洲精品久久久com| 91麻豆av在线| av专区在线播放| 久久精品国产亚洲av涩爱 | 欧美日韩乱码在线| 久久久久久伊人网av| 成年女人看的毛片在线观看| 人妻制服诱惑在线中文字幕| 亚洲欧美精品综合久久99| 真人一进一出gif抽搐免费| 久久亚洲真实| 搡老妇女老女人老熟妇| 欧美一区二区国产精品久久精品| 美女黄网站色视频| www日本黄色视频网| 欧美又色又爽又黄视频| 日本-黄色视频高清免费观看| 亚洲av五月六月丁香网| 日日摸夜夜添夜夜添av毛片 | 日韩精品青青久久久久久| 少妇被粗大猛烈的视频| 一级毛片久久久久久久久女| 久久九九热精品免费| 欧美色欧美亚洲另类二区| 久久久久精品国产欧美久久久| 色综合婷婷激情| 看十八女毛片水多多多| 久久欧美精品欧美久久欧美| 俄罗斯特黄特色一大片| 一本精品99久久精品77| 人妻夜夜爽99麻豆av| 综合色av麻豆| 久久草成人影院| 波野结衣二区三区在线| 免费av毛片视频| 亚洲欧美日韩高清专用| 久久精品国产亚洲av涩爱 | 国产精品久久久久久亚洲av鲁大| 日本 欧美在线| 一本精品99久久精品77| 人妻夜夜爽99麻豆av| 国产一区二区三区在线臀色熟女| 国产单亲对白刺激| 永久网站在线| 午夜亚洲福利在线播放| 久久久精品大字幕| 18禁黄网站禁片免费观看直播| 最近最新免费中文字幕在线| 国产精品人妻久久久久久| 日本色播在线视频| 亚洲国产色片| 黄色一级大片看看| 麻豆国产97在线/欧美| 久久久久性生活片| 女同久久另类99精品国产91| 深爱激情五月婷婷| 国产成人a区在线观看| 午夜福利高清视频| 欧美性猛交╳xxx乱大交人| 亚洲avbb在线观看| 亚洲精品亚洲一区二区| 99riav亚洲国产免费| 窝窝影院91人妻| 最好的美女福利视频网| 亚洲av免费高清在线观看| 亚洲午夜理论影院| 国产精品国产三级国产av玫瑰| 18禁黄网站禁片午夜丰满| 特级一级黄色大片| 日韩欧美一区二区三区在线观看| 非洲黑人性xxxx精品又粗又长| 内射极品少妇av片p| 香蕉av资源在线| 乱系列少妇在线播放| 欧美绝顶高潮抽搐喷水| 亚洲熟妇熟女久久| 亚洲第一电影网av| av中文乱码字幕在线| 亚洲aⅴ乱码一区二区在线播放| 九色国产91popny在线| 亚洲精品亚洲一区二区| 一本久久中文字幕| 亚洲最大成人中文| 丰满乱子伦码专区| 国产大屁股一区二区在线视频| 男女啪啪激烈高潮av片| 亚洲va在线va天堂va国产| 免费一级毛片在线播放高清视频| 久久这里只有精品中国| 日本黄色片子视频| 悠悠久久av| 搡老妇女老女人老熟妇| 国产激情偷乱视频一区二区| 久久99热这里只有精品18| 国产成人aa在线观看| 看免费成人av毛片| 国产一区二区激情短视频| 91麻豆av在线| 波多野结衣高清无吗| 欧美在线一区亚洲| 中文资源天堂在线| 天堂动漫精品| 亚洲不卡免费看| 中文字幕av成人在线电影| 日本a在线网址| 亚洲人成网站高清观看| 国产av在哪里看| 18禁裸乳无遮挡免费网站照片| 国产精品野战在线观看| 人人妻,人人澡人人爽秒播| 亚洲熟妇熟女久久| 精品无人区乱码1区二区| 99热精品在线国产| av.在线天堂| 国产精品久久久久久精品电影| 99国产精品一区二区蜜桃av| 日韩欧美在线二视频| 又黄又爽又免费观看的视频| 久久99热这里只有精品18| 99国产精品一区二区蜜桃av| 国产av一区在线观看免费| 国产毛片a区久久久久| 深夜a级毛片| 国产精品自产拍在线观看55亚洲| 91麻豆av在线| 99视频精品全部免费 在线| 最新中文字幕久久久久| 精品久久久久久久末码| 淫妇啪啪啪对白视频| 1000部很黄的大片| 午夜老司机福利剧场| 国产伦在线观看视频一区| 亚洲自偷自拍三级| 久久久久久久久久久丰满 | 亚洲国产色片| 国产精品久久久久久久电影| 夜夜爽天天搞| 免费看光身美女| 国产国拍精品亚洲av在线观看| 精品久久久久久成人av| 亚洲美女黄片视频| 无人区码免费观看不卡| 国产女主播在线喷水免费视频网站 | 我要搜黄色片| 欧美高清成人免费视频www| 国产精品免费一区二区三区在线| 亚洲av电影不卡..在线观看| 深夜a级毛片| 在线免费十八禁| 欧美不卡视频在线免费观看| 天美传媒精品一区二区| 在线观看舔阴道视频| 亚洲国产欧洲综合997久久,| 午夜久久久久精精品| 亚洲成人精品中文字幕电影| av女优亚洲男人天堂| 日韩欧美在线乱码| 色哟哟哟哟哟哟| 国产伦精品一区二区三区四那| 看十八女毛片水多多多| www.色视频.com| 一边摸一边抽搐一进一小说| 91精品国产九色| 国产精品乱码一区二三区的特点| 色噜噜av男人的天堂激情| 日韩高清综合在线| 欧美最黄视频在线播放免费| 大型黄色视频在线免费观看| 亚洲av第一区精品v没综合| 国产综合懂色| 亚洲真实伦在线观看| 搡老妇女老女人老熟妇| 国产色婷婷99| 国产精品亚洲美女久久久| 亚洲中文字幕一区二区三区有码在线看| 九九热线精品视视频播放| 校园春色视频在线观看| 亚洲色图av天堂| 国产黄a三级三级三级人| 久久草成人影院| 国产午夜福利久久久久久| 他把我摸到了高潮在线观看| 一区二区三区免费毛片| 人妻丰满熟妇av一区二区三区| 97人妻精品一区二区三区麻豆| 99久久无色码亚洲精品果冻| 日本撒尿小便嘘嘘汇集6| 国内精品一区二区在线观看| 欧美精品国产亚洲| 久久精品国产亚洲网站| 少妇高潮的动态图| 国产高清视频在线观看网站| 麻豆av噜噜一区二区三区| 一个人看视频在线观看www免费| 欧美精品国产亚洲| 中文字幕免费在线视频6| 亚洲va日本ⅴa欧美va伊人久久| 色噜噜av男人的天堂激情| 久99久视频精品免费| 黄色一级大片看看| 麻豆国产97在线/欧美| 亚洲人与动物交配视频| 中文字幕高清在线视频| 一区二区三区激情视频| 搡女人真爽免费视频火全软件 | 国产视频内射| 高清在线国产一区| 久久久久久久久久成人| 两个人的视频大全免费| eeuss影院久久| 一区福利在线观看| 男女边吃奶边做爰视频| 亚洲四区av| 亚洲av成人精品一区久久| 亚洲精品一卡2卡三卡4卡5卡| 亚洲中文字幕一区二区三区有码在线看| 欧美又色又爽又黄视频| 丰满人妻一区二区三区视频av| 乱人视频在线观看| 成年女人毛片免费观看观看9| 永久网站在线| 亚洲av一区综合| 精品99又大又爽又粗少妇毛片 | 国产精华一区二区三区| 国产亚洲精品综合一区在线观看| 少妇被粗大猛烈的视频| 国产精华一区二区三区| 免费在线观看日本一区| 夜夜看夜夜爽夜夜摸| 九九热线精品视视频播放| 欧美成人免费av一区二区三区| av福利片在线观看| 国产一区二区在线av高清观看| 性欧美人与动物交配| 久久香蕉精品热| 成人特级av手机在线观看| 亚洲avbb在线观看| 不卡一级毛片| 偷拍熟女少妇极品色| 又爽又黄a免费视频| 日韩一本色道免费dvd| 69av精品久久久久久| 大又大粗又爽又黄少妇毛片口| 人人妻人人澡欧美一区二区| 欧美xxxx黑人xx丫x性爽| 成人美女网站在线观看视频| 色在线成人网| 1024手机看黄色片| 成人国产综合亚洲| 欧美一区二区亚洲| 日本熟妇午夜| 少妇猛男粗大的猛烈进出视频 | av专区在线播放| 少妇的逼好多水| 成人美女网站在线观看视频| 九色成人免费人妻av| 大又大粗又爽又黄少妇毛片口| 88av欧美| 国产精品野战在线观看| 草草在线视频免费看| 精品一区二区三区视频在线| 日韩大尺度精品在线看网址| 欧美又色又爽又黄视频| 嫩草影院入口| 日韩欧美在线乱码| 国产主播在线观看一区二区| 乱系列少妇在线播放| 国产探花极品一区二区| 99热只有精品国产| 美女高潮的动态| 久久香蕉精品热| 偷拍熟女少妇极品色| 男女那种视频在线观看| 无遮挡黄片免费观看| 99riav亚洲国产免费| av在线亚洲专区| 免费电影在线观看免费观看| 亚洲成a人片在线一区二区| 全区人妻精品视频| 日本免费a在线| 国产又黄又爽又无遮挡在线| 嫩草影视91久久| 亚洲精品久久国产高清桃花| 啦啦啦观看免费观看视频高清| 日本 欧美在线| 国内精品久久久久精免费| 久久久久久久久大av| 欧美3d第一页| 在线天堂最新版资源| 欧美性感艳星| 免费一级毛片在线播放高清视频| 91麻豆精品激情在线观看国产| 特大巨黑吊av在线直播| 国产乱人视频| 久久天躁狠狠躁夜夜2o2o| 中文字幕高清在线视频| 哪里可以看免费的av片| 性欧美人与动物交配| 国产高清不卡午夜福利| 噜噜噜噜噜久久久久久91| 久久国产精品人妻蜜桃| 免费看av在线观看网站| 久久久精品大字幕| 18禁黄网站禁片免费观看直播| 日本欧美国产在线视频| 亚洲无线观看免费| 免费电影在线观看免费观看| 成人鲁丝片一二三区免费| 麻豆久久精品国产亚洲av| 免费人成视频x8x8入口观看| 免费av观看视频| 色综合亚洲欧美另类图片| 97超级碰碰碰精品色视频在线观看| 一级黄色大片毛片| 一级毛片久久久久久久久女| 一区二区三区高清视频在线| 日日干狠狠操夜夜爽| 黄色欧美视频在线观看| avwww免费| 夜夜夜夜夜久久久久| 国产精品综合久久久久久久免费| videossex国产| 成人永久免费在线观看视频| 亚洲自偷自拍三级| 国产精品电影一区二区三区| 国产一区二区三区在线臀色熟女| 亚洲国产精品合色在线| 男人的好看免费观看在线视频| 亚洲国产高清在线一区二区三| 久9热在线精品视频| 国内少妇人妻偷人精品xxx网站| 久久久久性生活片| 最新中文字幕久久久久| 久久久精品大字幕| av在线天堂中文字幕| 日韩在线高清观看一区二区三区 | 午夜福利视频1000在线观看| 精品人妻偷拍中文字幕| 久久精品国产99精品国产亚洲性色| 少妇丰满av| 乱码一卡2卡4卡精品| 少妇熟女aⅴ在线视频| 日韩欧美三级三区| 中文字幕高清在线视频| 久久久久久久精品吃奶| av.在线天堂| 国产麻豆成人av免费视频| 深夜a级毛片| 亚洲第一电影网av| 色哟哟·www| 国产成人一区二区在线| 国产精品久久视频播放| 久久久久久九九精品二区国产| 久久99热6这里只有精品| 久久精品久久久久久噜噜老黄 | 亚洲国产日韩欧美精品在线观看| 成人国产综合亚洲| 听说在线观看完整版免费高清| 99热这里只有是精品50| 男人舔女人下体高潮全视频| 少妇人妻精品综合一区二区 | 国产男人的电影天堂91| 日韩欧美国产在线观看| 在线观看一区二区三区| 成年免费大片在线观看| 精品人妻偷拍中文字幕| 精品一区二区三区视频在线| 欧美三级亚洲精品| 男女视频在线观看网站免费| 国产精品亚洲美女久久久| 嫩草影院精品99| 午夜视频国产福利| 亚洲精品粉嫩美女一区| 啦啦啦韩国在线观看视频| 亚洲18禁久久av| 精品一区二区三区av网在线观看| 美女黄网站色视频| 国产aⅴ精品一区二区三区波| 他把我摸到了高潮在线观看| 亚洲第一电影网av| 国产精品,欧美在线| 久久精品国产鲁丝片午夜精品 | 在线国产一区二区在线| 夜夜爽天天搞| 国产乱人视频| 国产精品久久电影中文字幕| 床上黄色一级片| 不卡视频在线观看欧美| 国产精华一区二区三区| 久久人人精品亚洲av| 级片在线观看| 久久久久久国产a免费观看| 成人高潮视频无遮挡免费网站| 性色avwww在线观看| 国产三级中文精品| 蜜桃久久精品国产亚洲av| 国产成人福利小说| 婷婷亚洲欧美| 久99久视频精品免费| 午夜福利欧美成人| 免费在线观看日本一区| 一区二区三区四区激情视频 | 亚洲色图av天堂| 国产真实乱freesex| 亚洲欧美日韩无卡精品| 老女人水多毛片| 乱人视频在线观看| 熟女电影av网| 狂野欧美白嫩少妇大欣赏| 午夜福利欧美成人| av专区在线播放| 亚洲国产精品成人综合色| 色哟哟哟哟哟哟| 蜜桃亚洲精品一区二区三区| 亚洲18禁久久av| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲精华国产精华精| 亚洲人成网站在线播放欧美日韩| 亚洲精品粉嫩美女一区| 毛片女人毛片| 亚洲成人免费电影在线观看| 国产成人影院久久av| 国产老妇女一区| a在线观看视频网站| 亚洲专区国产一区二区| 日本a在线网址| 美女被艹到高潮喷水动态| 黄色女人牲交| 国产单亲对白刺激| 亚洲 国产 在线| 国产亚洲精品久久久com| 中文字幕久久专区| 久久久国产成人精品二区| 婷婷精品国产亚洲av| 精品久久国产蜜桃| 级片在线观看| 一区福利在线观看| 精品午夜福利视频在线观看一区| 日日啪夜夜撸| 一进一出好大好爽视频| 在线免费观看不下载黄p国产 | 日韩精品青青久久久久久| 日韩亚洲欧美综合| 成年版毛片免费区| 国产一区二区亚洲精品在线观看| 美女高潮的动态| 一级黄色大片毛片| 免费人成视频x8x8入口观看| 少妇的逼水好多| 他把我摸到了高潮在线观看| 亚洲一区高清亚洲精品| 亚洲成人久久性| 欧美黑人巨大hd| 不卡视频在线观看欧美| 日韩欧美免费精品| 亚洲国产日韩欧美精品在线观看| 亚洲国产欧美人成| 国产探花在线观看一区二区| 99热网站在线观看| 中文字幕精品亚洲无线码一区| 在线天堂最新版资源| 99久久无色码亚洲精品果冻| 黄色配什么色好看| 国产国拍精品亚洲av在线观看| 极品教师在线视频| 久久久久久久久大av| 人妻丰满熟妇av一区二区三区| 久久久久久久午夜电影| 欧美绝顶高潮抽搐喷水| 99久久无色码亚洲精品果冻| 亚洲中文字幕一区二区三区有码在线看| 在线免费观看不下载黄p国产 | 国产亚洲欧美98| 久久婷婷人人爽人人干人人爱| 午夜精品一区二区三区免费看| 国产精品爽爽va在线观看网站| 国产精品国产高清国产av| 99riav亚洲国产免费| 老熟妇仑乱视频hdxx| 最新在线观看一区二区三区| 精品国产三级普通话版| 成人高潮视频无遮挡免费网站| 国产高清有码在线观看视频| av天堂中文字幕网| 国产成人aa在线观看| 国产 一区精品| 色播亚洲综合网| 国内久久婷婷六月综合欲色啪| 久久精品国产亚洲av香蕉五月| 国产高清三级在线| 国产精品三级大全|