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

    考慮等效擬合特性的輸配電網(wǎng)全局潮流計算方法

    2018-08-02 01:12:18俞龍飛丘國斌歐陽森
    現(xiàn)代電力 2018年4期
    關(guān)鍵詞:輸電網(wǎng)邊界點主從

    俞龍飛,曾 江,丘國斌,歐陽森

    (華南理工大學(xué)電力學(xué)院,廣東廣州 510640)

    0 引 言

    傳統(tǒng)的將輸電網(wǎng)與配電網(wǎng)割裂開來進行潮流計算的方法已無法滿足現(xiàn)代電力系統(tǒng)運行分析和規(guī)劃管理的需求,將輸電網(wǎng)和配電網(wǎng)潮流進行統(tǒng)一分析是當前的發(fā)展趨勢[1],尤其是在輸配網(wǎng)電壓無功的統(tǒng)一優(yōu)化、輸配網(wǎng)網(wǎng)架重構(gòu)等方面的研究與應(yīng)用,此時調(diào)控變量不斷變化,邊界點的功率未知,必須進行多次的輸配網(wǎng)聯(lián)合潮流計算進行尋優(yōu)。由于輸配網(wǎng)全局潮流計算側(cè)重點不同,規(guī)模龐大,通過減少計算量來提高收斂速度是頗具理論與工程價值的研究內(nèi)容。

    大規(guī)模電力系統(tǒng)潮流計算實質(zhì)上是求解一組非線性代數(shù)方程組,采用分布式潮流計算方法[2-4]將大電網(wǎng)分為多個子電網(wǎng),把大規(guī)模方程組的求解轉(zhuǎn)化為多個規(guī)模更小的方程組的求解,減少了潮流計算量。針對輸配網(wǎng)潮流計算的特點,文獻[5-10]提出基于主從分裂式的分布式潮流計算方法,將輸配網(wǎng)分為主系統(tǒng)和從系統(tǒng),主系統(tǒng)采用牛頓拉夫遜法,配電系統(tǒng)采用前推回代算法,加快了潮流收斂,減少了輸配網(wǎng)全局潮流計算量。文獻[11-13]對主從分裂法進行了改進,改善了輸配網(wǎng)在邊界點信息的匹配,進一步加快了潮流收斂。

    實際的輸配網(wǎng)全局電力系統(tǒng),配電系統(tǒng)節(jié)點數(shù)和支路數(shù)比相關(guān)的輸電系統(tǒng)的大一個甚至幾個數(shù)量級[5],因此,減少配電網(wǎng)潮流計算量是提升全局潮流計算速度的關(guān)鍵。為此,不少學(xué)者提出對配電網(wǎng)中的負荷與線路分支進行等效[14-15]的方法,極大地減少了計算量,但是直接進行配網(wǎng)模型的等效,會存在收斂性變差和計算結(jié)果不準確等問題。

    針對上述問題,本文提出一種考慮輸配網(wǎng)潮流相互影響的配電網(wǎng)等效擬合的潮流計算方法。通過兩次或三次主從迭代結(jié)果可對解配網(wǎng)潮流進行擬合,求解出配網(wǎng)根節(jié)點功率-電壓函數(shù)S(U)的擬合函數(shù)Sfit(U),由輸電網(wǎng)潮流計算得到的邊界點電壓U帶入擬合函數(shù)Sfit(U)進行配網(wǎng)等效擬合的輸電網(wǎng)潮流迭代,無需進行配電網(wǎng)子迭代?;谠摲椒?,通過傳統(tǒng)主從迭代與本文提出的配網(wǎng)等效擬合迭代交替進行,不斷修正擬合函數(shù)Sfit(U),縮小其與真實函數(shù)誤差,保證了潮流收斂的可靠性。同時,在很小的區(qū)間內(nèi)對S(U) 進行擬合誤差極小,潮流計算結(jié)果精確。本文方法極大地減少了配網(wǎng)潮流迭代次數(shù),加快了潮流收斂。

    本文研究的重點是所含節(jié)點為PQ節(jié)點類型的配網(wǎng)的等效擬合,對于部分配網(wǎng)中所含分布式電源可能是PI或PV等節(jié)點類型的情況,可采用常規(guī)方法進行配網(wǎng)的潮流迭代計算,本文暫不加以討論。

    1 配網(wǎng)等效擬合原理

    在配網(wǎng)潮流中,當所有節(jié)點的負荷大小確定時,根節(jié)點處的電壓幅值U與功率S呈非線性關(guān)系,假設(shè)該關(guān)系的函數(shù)表達式為S(U)。S(U)為復(fù)雜的非線性函數(shù),難以寫出其準確的表達式。但是可用多項式函數(shù)或其他形式函數(shù)對其進行擬合。因此全局潮流迭代計算時,可通過記錄前幾次各配網(wǎng)根節(jié)點處的電壓和功率作為離散點,從而對復(fù)雜函數(shù)S(U)進行擬合。

    進行曲線擬合最重要的是保證擬合精度,誤差太大則失去了擬合的意義。在某區(qū)間內(nèi),函數(shù)的極值越多,為保證擬合精度需采用更高階的多項式函數(shù)或者形式更復(fù)雜的函數(shù)進行擬合,此時參數(shù)難以求解,最理想的情況是需要被擬合的函數(shù)極值點個數(shù)為零,即為單調(diào)函數(shù)。因此需要分析函數(shù)S(U)的極值點個數(shù)以及單調(diào)性,以確定本文提出的設(shè)想是否具有實際意義以及可操作性。首先對配電網(wǎng)網(wǎng)絡(luò)模型進行簡化,等效電路圖如圖1所示。節(jié)點0為根節(jié)點,S0為根節(jié)點處總功率,節(jié)點1為與根節(jié)點相連的下一個節(jié)點,S1表示節(jié)點1處負荷與后續(xù)支路功率之和,Z1為支路阻抗。

    圖1 配電網(wǎng)等效電路圖Fig.1 Equivalent circuit diagram of distribution network

    由圖1可得根節(jié)點處電壓幅值與功率為

    (1)

    (2)

    由式(1)~(2)可得

    (3)

    (4)

    從節(jié)點0~1的數(shù)學(xué)模型可遞推得到從節(jié)點(i-1)~i的數(shù)學(xué)模型,進而由配電網(wǎng)的簡化模型可推廣到一般模型。假設(shè):h(Ui)= (Ui2-bi)-[(Ui2-bi)2-ci]1/2,fp(Ui)=arih(Ui),fq(Ui)=axih(Ui),Pi=Pil+Pis,Ui=Ui(U0),其中Pi為節(jié)點i處負荷功率Pil與后續(xù)支路功率Pis之和;Ui為節(jié)點i電壓幅值,fp(Ui)、fq(Ui)表示以節(jié)點i為末節(jié)點的支路有功損耗和無功損耗,由(4)可得出整個配網(wǎng)各節(jié)點電壓的前推關(guān)系:

    (5)

    由式(5)可知Ui=Ui(Ui-1)為單調(diào)遞增函數(shù),從而由遞推關(guān)系可知Ui=Ui(U0)為單調(diào)遞增函數(shù),fp(Ui)和fq(Ui)為單調(diào)遞減函數(shù),因此復(fù)合函數(shù)fp(Ui(U0))、fq(Ui(U0))也為單調(diào)遞減函數(shù)。對配電網(wǎng)所有節(jié)點列寫其有功功率方程為

    (6)

    類似的方法可得根節(jié)點處無功功率的數(shù)學(xué)模型,對式(6)進行整理可得配電網(wǎng)功率的數(shù)學(xué)模型為

    (7)

    式(7)可知,對于配電網(wǎng)的一般模型,根節(jié)點處功率與電壓幅值函數(shù)P0(U0)、Q0(U0)為單調(diào)遞減函數(shù)而非多極值函數(shù),單調(diào)函數(shù)可以用線性函數(shù)或者拋物線等函數(shù)進行擬合,且擬合函數(shù)的參數(shù)易于求解。因此對配電網(wǎng)潮流計算模型進行等效擬合的構(gòu)想具有理論上的可行性。

    2 配電網(wǎng)等效擬合的方法

    2.1 線性擬合

    由式(7)可知,根節(jié)點處功率為配網(wǎng)總負荷與線路損耗之和,總負荷功率不變,根節(jié)點處功率變化主要為電壓幅值變化引起的線路損耗的變化,在配電網(wǎng)正常運行情況下,根節(jié)點處電壓幅值的變化范圍在額定電壓附近,因此,根節(jié)點處功率變化不大,一定程度上可視為線性函數(shù)。因此復(fù)雜的非線性函數(shù)式(7)可用線性函數(shù)進行擬合,擬合形式如下:

    (8)

    式中:k1、k2和l1、l2為擬合參數(shù),通過兩次主從迭代之后求出U0(1)、P0(1)、Q0(1)、U0(2)、P0(2)、Q0(2)、從而求出擬合參數(shù)。第二次主從迭代后,邊界點的功率與電壓關(guān)系可直接由擬合函數(shù)近似計算出。

    2.2 非線性擬合

    采用線性函數(shù)擬合非線性函數(shù)誤差較大,在精度要求較高時,不一定能滿足要求。這時需要采用誤差更小的擬合方法,非線性函數(shù)擬合非線性函數(shù)具有小的擬合誤差。

    2.2.1拋物線擬合

    拋物線擬合非線性單調(diào)函數(shù),以其具有較小的擬合誤差,是最常用的一種擬合方式。拋物線擬合式(7)的形式如下:

    (9)

    式中:v1、v2、w1、w2和y1、y2為擬合參數(shù),參數(shù)求解方法和擬合方式與線性擬合類似。

    2.2.2冪函數(shù)函數(shù)擬合

    與線性擬合相比,拋物線擬合雖然具有較小的擬合誤差,但是拋物線函數(shù)有3個參數(shù),需要主從迭代3次后聯(lián)立3個方程求解。擬合函數(shù)中參數(shù)越少,擬合法的優(yōu)勢將越明顯。在輸配網(wǎng)全局潮流計算時,由于配電網(wǎng)節(jié)點數(shù)量遠多于輸電網(wǎng)節(jié)點數(shù)量,主要計算量為配網(wǎng)潮流的子迭代,每減少一次配網(wǎng)潮流子迭代,計算量將會大大減少。找到需要求解參數(shù)少、滿足精度要求的擬合函數(shù),是配網(wǎng)等效擬合發(fā)揮其優(yōu)勢的關(guān)鍵。因此可采用需要求解參數(shù)較少的冪函數(shù)擬合的方式擬合式(7),冪函數(shù)擬合的形式如下:

    (10)

    式中:d1、d2和t1、t2為擬合參數(shù),參數(shù)求解方法和擬合方式與線性擬合類似。

    2.2.3擬合方法減少計算量機理分析

    傳統(tǒng)主從分裂法在主從迭代計算的過程中,邊界點功率和電壓會不斷波動,因此需要反復(fù)進行輸、配網(wǎng)的潮流迭代,配網(wǎng)潮流迭代得到的功率是當前邊界點電壓下的精確結(jié)果,因此主從迭代可視為一種精確迭代。而通過配網(wǎng)等效擬合的方法求解配電網(wǎng)擬合功率,此時沒有進行配網(wǎng)潮流計算,邊界點功率為擬合值,可視為擬合迭代。

    本文所提的配網(wǎng)等效擬合的輸配潮流計算方法,在進行兩次或3次主從迭代之后,即可求解擬合函數(shù)對配網(wǎng)功率進行擬合,而無需通過配網(wǎng)迭代得到邊界點功率,減少了配網(wǎng)迭代次數(shù),在此基礎(chǔ)上通過精確迭代與擬合迭代交替進行,改善了采用主從分裂法每次都進行精確迭代計算量大的問題。

    圖2 迭代過程變化示意圖Fig.2 Schematic diagram of the iterative process of change

    如圖2所示,其中SBS為邊界點功率,UB為邊界點電壓,Sf(UB)為擬合函數(shù)。潮流迭代過程中邊界點電壓和功率的波動,實質(zhì)上是因為迭代過程中線損的波動,和負荷功率相比,線損占比很小,而線損的波動值和負荷功率相比則更小,因此邊界點電壓UB波動值很小,在極小的區(qū)間內(nèi)對邊界點電壓-功率函數(shù)S(UB)進行擬合誤差極小,并且可通過精確迭代不斷修正擬合函數(shù)Sf(UB)使擬合值更加準確。圖中虛線框中為潮流計算循環(huán)體,可見,與傳統(tǒng)主從分裂法相比,本文所提的擬合迭代的方法減少了輸配網(wǎng)信息交互,直接通過擬合函數(shù)Sf(UB)求解邊界點信息,無需每次邊界點功率傳遞都進行配電網(wǎng)潮流計算,增加了少量輸電網(wǎng)迭代次數(shù),減少了配網(wǎng)迭代次數(shù),由于配電網(wǎng)規(guī)模遠大于輸電網(wǎng),輸電網(wǎng)潮流計算量可忽略,因此本文方法能夠減少輸配網(wǎng)潮流計算量,加快潮流收斂。

    3 基于配網(wǎng)等效擬合協(xié)調(diào)潮流算法

    3.1 輸配網(wǎng)分解協(xié)調(diào)潮流模型

    主從分裂法將全局電網(wǎng)分解為輸電網(wǎng)系統(tǒng)、邊界系統(tǒng)和配電網(wǎng)系統(tǒng)[2]。其中,輸電網(wǎng)為主系統(tǒng),配電網(wǎng)為從系統(tǒng),邊界系統(tǒng)是由輸電網(wǎng)中與配電網(wǎng)相連的節(jié)點組成,在邊界系統(tǒng)進行輸電網(wǎng)與配電網(wǎng)信息的協(xié)調(diào),從而實現(xiàn)輸配網(wǎng)潮流的分解協(xié)調(diào)計算。輸配網(wǎng)全局模型簡化圖如圖3所示。

    圖3 輸配網(wǎng)全局模型簡化圖Fig.3 Simplified map of the global model of the transmission and distribution network

    其中,輸電網(wǎng)潮流計算數(shù)學(xué)模型為

    i∈CM,j∈(CM∪CB)

    (11)

    邊界點潮流計算的數(shù)學(xué)模型為

    i∈CB,j∈(CM∪CB),j∈(CD)

    (12)

    配電網(wǎng)潮流計算數(shù)學(xué)模型為

    i∈CD,j∈(CD∪CT)

    (13)

    式中:CM、CB分別為輸電網(wǎng)、邊界節(jié)點集合;CD為第i個邊界點所接配電網(wǎng)節(jié)點集合;CT為第i個邊界點;Si為節(jié)點i的注入功率。

    輸配網(wǎng)全局潮流計算的主要計算量是對數(shù)量眾多的配電網(wǎng)進行潮流計算,配網(wǎng)潮流等效擬合實質(zhì)是對配電網(wǎng)潮流計算模型進行簡化,將配電網(wǎng)模型與輸電網(wǎng)中的邊界系統(tǒng)模型進行合并,主從模型合并后,無需進行配網(wǎng)潮流子迭代,大大減少潮流計算量。配電網(wǎng)等效擬合后邊界點潮流計算的數(shù)學(xué)模型為

    i∈CB,j∈(CM∪CB)

    (14)

    式中:Ui0為第i個配電網(wǎng)根節(jié)點處電壓幅值;fis(Ui0)為第i個配電網(wǎng)潮流的擬合函數(shù)。

    傳統(tǒng)的主從分裂模型由式(11)~(13)構(gòu)成;基于配網(wǎng)潮流等效擬合的潮流計算模型由式(11)和(14)構(gòu)成,在邊界點處對輸電網(wǎng)模型和配電網(wǎng)模型進行合并,構(gòu)成主從合并的潮流計算模型。通過主從分裂和主從合并兩種潮流計算模型交替使用,進行輸配網(wǎng)全局潮流計算。

    本文輸電網(wǎng)潮流計算采用單相模型,配網(wǎng)潮流計算采用三相模型,對于配電系統(tǒng)三相不平衡問題的處理,本文采用文獻[5]所提方法進行邊界點功率、電壓交互,記錄每次主從迭代邊界點處功率和電壓結(jié)果,求解或修正擬合函數(shù),可根據(jù)實際需要對配電網(wǎng)ABC三相分別進行擬合或整體進行擬合,因此,對于配電網(wǎng)可能存在三相不平衡的情況,不影響本文所提配網(wǎng)等效擬合方法的使用。

    3.2 潮流計算步驟

    基于配電網(wǎng)等效擬合的輸配網(wǎng)全局潮流計算的計算步驟如下(假設(shè)m為擬合函數(shù)參數(shù)個數(shù)):

    ① 給輸配網(wǎng)所有節(jié)點電壓賦初值,k=0;

    4 算例分析

    用matlab編寫基于牛頓-拉夫遜法的輸電網(wǎng)潮流、基于前推回代法的配電網(wǎng)潮流和基于配地網(wǎng)等效擬合的協(xié)調(diào)等3個模塊,用于模擬輸配網(wǎng)的全局潮流計算。本文以IEEE118節(jié)點系統(tǒng)為輸電網(wǎng)系統(tǒng),在其節(jié)點1至節(jié)點60均接入290節(jié)點配電網(wǎng)作為配電網(wǎng)系統(tǒng),其中輸電網(wǎng)電壓等級為100kV,配網(wǎng)電壓等級為10kV。分別用本文所提主從迭代與擬合迭代交替進行的方法和文獻[6]所提的單一主從迭代不擬合的方法,對以上算例進行輸配網(wǎng)全局潮流計算(收斂精度為10-9,功率基準值取100MVA),并對計算結(jié)果和計算時間進行對比分析。

    表1給出了分別采用本文所提的3種等效擬合方法和文獻[6]所提單一主從迭代不擬合的方法進行輸配網(wǎng)全局潮流計算時,潮流收斂后計算結(jié)果對比。其中max|ΔU|與max|ΔSB|分別表示幾種擬合方式進行潮流計算時,相應(yīng)節(jié)點電壓差值的最大值和邊界點功率差值的最大值。

    由表1可知,在輸配網(wǎng)全局潮流計算中,采用擬合方法的計算結(jié)果與不擬合相比,邊界點電壓和功率的最大誤差很小,電壓最大之誤差為1.33×10-10,邊界點功率最大誤差為5.05×10-10,均小于10-9,本文所提方法進行潮流計算誤差很小,不存在邊界點電壓和功率的失配問題。

    表1 擬合與不擬合潮流計算結(jié)果對比

    表2為采用擬合與不擬合的輸配網(wǎng)潮流計算迭代次數(shù)對比。由表2可知,與傳統(tǒng)的單一主從迭代配網(wǎng)不擬合的方法相比,本文所提3種擬合方式潮流計算均減少了主從迭代次數(shù),不擬合的情況下主從迭代次數(shù)為10次,配網(wǎng)子迭代為11次,采用本文所提的線性擬合和冪函數(shù)擬合方法主從迭代次數(shù)為3次,配網(wǎng)子迭代次數(shù)為4次,主從迭代和配網(wǎng)子迭代次數(shù)均減少了7次。本文方法將部分主從迭代步驟采用配網(wǎng)等效擬合的輸電網(wǎng)迭代代替,因此輸電網(wǎng)迭代次數(shù)有少量增加,但是輸電網(wǎng)規(guī)模遠小于配網(wǎng)規(guī)模,增加的輸電網(wǎng)計算量很小可忽略,輸配網(wǎng)潮流計算量主要由主從迭代或配網(wǎng)迭代次數(shù)決定,本文方法能夠大大減少主從迭代和配網(wǎng)迭代次數(shù),因此,能夠顯著減少潮流計算量。

    表2 擬合與不擬合潮流計算迭代次數(shù)對比

    表3為擬合與不擬合方式全局潮流計算時間對比,其中節(jié)約時間是以不擬合潮流計算時間為基準,擬合方式節(jié)約的時間百分數(shù),計算結(jié)果為同一臺計算機仿真所得。由表3可知,在輸配網(wǎng)全局潮流計算中,由于配網(wǎng)的節(jié)點數(shù)遠超過輸電網(wǎng)節(jié)點數(shù),因此配電網(wǎng)潮流計算占了輸配網(wǎng)全局潮流計算的絕大多數(shù)時間,結(jié)合表2可知,配網(wǎng)迭代次數(shù)是影響潮流計算時間的關(guān)鍵。采用本文方法能夠減少配網(wǎng)迭代次數(shù),加快潮流收斂,因為線性函數(shù)與冪函數(shù)擬合時只需主從迭代兩次即可進行擬合迭代,而拋物線需要進行3次才能求解出擬合參數(shù),因此采用兩個參數(shù)進行擬合減少計算量的效果更佳,其中,線性擬合與冪函數(shù)擬合減少潮流計算時間約64%,拋物線擬合減少時間約為54.8%,可見,本文方法在保證潮流計算精度的情況下,計算速度顯著提升。

    表3 擬合與不擬合潮流計算時間對比

    5 結(jié) 論

    ① 本文首次提出了配電網(wǎng)等效擬合的概念,并據(jù)此構(gòu)建了基于等效擬合的輸配網(wǎng)潮流計算方法,并給出了線性函數(shù)與非線性函數(shù)擬合的形式。

    ② 在對配電網(wǎng)進行等效擬合的基礎(chǔ)上構(gòu)建了主從分裂與主從合并交替使用的潮流計算模型,通過主從分裂的精確迭代與配電等效的擬合迭代交替進行,可大幅減少主從迭代次數(shù)。經(jīng)仿真驗證,該方法計算結(jié)果精度可靠,計算速度較單一主從分裂法有明顯提升。

    ③ 本文所提方法在提高輸配網(wǎng)全局潮流計算速度上邁出了重要一步,可作為研究輸配網(wǎng)協(xié)同框架內(nèi)各類問題的基礎(chǔ),如研究輸配網(wǎng)電壓合格率提升措施,尋找輸配網(wǎng)最優(yōu)的電壓無功協(xié)調(diào)優(yōu)化控制策略,也可用于輸配網(wǎng)網(wǎng)架重構(gòu),加快網(wǎng)架重構(gòu)速度,提高系統(tǒng)運行的安全性和經(jīng)濟性,也可用于輸配網(wǎng)協(xié)同的配電網(wǎng)態(tài)勢感知,評估配電網(wǎng)靜態(tài)電壓穩(wěn)定態(tài)勢,也可用于輸配網(wǎng)故障分析、調(diào)度控制等。綜上所述,本文所提方法能夠提高輸配網(wǎng)系統(tǒng)的反應(yīng)速度,在輸配協(xié)同電網(wǎng)安全運行、預(yù)防控制和協(xié)調(diào)優(yōu)化等方面具有重要意義。

    猜你喜歡
    輸電網(wǎng)邊界點主從
    道路空間特征與測量距離相結(jié)合的LiDAR道路邊界點提取算法
    層次化點云邊界快速精確提取方法研究
    FANUC系統(tǒng)PROFIBUS主從功能應(yīng)用
    基于主從控制的微電網(wǎng)平滑切換控制策略研究
    計及多重不確定因素的輸電網(wǎng)隨機潮流計算
    含光伏電站的輸電網(wǎng)不對稱故障分析方法
    電測與儀表(2016年4期)2016-04-11 09:44:34
    基于飛行試驗數(shù)據(jù)的仿真模型主從一體化檢驗
    基于差分和聲搜索算法的輸電網(wǎng)差異化規(guī)劃
    一種去除掛網(wǎng)圖像鋸齒的方法及裝置
    電腦與電信(2014年6期)2014-03-22 13:21:06
    基于FPGA的機床數(shù)控系統(tǒng)主從控制器通信研究
    可以在线观看的亚洲视频| 好看av亚洲va欧美ⅴa在| 一区二区三区激情视频| 九九久久精品国产亚洲av麻豆 | 最近在线观看免费完整版| 法律面前人人平等表现在哪些方面| 韩国av一区二区三区四区| 精品国产超薄肉色丝袜足j| 日韩 欧美 亚洲 中文字幕| 叶爱在线成人免费视频播放| 亚洲人成伊人成综合网2020| 久久精品国产综合久久久| 老汉色∧v一级毛片| 婷婷六月久久综合丁香| 成人午夜高清在线视频| 久久久国产成人精品二区| 99国产精品99久久久久| 黄色 视频免费看| 丰满人妻一区二区三区视频av | 在线观看午夜福利视频| 久久午夜综合久久蜜桃| 老司机在亚洲福利影院| 1024香蕉在线观看| 青草久久国产| 久久午夜综合久久蜜桃| 我的老师免费观看完整版| 久久中文看片网| 99久久无色码亚洲精品果冻| 亚洲五月婷婷丁香| 淫妇啪啪啪对白视频| 麻豆成人av在线观看| 在线观看一区二区三区| 最近最新免费中文字幕在线| 国产 一区 欧美 日韩| 亚洲专区国产一区二区| 久久久成人免费电影| 国产成人精品久久二区二区91| 日本黄色片子视频| 亚洲va日本ⅴa欧美va伊人久久| 男女午夜视频在线观看| 日日干狠狠操夜夜爽| 成年女人毛片免费观看观看9| 亚洲国产欧美网| 18禁观看日本| 舔av片在线| 精品一区二区三区四区五区乱码| 亚洲无线观看免费| 国产av在哪里看| 亚洲无线在线观看| 女人被狂操c到高潮| 亚洲中文日韩欧美视频| 亚洲自拍偷在线| 亚洲狠狠婷婷综合久久图片| 又爽又黄无遮挡网站| 欧美日韩国产亚洲二区| 久久人人精品亚洲av| 露出奶头的视频| 欧美乱妇无乱码| 欧美大码av| 亚洲成人免费电影在线观看| 亚洲精品粉嫩美女一区| 夜夜夜夜夜久久久久| 黑人操中国人逼视频| 小说图片视频综合网站| 欧美又色又爽又黄视频| 久久久久久九九精品二区国产| 亚洲精品在线观看二区| 巨乳人妻的诱惑在线观看| 国产精品免费一区二区三区在线| 国产精品九九99| 51午夜福利影视在线观看| 夜夜爽天天搞| 真人一进一出gif抽搐免费| x7x7x7水蜜桃| 午夜成年电影在线免费观看| 久久国产精品人妻蜜桃| 欧美最黄视频在线播放免费| 中文资源天堂在线| 精品国产亚洲在线| 免费在线观看影片大全网站| 日韩欧美国产一区二区入口| 草草在线视频免费看| 久久久精品欧美日韩精品| 国产欧美日韩一区二区精品| 欧美乱码精品一区二区三区| 精品国产三级普通话版| 两性夫妻黄色片| 一本精品99久久精品77| 亚洲电影在线观看av| 久久久久久久久久黄片| 热99re8久久精品国产| 国产亚洲精品综合一区在线观看| 一级a爱片免费观看的视频| а√天堂www在线а√下载| 人妻久久中文字幕网| 欧美不卡视频在线免费观看| 岛国视频午夜一区免费看| 好男人在线观看高清免费视频| 免费在线观看日本一区| 最好的美女福利视频网| 亚洲七黄色美女视频| 熟女少妇亚洲综合色aaa.| 激情在线观看视频在线高清| avwww免费| 夜夜爽天天搞| 国产精品电影一区二区三区| 美女免费视频网站| 色精品久久人妻99蜜桃| 啦啦啦观看免费观看视频高清| 18禁观看日本| 亚洲aⅴ乱码一区二区在线播放| 国产精品 欧美亚洲| 精品不卡国产一区二区三区| 亚洲无线在线观看| 午夜成年电影在线免费观看| 国产爱豆传媒在线观看| 黄色视频,在线免费观看| 首页视频小说图片口味搜索| 亚洲国产精品合色在线| 亚洲,欧美精品.| 午夜精品久久久久久毛片777| 亚洲中文日韩欧美视频| 久久精品aⅴ一区二区三区四区| 这个男人来自地球电影免费观看| 一级毛片高清免费大全| 一a级毛片在线观看| 一边摸一边抽搐一进一小说| 夜夜躁狠狠躁天天躁| 日本免费一区二区三区高清不卡| 亚洲精品中文字幕一二三四区| 国产精品永久免费网站| 国产乱人视频| 亚洲国产欧美一区二区综合| 日日夜夜操网爽| 欧美又色又爽又黄视频| 欧美激情在线99| 国产综合懂色| 国产男靠女视频免费网站| av视频在线观看入口| 高清在线国产一区| 国产精品影院久久| 十八禁人妻一区二区| 免费观看精品视频网站| 亚洲国产看品久久| 午夜精品在线福利| 黑人操中国人逼视频| 亚洲国产看品久久| 18禁黄网站禁片午夜丰满| 欧美乱码精品一区二区三区| 最新中文字幕久久久久 | 在线观看舔阴道视频| 午夜福利欧美成人| 别揉我奶头~嗯~啊~动态视频| 午夜日韩欧美国产| 国产一区在线观看成人免费| 99精品久久久久人妻精品| 国产亚洲精品av在线| www日本黄色视频网| 免费看十八禁软件| 国产伦一二天堂av在线观看| 精品乱码久久久久久99久播| 狂野欧美激情性xxxx| 在线a可以看的网站| 天堂√8在线中文| 色视频www国产| 黄色成人免费大全| 国产av不卡久久| av片东京热男人的天堂| 精品国产超薄肉色丝袜足j| 九九热线精品视视频播放| 精品久久蜜臀av无| 欧美一区二区精品小视频在线| 制服丝袜大香蕉在线| 亚洲国产中文字幕在线视频| 亚洲中文字幕日韩| 18禁美女被吸乳视频| 亚洲欧洲精品一区二区精品久久久| 欧美一区二区精品小视频在线| 男女之事视频高清在线观看| 观看免费一级毛片| 搡老岳熟女国产| 性色avwww在线观看| 少妇的逼水好多| 国产一区二区在线观看日韩 | 女警被强在线播放| 变态另类成人亚洲欧美熟女| 亚洲狠狠婷婷综合久久图片| 国产精品一区二区免费欧美| 国产成人系列免费观看| 久99久视频精品免费| 欧美精品啪啪一区二区三区| 长腿黑丝高跟| 在线看三级毛片| 欧美绝顶高潮抽搐喷水| АⅤ资源中文在线天堂| 国产久久久一区二区三区| 免费在线观看影片大全网站| 精品久久久久久久末码| 看免费av毛片| 午夜福利欧美成人| 别揉我奶头~嗯~啊~动态视频| 老司机深夜福利视频在线观看| 在线永久观看黄色视频| 色视频www国产| 亚洲成人免费电影在线观看| 91在线精品国自产拍蜜月 | 动漫黄色视频在线观看| 日本免费a在线| 日本一本二区三区精品| 五月伊人婷婷丁香| 18禁黄网站禁片午夜丰满| 老司机在亚洲福利影院| 俄罗斯特黄特色一大片| 夜夜躁狠狠躁天天躁| 国产成人系列免费观看| 久久这里只有精品19| 黑人巨大精品欧美一区二区mp4| or卡值多少钱| 久久久精品大字幕| av欧美777| 国产美女午夜福利| 国产精品久久久久久久电影 | 日本黄色片子视频| 999久久久国产精品视频| 久久亚洲真实| 亚洲一区二区三区不卡视频| АⅤ资源中文在线天堂| 亚洲天堂国产精品一区在线| 亚洲成人中文字幕在线播放| 三级毛片av免费| 在线a可以看的网站| 亚洲狠狠婷婷综合久久图片| 怎么达到女性高潮| 色综合亚洲欧美另类图片| 97超级碰碰碰精品色视频在线观看| 亚洲人成网站在线播放欧美日韩| 日日干狠狠操夜夜爽| 一级毛片精品| 免费电影在线观看免费观看| 亚洲国产中文字幕在线视频| 少妇人妻一区二区三区视频| 成人av在线播放网站| 成年版毛片免费区| 黄频高清免费视频| 亚洲午夜理论影院| 宅男免费午夜| 99久国产av精品| 一区二区三区高清视频在线| 丝袜人妻中文字幕| 曰老女人黄片| 日本a在线网址| 看片在线看免费视频| 亚洲美女视频黄频| 亚洲第一欧美日韩一区二区三区| 亚洲真实伦在线观看| 亚洲五月婷婷丁香| 小说图片视频综合网站| 免费人成视频x8x8入口观看| 搡老熟女国产l中国老女人| 美女黄网站色视频| 国产av一区在线观看免费| 免费看a级黄色片| 一本一本综合久久| 日本成人三级电影网站| 欧美zozozo另类| 国产成人精品久久二区二区免费| 88av欧美| 国产高清videossex| 亚洲精品一区av在线观看| 欧美乱妇无乱码| 网址你懂的国产日韩在线| 舔av片在线| 日韩大尺度精品在线看网址| 国产高清激情床上av| 在线看三级毛片| 1000部很黄的大片| 青草久久国产| 成人特级黄色片久久久久久久| 两个人视频免费观看高清| 成人精品一区二区免费| АⅤ资源中文在线天堂| 中文字幕最新亚洲高清| 99久久精品热视频| 亚洲精品乱码久久久v下载方式 | 少妇裸体淫交视频免费看高清| 99精品在免费线老司机午夜| 97超视频在线观看视频| 蜜桃久久精品国产亚洲av| 日本与韩国留学比较| tocl精华| 成人特级黄色片久久久久久久| 亚洲aⅴ乱码一区二区在线播放| 91在线观看av| 舔av片在线| 窝窝影院91人妻| 黄片小视频在线播放| 婷婷丁香在线五月| 少妇人妻一区二区三区视频| 亚洲一区高清亚洲精品| 日本五十路高清| 好男人在线观看高清免费视频| av天堂中文字幕网| 欧美激情在线99| 国产亚洲精品一区二区www| 天堂√8在线中文| 亚洲天堂国产精品一区在线| 精品国产美女av久久久久小说| 首页视频小说图片口味搜索| 亚洲美女黄片视频| 99久久精品国产亚洲精品| 成人av一区二区三区在线看| 亚洲欧美日韩高清专用| 不卡一级毛片| 亚洲色图 男人天堂 中文字幕| 麻豆成人av在线观看| 国产成人精品无人区| 久久精品亚洲精品国产色婷小说| 亚洲九九香蕉| 欧美日本视频| 中文字幕熟女人妻在线| 黑人操中国人逼视频| 久久香蕉国产精品| 国产激情欧美一区二区| 国产一区二区激情短视频| 精品一区二区三区av网在线观看| 国产高清视频在线观看网站| 亚洲av成人精品一区久久| 丰满的人妻完整版| 亚洲天堂国产精品一区在线| 丰满人妻一区二区三区视频av | 免费在线观看日本一区| 男女床上黄色一级片免费看| 高潮久久久久久久久久久不卡| 国产精品一及| 国产精品一区二区免费欧美| 搡老妇女老女人老熟妇| www日本在线高清视频| 国产精品自产拍在线观看55亚洲| 国产亚洲精品综合一区在线观看| 成在线人永久免费视频| 中文字幕av在线有码专区| 国产一区二区在线av高清观看| 真实男女啪啪啪动态图| 五月玫瑰六月丁香| 国产探花在线观看一区二区| 亚洲精品国产精品久久久不卡| 国产 一区 欧美 日韩| 国产欧美日韩一区二区三| 国产精品 欧美亚洲| av欧美777| 亚洲av电影不卡..在线观看| 热99在线观看视频| 中文字幕最新亚洲高清| 亚洲精品色激情综合| 亚洲第一电影网av| 色吧在线观看| 日本一本二区三区精品| 老汉色av国产亚洲站长工具| 看免费av毛片| 欧美一区二区国产精品久久精品| netflix在线观看网站| 成人特级av手机在线观看| 免费一级毛片在线播放高清视频| 成人特级av手机在线观看| 天堂av国产一区二区熟女人妻| 国产一区二区在线av高清观看| 看免费av毛片| 九九热线精品视视频播放| 99国产精品99久久久久| 极品教师在线免费播放| 免费看光身美女| 日日干狠狠操夜夜爽| 91在线观看av| 男女之事视频高清在线观看| 曰老女人黄片| 丝袜人妻中文字幕| 成人特级黄色片久久久久久久| 日本一二三区视频观看| 久久香蕉精品热| 午夜福利在线观看吧| 在线a可以看的网站| 又大又爽又粗| 成年女人看的毛片在线观看| 好男人在线观看高清免费视频| 麻豆久久精品国产亚洲av| 久久久色成人| 国内少妇人妻偷人精品xxx网站 | 男人舔女人下体高潮全视频| 99久久成人亚洲精品观看| 草草在线视频免费看| 国内久久婷婷六月综合欲色啪| 国产伦精品一区二区三区视频9 | 人人妻人人看人人澡| 露出奶头的视频| 欧美乱码精品一区二区三区| 欧美成人性av电影在线观看| 精品福利观看| 嫁个100分男人电影在线观看| 国产精品久久久久久久电影 | 欧美一级a爱片免费观看看| 五月玫瑰六月丁香| 一级毛片高清免费大全| 精品国产乱码久久久久久男人| 亚洲自偷自拍图片 自拍| 18禁美女被吸乳视频| 岛国在线观看网站| 欧美成人一区二区免费高清观看 | 在线免费观看不下载黄p国产 | 精品日产1卡2卡| 美女午夜性视频免费| 国产精品久久电影中文字幕| 18禁国产床啪视频网站| 窝窝影院91人妻| 欧美黑人欧美精品刺激| svipshipincom国产片| 搡老岳熟女国产| 蜜桃久久精品国产亚洲av| 精品人妻1区二区| 亚洲欧美日韩高清专用| 大型黄色视频在线免费观看| 国产精品亚洲av一区麻豆| av国产免费在线观看| 法律面前人人平等表现在哪些方面| 国产伦精品一区二区三区视频9 | 国产亚洲精品av在线| 午夜亚洲福利在线播放| av天堂中文字幕网| 麻豆国产av国片精品| 亚洲专区中文字幕在线| 成人鲁丝片一二三区免费| 国产不卡一卡二| av黄色大香蕉| 麻豆国产97在线/欧美| 99国产精品99久久久久| 性色av乱码一区二区三区2| 久久久久国产精品人妻aⅴ院| 最近最新免费中文字幕在线| 日本黄大片高清| 亚洲人与动物交配视频| 在线视频色国产色| 亚洲欧洲精品一区二区精品久久久| 国产成年人精品一区二区| 一级毛片精品| 性色av乱码一区二区三区2| 网址你懂的国产日韩在线| 桃红色精品国产亚洲av| 好男人电影高清在线观看| 欧美丝袜亚洲另类 | 日本在线视频免费播放| 男女床上黄色一级片免费看| 亚洲国产欧美网| а√天堂www在线а√下载| 一二三四在线观看免费中文在| 给我免费播放毛片高清在线观看| 波多野结衣高清无吗| 亚洲第一欧美日韩一区二区三区| 香蕉av资源在线| 欧美一级a爱片免费观看看| 狂野欧美激情性xxxx| 国产又色又爽无遮挡免费看| 日本一本二区三区精品| 精品久久久久久久毛片微露脸| 在线十欧美十亚洲十日本专区| 夜夜爽天天搞| 1000部很黄的大片| 国产精品亚洲av一区麻豆| 草草在线视频免费看| 国内精品久久久久精免费| 舔av片在线| 曰老女人黄片| 最近在线观看免费完整版| 日本五十路高清| 99热这里只有精品一区 | 一卡2卡三卡四卡精品乱码亚洲| 国产免费av片在线观看野外av| 国产精品一区二区免费欧美| 噜噜噜噜噜久久久久久91| 国内精品一区二区在线观看| 一本久久中文字幕| 亚洲精品在线观看二区| 亚洲国产精品久久男人天堂| 国产成人福利小说| 在线十欧美十亚洲十日本专区| 99riav亚洲国产免费| 欧美一区二区精品小视频在线| 国产精品亚洲av一区麻豆| 国产精品久久久人人做人人爽| 成人一区二区视频在线观看| 精品免费久久久久久久清纯| 亚洲 欧美 日韩 在线 免费| 99热6这里只有精品| 五月伊人婷婷丁香| 日本一本二区三区精品| 国产精品av久久久久免费| 欧美不卡视频在线免费观看| 麻豆国产av国片精品| 亚洲成人精品中文字幕电影| 国产又黄又爽又无遮挡在线| 天天添夜夜摸| 国产黄片美女视频| 久99久视频精品免费| 一级a爱片免费观看的视频| 观看免费一级毛片| 国产精品精品国产色婷婷| 亚洲欧美精品综合一区二区三区| 亚洲成人久久爱视频| 欧美在线黄色| 日韩大尺度精品在线看网址| 18禁美女被吸乳视频| 一区二区三区激情视频| e午夜精品久久久久久久| 亚洲美女视频黄频| 精品久久久久久久人妻蜜臀av| 国产成人精品无人区| 身体一侧抽搐| 又爽又黄无遮挡网站| 黄片小视频在线播放| netflix在线观看网站| 最近最新免费中文字幕在线| 伊人久久大香线蕉亚洲五| 人妻久久中文字幕网| 一二三四社区在线视频社区8| 特级一级黄色大片| 久久性视频一级片| 亚洲国产精品sss在线观看| 久久热在线av| 一个人免费在线观看电影 | aaaaa片日本免费| 国产免费av片在线观看野外av| 天堂av国产一区二区熟女人妻| 69av精品久久久久久| 免费在线观看影片大全网站| 人妻久久中文字幕网| 亚洲熟妇熟女久久| 亚洲午夜精品一区,二区,三区| 亚洲色图 男人天堂 中文字幕| 少妇人妻一区二区三区视频| 亚洲在线观看片| 18禁美女被吸乳视频| 久久精品国产亚洲av香蕉五月| 可以在线观看毛片的网站| 久久久久久久精品吃奶| 精品电影一区二区在线| 亚洲国产精品久久男人天堂| 久久久国产精品麻豆| 变态另类成人亚洲欧美熟女| 热99re8久久精品国产| 久久99热这里只有精品18| 色综合婷婷激情| 亚洲精品久久国产高清桃花| 婷婷精品国产亚洲av在线| 伦理电影免费视频| 久久久精品欧美日韩精品| 99视频精品全部免费 在线 | 欧美午夜高清在线| 久久久久免费精品人妻一区二区| 十八禁人妻一区二区| 亚洲美女黄片视频| 欧美一级a爱片免费观看看| 亚洲七黄色美女视频| 国产精品久久久久久久电影 | 一进一出抽搐动态| 网址你懂的国产日韩在线| 俺也久久电影网| 亚洲乱码一区二区免费版| 精品久久久久久久久久免费视频| 亚洲av第一区精品v没综合| 每晚都被弄得嗷嗷叫到高潮| 欧美日韩亚洲国产一区二区在线观看| 亚洲国产欧美人成| 在线观看66精品国产| 两人在一起打扑克的视频| 99久久精品热视频| 一级毛片高清免费大全| 狂野欧美白嫩少妇大欣赏| 两性午夜刺激爽爽歪歪视频在线观看| 青草久久国产| 国产一区二区在线观看日韩 | 好看av亚洲va欧美ⅴa在| 久久久久九九精品影院| 美女大奶头视频| 窝窝影院91人妻| 国产精品久久视频播放| 亚洲乱码一区二区免费版| 中文字幕久久专区| 嫩草影院入口| 午夜影院日韩av| 啦啦啦观看免费观看视频高清| 色播亚洲综合网| av天堂在线播放| 国产高清三级在线| 精品免费久久久久久久清纯| 国产爱豆传媒在线观看| 天堂影院成人在线观看| 少妇裸体淫交视频免费看高清| 免费在线观看影片大全网站| 亚洲国产精品久久男人天堂| 成在线人永久免费视频| 亚洲一区二区三区色噜噜| 欧美中文综合在线视频| 极品教师在线免费播放| 1000部很黄的大片| 成人18禁在线播放| 少妇裸体淫交视频免费看高清| 午夜视频精品福利| 中文字幕最新亚洲高清| 91在线精品国自产拍蜜月 | 久久久久久大精品| 国产亚洲av高清不卡| 91麻豆av在线| 欧美性猛交黑人性爽| 亚洲精品在线美女| 日韩欧美在线二视频| 我要搜黄色片| 日本免费一区二区三区高清不卡|