• <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
    小蜜桃在线观看免费完整版高清| www.色视频.com| 一夜夜www| 人妻少妇偷人精品九色| 真实男女啪啪啪动态图| 一级黄色大片毛片| 国产高清激情床上av| 深爱激情五月婷婷| 男人狂女人下面高潮的视频| 久久久久免费精品人妻一区二区| 日韩av不卡免费在线播放| 一进一出好大好爽视频| 一本久久中文字幕| 欧美一级a爱片免费观看看| .国产精品久久| 国产三级在线视频| 午夜亚洲福利在线播放| 精品日产1卡2卡| 亚洲成av人片在线播放无| 我的老师免费观看完整版| 久久久久久伊人网av| 男女视频在线观看网站免费| 一本精品99久久精品77| 免费观看在线日韩| 国产精品一区二区三区四区久久| 麻豆成人午夜福利视频| 美女被艹到高潮喷水动态| 久久精品91蜜桃| 久久婷婷人人爽人人干人人爱| 亚洲av第一区精品v没综合| 国产亚洲欧美98| 亚洲人成网站高清观看| 免费在线观看影片大全网站| 免费不卡的大黄色大毛片视频在线观看 | 日韩在线高清观看一区二区三区| 久久久久国产网址| 免费看a级黄色片| 色噜噜av男人的天堂激情| 一a级毛片在线观看| 国产伦精品一区二区三区视频9| 精品免费久久久久久久清纯| 干丝袜人妻中文字幕| av在线亚洲专区| 人人妻人人澡欧美一区二区| 欧美性猛交╳xxx乱大交人| 美女 人体艺术 gogo| 99热精品在线国产| av中文乱码字幕在线| 久久久久久大精品| 97人妻精品一区二区三区麻豆| 国产精品日韩av在线免费观看| 在线免费十八禁| 人妻丰满熟妇av一区二区三区| 亚洲电影在线观看av| 国产私拍福利视频在线观看| 亚洲国产精品成人久久小说 | 女同久久另类99精品国产91| 在线免费观看的www视频| 性色avwww在线观看| 精品人妻一区二区三区麻豆 | 免费观看人在逋| 成年免费大片在线观看| 国内精品一区二区在线观看| 免费av毛片视频| 一边摸一边抽搐一进一小说| 国产成人a区在线观看| .国产精品久久| 村上凉子中文字幕在线| av在线观看视频网站免费| 国产午夜精品久久久久久一区二区三区 | 免费不卡的大黄色大毛片视频在线观看 | 午夜福利在线在线| 久久欧美精品欧美久久欧美| 欧美日韩精品成人综合77777| 九九久久精品国产亚洲av麻豆| www.色视频.com| 亚洲精品粉嫩美女一区| 国产精品女同一区二区软件| 国产成人a∨麻豆精品| 久久国内精品自在自线图片| 亚洲欧美中文字幕日韩二区| 99久国产av精品| 久久久久久久亚洲中文字幕| 此物有八面人人有两片| 69人妻影院| 人妻久久中文字幕网| 三级经典国产精品| www.色视频.com| 最近手机中文字幕大全| 99久久精品一区二区三区| 亚洲激情五月婷婷啪啪| 深爱激情五月婷婷| 欧美国产日韩亚洲一区| 成人av一区二区三区在线看| 亚洲av一区综合| 晚上一个人看的免费电影| 亚洲色图av天堂| 一区福利在线观看| 免费观看的影片在线观看| 淫妇啪啪啪对白视频| av在线蜜桃| 亚洲精品粉嫩美女一区| 晚上一个人看的免费电影| 国产黄片美女视频| АⅤ资源中文在线天堂| 成年免费大片在线观看| 99视频精品全部免费 在线| 国产成年人精品一区二区| 国产人妻一区二区三区在| 小蜜桃在线观看免费完整版高清| 在线观看一区二区三区| 日本欧美国产在线视频| 欧美日本视频| 91午夜精品亚洲一区二区三区| 22中文网久久字幕| 欧美极品一区二区三区四区| 精品久久久久久成人av| 男女那种视频在线观看| 欧美绝顶高潮抽搐喷水| 日本黄大片高清| 日韩制服骚丝袜av| 欧美zozozo另类| 一卡2卡三卡四卡精品乱码亚洲| 人妻制服诱惑在线中文字幕| 精华霜和精华液先用哪个| 日本a在线网址| 成年女人毛片免费观看观看9| 熟妇人妻久久中文字幕3abv| 一级a爱片免费观看的视频| 最近中文字幕高清免费大全6| 极品教师在线视频| 不卡视频在线观看欧美| 国产不卡一卡二| 久久天躁狠狠躁夜夜2o2o| 特大巨黑吊av在线直播| 丝袜喷水一区| 又黄又爽又免费观看的视频| 中文字幕av成人在线电影| 长腿黑丝高跟| 一个人观看的视频www高清免费观看| 免费观看的影片在线观看| 最近的中文字幕免费完整| 美女被艹到高潮喷水动态| 在线观看66精品国产| 日本与韩国留学比较| 久久久久久久午夜电影| av视频在线观看入口| 欧美色欧美亚洲另类二区| 精品国内亚洲2022精品成人| 亚洲专区国产一区二区| 成人三级黄色视频| 亚洲熟妇中文字幕五十中出| 一区二区三区免费毛片| 国产69精品久久久久777片| 尤物成人国产欧美一区二区三区| 婷婷亚洲欧美| 日韩人妻高清精品专区| 国产精品免费一区二区三区在线| 最后的刺客免费高清国语| 韩国av在线不卡| 欧美最黄视频在线播放免费| 看十八女毛片水多多多| 国产男靠女视频免费网站| 黑人高潮一二区| 18禁在线无遮挡免费观看视频 | 搡老熟女国产l中国老女人| .国产精品久久| 国产精品福利在线免费观看| 高清日韩中文字幕在线| 在线播放无遮挡| 人妻制服诱惑在线中文字幕| 又粗又爽又猛毛片免费看| 婷婷色综合大香蕉| 色哟哟·www| 久久久久国内视频| 免费看av在线观看网站| 啦啦啦韩国在线观看视频| 婷婷精品国产亚洲av| 免费看a级黄色片| 91久久精品国产一区二区三区| 亚洲欧美清纯卡通| 色av中文字幕| 国产乱人视频| 99riav亚洲国产免费| 小蜜桃在线观看免费完整版高清| 久久精品久久久久久噜噜老黄 | 成人午夜高清在线视频| 成人精品一区二区免费| 蜜臀久久99精品久久宅男| 香蕉av资源在线| 国产精品一区二区三区四区久久| 在线观看午夜福利视频| 高清毛片免费看| 亚洲第一区二区三区不卡| АⅤ资源中文在线天堂| 狂野欧美白嫩少妇大欣赏| 99热精品在线国产| 夜夜看夜夜爽夜夜摸| 午夜影院日韩av| 九色成人免费人妻av| 99久久无色码亚洲精品果冻| 久久久国产成人精品二区| 人人妻人人澡人人爽人人夜夜 | 99久国产av精品| 身体一侧抽搐| 精品久久久久久久末码| 国产成人一区二区在线| 免费不卡的大黄色大毛片视频在线观看 | 成年版毛片免费区| 听说在线观看完整版免费高清| 亚洲图色成人| 国产三级中文精品| 国产亚洲精品久久久com| 国产亚洲精品av在线| 久久久成人免费电影| 午夜激情欧美在线| 色哟哟哟哟哟哟| 精品免费久久久久久久清纯| 99久久精品一区二区三区| avwww免费| 午夜爱爱视频在线播放| 欧美日韩一区二区视频在线观看视频在线 | 久久精品国产亚洲av香蕉五月| 欧美xxxx性猛交bbbb| 白带黄色成豆腐渣| 免费大片18禁| 精品久久久久久久人妻蜜臀av| 婷婷六月久久综合丁香| 啦啦啦啦在线视频资源| 亚洲第一区二区三区不卡| 亚洲av第一区精品v没综合| 免费观看精品视频网站| av免费在线看不卡| 午夜精品国产一区二区电影 | 国产真实伦视频高清在线观看| 国产亚洲av嫩草精品影院| 成人永久免费在线观看视频| 亚洲熟妇中文字幕五十中出| 日韩欧美国产在线观看| 精品久久久久久久久亚洲| 麻豆国产av国片精品| 成年版毛片免费区| 搡老妇女老女人老熟妇| 18+在线观看网站| 国产精品免费一区二区三区在线| 亚洲美女搞黄在线观看 | 免费黄网站久久成人精品| 天天躁夜夜躁狠狠久久av| 一区二区三区免费毛片| 国产精品国产高清国产av| 亚洲精品国产成人久久av| 老师上课跳d突然被开到最大视频| 2021天堂中文幕一二区在线观| 日本欧美国产在线视频| 亚洲人成网站在线播放欧美日韩| 91久久精品国产一区二区成人| 日日摸夜夜添夜夜添小说| 国产精品久久久久久精品电影| 美女高潮的动态| 不卡一级毛片| 国模一区二区三区四区视频| 久久鲁丝午夜福利片| 国产一区二区亚洲精品在线观看| 观看免费一级毛片| 婷婷色综合大香蕉| 欧美3d第一页| 人人妻,人人澡人人爽秒播| 国产亚洲av嫩草精品影院| 国产毛片a区久久久久| 中文字幕熟女人妻在线| 99久久精品一区二区三区| 一进一出好大好爽视频| 亚洲精品乱码久久久v下载方式| 亚洲精品亚洲一区二区| 久99久视频精品免费| 九九在线视频观看精品| 亚洲精品影视一区二区三区av| 搡老熟女国产l中国老女人| 亚洲熟妇熟女久久| 欧美一级a爱片免费观看看| 搡女人真爽免费视频火全软件 | 精品人妻偷拍中文字幕| 欧美色视频一区免费| 婷婷精品国产亚洲av| 在线观看美女被高潮喷水网站| 国产av麻豆久久久久久久| 国产精品一及| 精品久久久久久久久久久久久| 18禁黄网站禁片免费观看直播| 最近的中文字幕免费完整| 永久网站在线| 禁无遮挡网站| 高清午夜精品一区二区三区 | 欧美成人精品欧美一级黄| 亚洲av中文av极速乱| 午夜福利视频1000在线观看| 熟女电影av网| 久久久久九九精品影院| 插逼视频在线观看| 99riav亚洲国产免费| 伦理电影大哥的女人| 我要搜黄色片| 欧美人与善性xxx| 国产黄a三级三级三级人| 老师上课跳d突然被开到最大视频| 久久精品国产自在天天线| 亚洲成a人片在线一区二区| 国产爱豆传媒在线观看| 国产精品一区www在线观看| 国产亚洲精品久久久久久毛片| av在线播放精品| 日韩成人伦理影院| 亚洲性夜色夜夜综合| 免费观看精品视频网站| 免费高清视频大片| 国产精品一区二区免费欧美| 国产日本99.免费观看| 夜夜看夜夜爽夜夜摸| 亚洲高清免费不卡视频| 亚洲经典国产精华液单| 美女内射精品一级片tv| 黄色日韩在线| 18禁在线无遮挡免费观看视频 | 亚洲高清免费不卡视频| 亚洲经典国产精华液单| 卡戴珊不雅视频在线播放| 日日摸夜夜添夜夜添av毛片| 国产精品精品国产色婷婷| 可以在线观看毛片的网站| 免费看a级黄色片| 岛国在线免费视频观看| 亚洲在线自拍视频| 亚洲专区国产一区二区| 欧美人与善性xxx| av在线亚洲专区| 国产av不卡久久| 熟女人妻精品中文字幕| 小说图片视频综合网站| 男女边吃奶边做爰视频| 久久精品国产99精品国产亚洲性色| 三级男女做爰猛烈吃奶摸视频| 一级毛片久久久久久久久女| 亚洲av成人精品一区久久| 亚洲真实伦在线观看| 天天躁日日操中文字幕| 真人做人爱边吃奶动态| 精品一区二区免费观看| 欧美xxxx性猛交bbbb| 欧美激情国产日韩精品一区| 听说在线观看完整版免费高清| 午夜免费男女啪啪视频观看 | 青春草视频在线免费观看| 国产精品久久久久久精品电影| 亚洲图色成人| 日日撸夜夜添| 中文字幕人妻熟人妻熟丝袜美| 久久亚洲精品不卡| av黄色大香蕉| 五月玫瑰六月丁香| 亚洲丝袜综合中文字幕| 久久久久久大精品| 18+在线观看网站| a级毛片a级免费在线| 亚洲中文日韩欧美视频| 99九九线精品视频在线观看视频| 国产成人一区二区在线| 两个人视频免费观看高清| 内射极品少妇av片p| 日韩欧美免费精品| 成人高潮视频无遮挡免费网站| 噜噜噜噜噜久久久久久91| 国产精品人妻久久久久久| 日韩强制内射视频| 日韩成人伦理影院| 狠狠狠狠99中文字幕| 伊人久久精品亚洲午夜| 欧美zozozo另类| 色综合色国产| 精品久久久久久久末码| 久久精品91蜜桃| 禁无遮挡网站| 老熟妇仑乱视频hdxx| 99在线视频只有这里精品首页| 亚洲一区二区三区色噜噜| 好男人在线观看高清免费视频| 精品一区二区三区人妻视频| 51国产日韩欧美| 日韩一本色道免费dvd| 狂野欧美白嫩少妇大欣赏| 日韩精品青青久久久久久| 热99re8久久精品国产| 欧美又色又爽又黄视频| 久久精品国产99精品国产亚洲性色| 国产 一区精品| 成熟少妇高潮喷水视频| 久久久久久久久大av| 观看免费一级毛片| 精品久久久久久久久av| 日韩成人伦理影院| 亚洲欧美日韩卡通动漫| 精品不卡国产一区二区三区| 亚洲欧美日韩高清专用| 波多野结衣高清无吗| 夜夜夜夜夜久久久久| 亚洲美女黄片视频| 国产精品久久久久久久久免| 日韩制服骚丝袜av| 精品久久久噜噜| 国产亚洲91精品色在线| 欧美高清成人免费视频www| 美女xxoo啪啪120秒动态图| 精品少妇黑人巨大在线播放 | 中文字幕熟女人妻在线| 亚洲欧美日韩东京热| 免费人成视频x8x8入口观看| 成人亚洲欧美一区二区av| 在线观看午夜福利视频| 亚洲成人久久性| 最近的中文字幕免费完整| 蜜臀久久99精品久久宅男| 久久精品国产99精品国产亚洲性色| 精品久久久久久久久亚洲| 亚洲美女搞黄在线观看 | 久久久欧美国产精品| 国产精品女同一区二区软件| 亚洲色图av天堂| 插阴视频在线观看视频| 日韩欧美 国产精品| 一级黄片播放器| 亚洲自拍偷在线| 精品99又大又爽又粗少妇毛片| 国产成人91sexporn| 午夜精品一区二区三区免费看| 哪里可以看免费的av片| 中国美白少妇内射xxxbb| 一夜夜www| .国产精品久久| 国产亚洲精品久久久com| 99久久九九国产精品国产免费| 日本精品一区二区三区蜜桃| 国产高清三级在线| 全区人妻精品视频| 成人无遮挡网站| 国产在线男女| 日韩三级伦理在线观看| 久久精品国产亚洲av天美| 一级a爱片免费观看的视频| 欧美丝袜亚洲另类| 午夜福利高清视频| 国产成人aa在线观看| 99久久久亚洲精品蜜臀av| 麻豆成人午夜福利视频| 亚洲中文字幕一区二区三区有码在线看| 99riav亚洲国产免费| 国产精品1区2区在线观看.| av卡一久久| 国产午夜精品论理片| 精品久久久久久久久av| 黄色欧美视频在线观看| 久久中文看片网| 99热这里只有精品一区| av在线观看视频网站免费| 51国产日韩欧美| 大型黄色视频在线免费观看| 最近手机中文字幕大全| 欧美区成人在线视频| 俺也久久电影网| 黄色一级大片看看| or卡值多少钱| 看黄色毛片网站| 天堂√8在线中文| 国产综合懂色| 男女做爰动态图高潮gif福利片| 色视频www国产| 嫩草影视91久久| 日本欧美国产在线视频| 亚洲国产高清在线一区二区三| 亚洲精品久久国产高清桃花| 国产午夜精品论理片| a级毛片免费高清观看在线播放| 精品久久久噜噜| 亚洲国产欧洲综合997久久,| 真实男女啪啪啪动态图| 久久精品国产清高在天天线| 亚洲国产精品国产精品| 高清日韩中文字幕在线| 午夜久久久久精精品| 免费电影在线观看免费观看| 日本一本二区三区精品| 丰满人妻一区二区三区视频av| 亚洲久久久久久中文字幕| 午夜福利在线在线| 亚洲四区av| 人人妻人人澡人人爽人人夜夜 | 1000部很黄的大片| 成人无遮挡网站| h日本视频在线播放| АⅤ资源中文在线天堂| 淫妇啪啪啪对白视频| 亚洲一级一片aⅴ在线观看| 老司机福利观看| 日韩av不卡免费在线播放| 久久久久九九精品影院| 人妻制服诱惑在线中文字幕| 日本免费a在线| 精品久久久久久久人妻蜜臀av| 天堂√8在线中文| 国产乱人偷精品视频| 村上凉子中文字幕在线| 亚洲国产精品成人久久小说 | 久久久久国产精品人妻aⅴ院| 五月玫瑰六月丁香| 一级av片app| 国产精品一区二区性色av| 天美传媒精品一区二区| 欧美潮喷喷水| 99在线人妻在线中文字幕| 国产av一区在线观看免费| 亚洲欧美日韩无卡精品| 美女 人体艺术 gogo| 老熟妇仑乱视频hdxx| 久久精品国产亚洲网站| 国产三级在线视频| 成年女人毛片免费观看观看9| 亚洲精品一区av在线观看| 亚洲av免费高清在线观看| 三级经典国产精品| 综合色丁香网| 日韩一区二区视频免费看| 亚洲欧美成人综合另类久久久 | 熟女电影av网| 亚洲四区av| 男女啪啪激烈高潮av片| 久99久视频精品免费| 亚洲成人久久性| 俺也久久电影网| 亚洲经典国产精华液单| 一夜夜www| 91在线观看av| 99在线人妻在线中文字幕| 国产大屁股一区二区在线视频| 美女大奶头视频| 久久久久性生活片| 欧美一区二区精品小视频在线| 老司机影院成人| h日本视频在线播放| 日本爱情动作片www.在线观看 | 国产精品av视频在线免费观看| 少妇人妻精品综合一区二区 | 欧洲精品卡2卡3卡4卡5卡区| 久久久成人免费电影| 亚洲图色成人| 久久亚洲精品不卡| 国产91av在线免费观看| 精华霜和精华液先用哪个| 99热这里只有精品一区| 高清日韩中文字幕在线| 亚洲美女搞黄在线观看 | 高清毛片免费看| 国产蜜桃级精品一区二区三区| 在线播放国产精品三级| 菩萨蛮人人尽说江南好唐韦庄 | 国产亚洲精品久久久久久毛片| 日日摸夜夜添夜夜爱| 熟妇人妻久久中文字幕3abv| a级毛片a级免费在线| 日韩一区二区视频免费看| 国产亚洲91精品色在线| 老女人水多毛片| 婷婷六月久久综合丁香| 色av中文字幕| 欧美日韩精品成人综合77777| 欧美3d第一页| 自拍偷自拍亚洲精品老妇| 精品无人区乱码1区二区| 中文字幕av成人在线电影| 亚洲国产色片| 1000部很黄的大片| 一本一本综合久久| 日本-黄色视频高清免费观看| 日韩高清综合在线| 小说图片视频综合网站| 六月丁香七月| 国产精品1区2区在线观看.| 国产成人a∨麻豆精品| 少妇高潮的动态图| 亚洲七黄色美女视频| 日本熟妇午夜| 狠狠狠狠99中文字幕| 国产精品乱码一区二三区的特点| 亚洲av一区综合| 伦精品一区二区三区| 久久6这里有精品| 欧美最新免费一区二区三区| 午夜精品在线福利| 亚洲中文字幕一区二区三区有码在线看| 亚洲欧美清纯卡通| 欧美日韩一区二区视频在线观看视频在线 | 丰满乱子伦码专区| a级一级毛片免费在线观看| 国国产精品蜜臀av免费| 午夜视频国产福利| 丰满乱子伦码专区| 免费av观看视频| 俄罗斯特黄特色一大片| 欧美日韩国产亚洲二区| 日韩欧美免费精品| 伊人久久精品亚洲午夜| 亚洲av熟女| 国语自产精品视频在线第100页| 一个人观看的视频www高清免费观看| 老熟妇乱子伦视频在线观看| 欧美成人a在线观看| 日本成人三级电影网站| 高清毛片免费观看视频网站|