• <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
    免费观看人在逋| 国内精品一区二区在线观看| 国产男靠女视频免费网站| 制服丝袜大香蕉在线| 中文字幕人成人乱码亚洲影| 少妇熟女aⅴ在线视频| 夜夜躁狠狠躁天天躁| 日本三级黄在线观看| 1024手机看黄色片| 狠狠狠狠99中文字幕| 亚洲精品色激情综合| 亚洲国产精品sss在线观看| 精品国产亚洲在线| 12—13女人毛片做爰片一| 国产成人精品久久二区二区免费| 欧美另类亚洲清纯唯美| 国产av一区在线观看免费| 麻豆一二三区av精品| 少妇的丰满在线观看| 老司机福利观看| 黄色片一级片一级黄色片| 国产亚洲av嫩草精品影院| 欧美成人性av电影在线观看| 舔av片在线| 国产主播在线观看一区二区| 在线视频色国产色| 国产精品乱码一区二三区的特点| 成人一区二区视频在线观看| 成年免费大片在线观看| 99riav亚洲国产免费| 琪琪午夜伦伦电影理论片6080| 欧美在线黄色| 我的老师免费观看完整版| 无限看片的www在线观看| 国产高清videossex| 成人性生交大片免费视频hd| 性色av乱码一区二区三区2| 国产91精品成人一区二区三区| 好男人电影高清在线观看| 欧美激情久久久久久爽电影| 高潮久久久久久久久久久不卡| 男人舔奶头视频| 免费av不卡在线播放| 国产熟女xx| 夜夜躁狠狠躁天天躁| 成人亚洲精品av一区二区| 别揉我奶头~嗯~啊~动态视频| 亚洲av第一区精品v没综合| www.999成人在线观看| 波多野结衣高清作品| 久久婷婷人人爽人人干人人爱| 波多野结衣高清无吗| 国产爱豆传媒在线观看| 一级黄色大片毛片| 亚洲欧美日韩东京热| 亚洲国产日韩欧美精品在线观看 | 男女下面进入的视频免费午夜| 国产午夜福利久久久久久| 亚洲va日本ⅴa欧美va伊人久久| 日韩有码中文字幕| 免费人成视频x8x8入口观看| 丰满人妻一区二区三区视频av | 美女扒开内裤让男人捅视频| 天天躁日日操中文字幕| or卡值多少钱| 亚洲av成人av| 婷婷六月久久综合丁香| 嫩草影视91久久| 久久国产精品影院| 搞女人的毛片| 天堂av国产一区二区熟女人妻| 欧美一级a爱片免费观看看| 美女扒开内裤让男人捅视频| 在线观看日韩欧美| 制服丝袜大香蕉在线| 国产高清有码在线观看视频| 久久久久久久午夜电影| 99精品欧美一区二区三区四区| 亚洲国产精品999在线| 午夜福利高清视频| 国产精品1区2区在线观看.| 国产亚洲精品一区二区www| 成人一区二区视频在线观看| 中文字幕久久专区| 精品久久久久久久久久免费视频| 午夜成年电影在线免费观看| 伦理电影免费视频| 精品午夜福利视频在线观看一区| 最近最新中文字幕大全免费视频| 国产精品,欧美在线| 草草在线视频免费看| 精品人妻1区二区| 无遮挡黄片免费观看| 久久久久亚洲av毛片大全| 日本黄色视频三级网站网址| 不卡av一区二区三区| 日日摸夜夜添夜夜添小说| 99热精品在线国产| 国产高清三级在线| 亚洲欧洲精品一区二区精品久久久| 中出人妻视频一区二区| 中文字幕精品亚洲无线码一区| 1024手机看黄色片| 给我免费播放毛片高清在线观看| 母亲3免费完整高清在线观看| 国产成人精品久久二区二区91| 极品教师在线免费播放| 欧美性猛交╳xxx乱大交人| 波多野结衣高清作品| 一进一出好大好爽视频| 一a级毛片在线观看| 久久国产乱子伦精品免费另类| 伊人久久大香线蕉亚洲五| 国产99白浆流出| 欧美日韩综合久久久久久 | 欧美黄色淫秽网站| 三级毛片av免费| 性色av乱码一区二区三区2| 久久久久久久精品吃奶| 国产高清有码在线观看视频| 欧美成人性av电影在线观看| 午夜精品在线福利| 狂野欧美白嫩少妇大欣赏| 国内精品久久久久久久电影| 亚洲国产色片| 亚洲国产精品久久男人天堂| 亚洲 欧美一区二区三区| 韩国av一区二区三区四区| 无限看片的www在线观看| 国产精华一区二区三区| 久久久久九九精品影院| 可以在线观看的亚洲视频| av黄色大香蕉| 国产av麻豆久久久久久久| 99精品在免费线老司机午夜| www日本在线高清视频| 一区二区三区激情视频| 国产精品综合久久久久久久免费| 搡老岳熟女国产| 91麻豆精品激情在线观看国产| 老汉色av国产亚洲站长工具| 久久久国产精品麻豆| 亚洲精品中文字幕一二三四区| 黄频高清免费视频| 18美女黄网站色大片免费观看| 麻豆成人午夜福利视频| 久久精品国产99精品国产亚洲性色| 欧洲精品卡2卡3卡4卡5卡区| 最新美女视频免费是黄的| 亚洲午夜理论影院| 全区人妻精品视频| 97碰自拍视频| 中文字幕人成人乱码亚洲影| 亚洲av美国av| 露出奶头的视频| a级毛片在线看网站| 一本一本综合久久| 很黄的视频免费| 色视频www国产| www.熟女人妻精品国产| 国内久久婷婷六月综合欲色啪| 欧美大码av| 麻豆一二三区av精品| 欧美不卡视频在线免费观看| 夜夜看夜夜爽夜夜摸| 国产三级在线视频| 国产91精品成人一区二区三区| av在线天堂中文字幕| 日韩av在线大香蕉| 麻豆成人午夜福利视频| www.熟女人妻精品国产| 国产1区2区3区精品| 老司机深夜福利视频在线观看| 国产精品98久久久久久宅男小说| 日韩免费av在线播放| 久久人人精品亚洲av| 国产av不卡久久| 亚洲精华国产精华精| 亚洲国产色片| 99久久精品国产亚洲精品| 欧美丝袜亚洲另类 | 欧美日韩福利视频一区二区| 两人在一起打扑克的视频| 精品久久久久久成人av| 中文字幕高清在线视频| 91老司机精品| 看片在线看免费视频| 国产视频一区二区在线看| 亚洲自拍偷在线| 无遮挡黄片免费观看| av欧美777| 久久久久久久久久黄片| 天堂√8在线中文| 噜噜噜噜噜久久久久久91| 在线观看美女被高潮喷水网站 | 亚洲成a人片在线一区二区| 精品午夜福利视频在线观看一区| 蜜桃久久精品国产亚洲av| 久久久久久久久免费视频了| 亚洲 欧美一区二区三区| 久久精品国产综合久久久| 欧美激情久久久久久爽电影| 在线a可以看的网站| 久久久久久人人人人人| 一进一出好大好爽视频| 嫩草影院精品99| 久9热在线精品视频| 成年女人毛片免费观看观看9| 婷婷丁香在线五月| 99精品在免费线老司机午夜| 网址你懂的国产日韩在线| 身体一侧抽搐| 91字幕亚洲| 一a级毛片在线观看| 变态另类丝袜制服| 欧美黑人巨大hd| 日本黄色片子视频| 中文字幕最新亚洲高清| 99国产精品99久久久久| 亚洲国产精品合色在线| 老汉色∧v一级毛片| АⅤ资源中文在线天堂| 亚洲黑人精品在线| 女警被强在线播放| 精品久久久久久久久久久久久| 中文字幕av在线有码专区| 国产精品1区2区在线观看.| 亚洲 国产 在线| 日本 av在线| 亚洲熟妇熟女久久| 国产精品久久久久久人妻精品电影| 国产精品美女特级片免费视频播放器 | cao死你这个sao货| 午夜精品在线福利| 久久精品国产99精品国产亚洲性色| 99热6这里只有精品| 久99久视频精品免费| 成年免费大片在线观看| 亚洲精品乱码久久久v下载方式 | 国产v大片淫在线免费观看| 国产av一区在线观看免费| 中文字幕高清在线视频| 色老头精品视频在线观看| 久久久久久人人人人人| av国产免费在线观看| 久久久久久国产a免费观看| 最近最新免费中文字幕在线| 欧美丝袜亚洲另类 | 国产成人啪精品午夜网站| 成人一区二区视频在线观看| 色老头精品视频在线观看| 日韩免费av在线播放| 日韩欧美一区二区三区在线观看| 国产蜜桃级精品一区二区三区| 亚洲av电影在线进入| a在线观看视频网站| 国产成人欧美在线观看| 伦理电影免费视频| 国产精品亚洲美女久久久| 久久欧美精品欧美久久欧美| 美女大奶头视频| 国产成年人精品一区二区| 19禁男女啪啪无遮挡网站| 色精品久久人妻99蜜桃| 亚洲人成伊人成综合网2020| 成人亚洲精品av一区二区| 国产成人av激情在线播放| 亚洲中文日韩欧美视频| 网址你懂的国产日韩在线| 日本免费一区二区三区高清不卡| 午夜a级毛片| a级毛片a级免费在线| 啪啪无遮挡十八禁网站| 欧美成狂野欧美在线观看| 亚洲一区高清亚洲精品| 国产1区2区3区精品| 亚洲精品一卡2卡三卡4卡5卡| 最新中文字幕久久久久 | 舔av片在线| 亚洲国产看品久久| 精品国产乱码久久久久久男人| 亚洲人成电影免费在线| 亚洲 欧美一区二区三区| 国产精品影院久久| 国产高清视频在线观看网站| 69av精品久久久久久| 国产单亲对白刺激| 欧美乱码精品一区二区三区| 一进一出抽搐gif免费好疼| 日韩 欧美 亚洲 中文字幕| av视频在线观看入口| 老司机福利观看| 成人国产综合亚洲| 国产高清激情床上av| 国产黄色小视频在线观看| 国产成人av教育| 99热这里只有是精品50| 欧美日韩综合久久久久久 | 看黄色毛片网站| 动漫黄色视频在线观看| 久久精品aⅴ一区二区三区四区| 亚洲九九香蕉| 免费在线观看成人毛片| 色综合婷婷激情| 好男人电影高清在线观看| 99热只有精品国产| 久久九九热精品免费| 亚洲人成网站高清观看| 欧美一区二区精品小视频在线| 黄色视频,在线免费观看| 免费观看的影片在线观看| 亚洲第一欧美日韩一区二区三区| 亚洲中文日韩欧美视频| 日本黄大片高清| 欧美黄色片欧美黄色片| 男女视频在线观看网站免费| 日本撒尿小便嘘嘘汇集6| 精品国产乱子伦一区二区三区| 精品不卡国产一区二区三区| 波多野结衣高清无吗| 久久香蕉国产精品| 欧美成狂野欧美在线观看| 丰满人妻熟妇乱又伦精品不卡| 无限看片的www在线观看| 欧美又色又爽又黄视频| 国产亚洲欧美在线一区二区| 午夜福利视频1000在线观看| svipshipincom国产片| 欧美精品啪啪一区二区三区| 亚洲成人精品中文字幕电影| 九九热线精品视视频播放| 午夜精品久久久久久毛片777| 男女那种视频在线观看| 一二三四在线观看免费中文在| 亚洲欧美日韩无卡精品| 久久精品aⅴ一区二区三区四区| 久久中文字幕一级| 国产麻豆成人av免费视频| 亚洲成av人片在线播放无| av欧美777| 国产精品美女特级片免费视频播放器 | 国产成人精品久久二区二区免费| 两个人的视频大全免费| 精品乱码久久久久久99久播| 1024手机看黄色片| 可以在线观看毛片的网站| 嫩草影视91久久| 欧美一区二区国产精品久久精品| 丝袜人妻中文字幕| 免费在线观看日本一区| 麻豆久久精品国产亚洲av| 宅男免费午夜| av视频在线观看入口| 国产久久久一区二区三区| 在线a可以看的网站| 久久这里只有精品中国| 亚洲精华国产精华精| 亚洲国产精品久久男人天堂| 婷婷亚洲欧美| 久久精品国产亚洲av香蕉五月| 日本一二三区视频观看| 亚洲人成网站高清观看| 人妻夜夜爽99麻豆av| 日韩三级视频一区二区三区| 成年人黄色毛片网站| 亚洲av成人精品一区久久| 两个人的视频大全免费| 国产精品日韩av在线免费观看| 精品熟女少妇八av免费久了| 欧美日本视频| 午夜福利18| 99国产精品一区二区蜜桃av| 欧美日韩亚洲国产一区二区在线观看| 精品福利观看| 长腿黑丝高跟| 国产成人欧美在线观看| 99国产精品一区二区三区| 精品国产三级普通话版| 好男人在线观看高清免费视频| 久久精品国产综合久久久| 少妇丰满av| 国产亚洲精品久久久com| 啦啦啦观看免费观看视频高清| 国内少妇人妻偷人精品xxx网站 | 欧美大码av| 精品电影一区二区在线| 淫秽高清视频在线观看| 国产激情久久老熟女| 99视频精品全部免费 在线 | 级片在线观看| 每晚都被弄得嗷嗷叫到高潮| 极品教师在线免费播放| 成人18禁在线播放| 伊人久久大香线蕉亚洲五| 亚洲av电影在线进入| 亚洲欧美精品综合久久99| 老熟妇乱子伦视频在线观看| 色综合亚洲欧美另类图片| 日本一二三区视频观看| 精品国内亚洲2022精品成人| 好男人电影高清在线观看| 久久中文字幕人妻熟女| 色综合欧美亚洲国产小说| 婷婷亚洲欧美| 中文字幕人妻丝袜一区二区| 成人av在线播放网站| 一本综合久久免费| 亚洲av美国av| 特大巨黑吊av在线直播| 亚洲国产高清在线一区二区三| 男女之事视频高清在线观看| 中出人妻视频一区二区| 一级黄色大片毛片| 亚洲国产欧美人成| 亚洲av电影不卡..在线观看| 麻豆久久精品国产亚洲av| 中文字幕熟女人妻在线| 99久国产av精品| 精品国产超薄肉色丝袜足j| 午夜亚洲福利在线播放| 日本一本二区三区精品| 999久久久精品免费观看国产| 舔av片在线| 99热这里只有精品一区 | 99re在线观看精品视频| 俺也久久电影网| 看免费av毛片| 亚洲国产看品久久| 久久久色成人| 成人永久免费在线观看视频| 丰满的人妻完整版| 国产美女午夜福利| 国语自产精品视频在线第100页| 麻豆国产97在线/欧美| www.精华液| 国产精品综合久久久久久久免费| 亚洲黑人精品在线| 窝窝影院91人妻| 搡老岳熟女国产| 亚洲色图 男人天堂 中文字幕| 日韩欧美免费精品| 国产97色在线日韩免费| 亚洲中文字幕一区二区三区有码在线看 | 丝袜人妻中文字幕| 亚洲成av人片在线播放无| 丰满人妻一区二区三区视频av | 九九热线精品视视频播放| 亚洲欧洲精品一区二区精品久久久| 又粗又爽又猛毛片免费看| 男女视频在线观看网站免费| 一级毛片高清免费大全| 99久久精品热视频| 在线观看免费午夜福利视频| 亚洲最大成人中文| 午夜影院日韩av| 美女cb高潮喷水在线观看 | 亚洲精品美女久久久久99蜜臀| 午夜激情欧美在线| 免费av不卡在线播放| 久久久国产精品麻豆| 级片在线观看| 神马国产精品三级电影在线观看| 亚洲黑人精品在线| 中文资源天堂在线| 九色成人免费人妻av| 亚洲av电影不卡..在线观看| 国产精品综合久久久久久久免费| 观看美女的网站| 99精品在免费线老司机午夜| 在线观看一区二区三区| 久久这里只有精品19| 成年女人看的毛片在线观看| 制服丝袜大香蕉在线| 中文在线观看免费www的网站| 久久久国产欧美日韩av| 麻豆国产av国片精品| 成人高潮视频无遮挡免费网站| www.999成人在线观看| 成人三级做爰电影| 精品电影一区二区在线| 性色avwww在线观看| 日本熟妇午夜| 一个人免费在线观看的高清视频| 天堂影院成人在线观看| 午夜福利在线在线| 91在线精品国自产拍蜜月 | 亚洲国产中文字幕在线视频| 一进一出抽搐gif免费好疼| 久久久久久久精品吃奶| 一区二区三区国产精品乱码| 国产蜜桃级精品一区二区三区| 亚洲人成伊人成综合网2020| 精品不卡国产一区二区三区| 亚洲九九香蕉| 亚洲在线观看片| 少妇丰满av| 中文字幕精品亚洲无线码一区| 小说图片视频综合网站| 久久中文字幕人妻熟女| 人妻丰满熟妇av一区二区三区| 99久久久亚洲精品蜜臀av| 国产精品久久久久久人妻精品电影| 窝窝影院91人妻| 他把我摸到了高潮在线观看| 国产97色在线日韩免费| 老汉色∧v一级毛片| 这个男人来自地球电影免费观看| 亚洲国产精品合色在线| 国产三级黄色录像| 高清在线国产一区| 五月伊人婷婷丁香| 成人鲁丝片一二三区免费| 成人av一区二区三区在线看| 欧美成人一区二区免费高清观看 | 中文字幕人妻丝袜一区二区| 成人亚洲精品av一区二区| 日本五十路高清| 99在线视频只有这里精品首页| 国产毛片a区久久久久| 欧美性猛交╳xxx乱大交人| 亚洲av成人一区二区三| 国产一区二区在线av高清观看| 欧洲精品卡2卡3卡4卡5卡区| 2021天堂中文幕一二区在线观| 男人和女人高潮做爰伦理| 在线a可以看的网站| 69av精品久久久久久| 国产免费男女视频| 嫩草影院精品99| 99久久精品一区二区三区| 88av欧美| 国产 一区 欧美 日韩| 久久99热这里只有精品18| 欧美另类亚洲清纯唯美| 成人三级黄色视频| x7x7x7水蜜桃| 男人舔女人的私密视频| tocl精华| 精品国产超薄肉色丝袜足j| a在线观看视频网站| 成人国产综合亚洲| 18禁美女被吸乳视频| 国产精品亚洲美女久久久| 搡老妇女老女人老熟妇| 国内精品一区二区在线观看| 国产精品av久久久久免费| 成人精品一区二区免费| 99在线人妻在线中文字幕| 人妻丰满熟妇av一区二区三区| 美女高潮喷水抽搐中文字幕| 观看免费一级毛片| 欧美午夜高清在线| 中文字幕最新亚洲高清| 日韩欧美一区二区三区在线观看| 日本在线视频免费播放| 在线永久观看黄色视频| 国产又黄又爽又无遮挡在线| 国产伦精品一区二区三区视频9 | a级毛片在线看网站| 一二三四在线观看免费中文在| 最近最新免费中文字幕在线| 欧美日韩综合久久久久久 | 亚洲中文av在线| 中文字幕av在线有码专区| 综合色av麻豆| 国产亚洲精品av在线| 久久久精品欧美日韩精品| 亚洲av成人一区二区三| 99国产综合亚洲精品| 国产1区2区3区精品| 90打野战视频偷拍视频| 欧美高清成人免费视频www| 国产精品爽爽va在线观看网站| 91在线观看av| 香蕉丝袜av| 国产精品乱码一区二三区的特点| 亚洲精品中文字幕一二三四区| 亚洲国产精品合色在线| 99国产精品一区二区蜜桃av| 黄色日韩在线| a在线观看视频网站| 国产午夜福利久久久久久| 国产淫片久久久久久久久 | 国产精品精品国产色婷婷| 亚洲国产日韩欧美精品在线观看 | 亚洲精品国产精品久久久不卡| 国产成人av激情在线播放| 久久久久久久久久黄片| 国产又黄又爽又无遮挡在线| 丰满人妻一区二区三区视频av | 午夜免费激情av| 亚洲中文字幕日韩| 美女cb高潮喷水在线观看 | 亚洲精品久久国产高清桃花| 国产精品亚洲一级av第二区| 俺也久久电影网| 午夜福利在线观看免费完整高清在 | 美女高潮喷水抽搐中文字幕| 国内精品久久久久精免费| 日韩三级视频一区二区三区| 真人做人爱边吃奶动态| a级毛片a级免费在线| 色综合婷婷激情| av在线天堂中文字幕| 一个人看视频在线观看www免费 | 精品久久久久久久毛片微露脸| 成人精品一区二区免费| 淫秽高清视频在线观看| 搡老熟女国产l中国老女人| 成人午夜高清在线视频| 黄色日韩在线| 国产淫片久久久久久久久 | 亚洲无线在线观看| 搡老熟女国产l中国老女人| 一夜夜www|