《電子技術(shù)應(yīng)用》
您所在的位置:首頁 > 電源技術(shù) > 設(shè)計應(yīng)用 > 輸電斷面的多角度認識與搜索方法研究
輸電斷面的多角度認識與搜索方法研究
2015年電子技術(shù)應(yīng)用第12期
楊萬清,尚亞男,姜學樸,張葆剛,劉 林,劉 冰
國網(wǎng)大連供電公司,遼寧 大連116024
摘要: 輸電斷面一直是電力系統(tǒng)安全分析和事故應(yīng)急處理的重要監(jiān)控和分析控制對象。從潮流轉(zhuǎn)移斷面、輸電通道以及電磁環(huán)網(wǎng)三種類型進行輸電斷面的多角度認識分析工作,并在此基礎(chǔ)上進行搜索方法研究。首先闡述了潮流轉(zhuǎn)移因子概念以及基于此定義的潮流轉(zhuǎn)移斷面,并提出利用割點和前k路徑搜索潮流轉(zhuǎn)移斷面;其次提出聯(lián)絡(luò)線的確定方法,并基于聯(lián)絡(luò)線開展節(jié)點收縮分區(qū)研究,最終根據(jù)潮流同向性原則確定輸電通道;最后指出電磁環(huán)網(wǎng)作為輸電通道的一個特例,能夠表明輸電通道的多電壓耦合關(guān)系。以新英格蘭39和IEEE118節(jié)點系統(tǒng)進行的輸電斷面搜索方法實例研究證明了本文搜索策略的準確性和有效性。
中圖分類號: TM744
文獻標識碼: A
DOI:10.16157/j.issn.0258-7998.2015.12.039

中文引用格式: 楊萬清,尚亞男,姜學樸,等. 輸電斷面的多角度認識與搜索方法研究[J].電子技術(shù)應(yīng)用,2015,41(12):148-151,155.
英文引用格式: Yang Wanqing,Shang Yanan,Jiang Xuepu,et al. Multi-angle understanding of transmission section and the research on searching method[J].Application of Electronic Technique,2015,41(12):148-151,155.
Multi-angle understanding of transmission section and the research on searching method
Yang Wanqing,Shang Yanan,Jiang Xuepu,Zhang Baogang,Liu Lin,Liu Bing
Dalian Power Supply Company,Dalian 116024,China
Abstract: The transmission section has been regarded as the important object to be monitored, analyzed and controlled in the power system security analysis and emergency response for a long time. This article presents a deep and complete transmission section analysis from three different perspectives of power flow transferring section, transmission corridor and electromagnetic loop, then proposes the searching method of transmission section.The concept of flow transferring relativity factor and the power flow transferring section based on this concept is described, then proposes the searching method based on cut-vertex and K shortest patchs; It proposes the method for determining tie-line and the node contraction partition. With the principle of consistency of power flow, transmission corridor is determined; It also points out that electromagnetic loop network, as a special case, can show the multi-voltage coupling relations among transmission corridors. The New England 39-node system and IEEE 118-node system is used as examples to illustrate the accuracy and efficiency of the proposed searching method.
Key words : transmisssion section;searching method;power flow transferring section;transmission corridor;electromagnetic loop

    

0 引言

    輸電斷面作為電網(wǎng)安全上關(guān)聯(lián)緊密的一組輸電元件集合,一直是電力系統(tǒng)安全分析和事故應(yīng)急處理的重要監(jiān)控和分析控制對象[1-2]

    現(xiàn)有文獻對于輸電斷面的認識與搜索方法研究主要從兩方面展開:(1)分析輸電網(wǎng)絡(luò)某一支路過載切除后,原有潮流的轉(zhuǎn)移分配模式,尋找其余回路中潮流變化較大的支路形成輸電斷面[3-7]。(2)在廠站作為單元節(jié)點的輸電網(wǎng)絡(luò)中結(jié)合電氣距離進行節(jié)點收縮與分區(qū),通過割集搜索和潮流同向性判斷形成輸電斷面[8-9]

    本文在上述文獻基礎(chǔ)上進一步深化研究,從不同角度全面認識評價輸電斷面,建立輸電斷面認識體系。綜合電網(wǎng)實際運行模式和已有研究成果,本文認為輸電斷面認識體系應(yīng)包含三類研究對象:以網(wǎng)絡(luò)潮流轉(zhuǎn)移變化為特征的潮流轉(zhuǎn)移斷面;以聯(lián)絡(luò)線潮流輸送同向性為特征的輸電通道;以多電壓等級線路拓撲耦合為特征的電磁環(huán)網(wǎng)

    輸電斷面的三種研究對象都立足于電網(wǎng)安全分析。潮流轉(zhuǎn)移斷面關(guān)注局部網(wǎng)絡(luò)的潮流變化情況,輸電斷面與切除支路密切相關(guān);輸電通道關(guān)注區(qū)域間潮流輸送情況,輸電斷面與分區(qū)聯(lián)絡(luò)線密切相關(guān);電磁環(huán)網(wǎng)關(guān)注線路間的拓撲耦合關(guān)系,輸電斷面與環(huán)網(wǎng)拓撲和電壓層級密切相關(guān)。

    在輸電斷面的認識評價基礎(chǔ)上,通過圖論分析,本文設(shè)計了利用割點和前k路徑的潮流轉(zhuǎn)移斷面搜索方法以及基于聯(lián)絡(luò)線節(jié)點收縮分區(qū)和潮流同向性判斷的輸電通道搜索方法。與現(xiàn)有研究成果相比,利用割點成塊可以降低網(wǎng)絡(luò)模型維度,縮短搜索時間,前k路徑法可以消除最短路徑法[3-4]產(chǎn)生的斷面遺漏問題,也可以防止出現(xiàn)系統(tǒng)狀態(tài)圖分區(qū)法[6-7]產(chǎn)生的分區(qū)繁多現(xiàn)象,同時設(shè)計的k值選取公式相比定值原則[5]進一步提升了算法適應(yīng)性。根據(jù)電氣距離選定聯(lián)絡(luò)線并在此基礎(chǔ)上開展節(jié)點收縮分區(qū)符合實際電網(wǎng)分層分區(qū)運行經(jīng)驗[11],同時可以消除先分區(qū)再確定聯(lián)絡(luò)線出現(xiàn)的關(guān)鍵支路收縮問題[9]

1 潮流轉(zhuǎn)移斷面

1.1 潮流轉(zhuǎn)移模型分析

    過載支路的跳閘勢必影響全網(wǎng)的潮流變化,但不同支路的受影響程度是大不相同的,距切除支路電氣距離較近的區(qū)域潮流變化較為明顯,其他區(qū)域變化較小或者無變化[6]。根據(jù)這一特點能夠分析過載支路切除之后,網(wǎng)絡(luò)潮流轉(zhuǎn)移變化的影響范圍與程度。

    假設(shè)電力網(wǎng)絡(luò)任一支路li-j切除后,各發(fā)電機出力和負荷均不發(fā)生變化,則可以認為支路切除前后從節(jié)點i流向節(jié)點j的有功功率并未改變,僅是傳輸路徑發(fā)生了變化。根據(jù)疊加原理可以得到如下關(guān)系式:

    li-j切除后各支路潮流=li-j切除前的潮流+li-j切除后的轉(zhuǎn)移潮流。

    li-j切除后的轉(zhuǎn)移潮流網(wǎng)絡(luò)是僅保留唯一電流源的網(wǎng)絡(luò)模型,且該電流源與切除支路電流大小相同方向相反。

1.2 潮流轉(zhuǎn)移因子與輸電斷面

    如果忽略電力網(wǎng)絡(luò)中各支路電阻且不計電力電子等非線性元件。根據(jù)電路知識可知:各支路轉(zhuǎn)移電流大小與電流源成一固定比例,即有關(guān)系式:

    zndw3-gs1.gif

式中,Iλ.m-t表示支路lm-t的轉(zhuǎn)移電流大小;Is表示電流源大小,即切除支路li-j上的原始電流大小;比例系數(shù)λ為潮流轉(zhuǎn)移因子。

    潮流轉(zhuǎn)移因子用以衡量某一支路斷開對于周圍支路的潮流影響程度,且該值僅和電力網(wǎng)絡(luò)結(jié)構(gòu)與參數(shù)相關(guān)。在此可根據(jù)該參數(shù)定義輸電斷面:如果某一支路斷開,網(wǎng)絡(luò)中潮流轉(zhuǎn)移因子大于閾值λo(一般取0.2-0.3)的全部支路集合稱為斷開支路的輸電斷面,這種輸電斷面一般也稱作潮流轉(zhuǎn)移斷面。

1.3 潮流轉(zhuǎn)移斷面的快速搜索

    根據(jù)電路知識,轉(zhuǎn)移電流僅會流向與斷開支路構(gòu)成回路的各條支路,結(jié)合圖論中割點與塊的性質(zhì)[12]可知轉(zhuǎn)移電流僅分布在斷開支路所在的塊中;根據(jù)歐姆定理,回路阻抗越小,流經(jīng)的轉(zhuǎn)移電流越大,即潮流轉(zhuǎn)移因子越大。因此可以根據(jù)割點性質(zhì)對全網(wǎng)分塊,在斷開支路所在的塊中搜索該支路兩端的前k次最短路徑獲得潮流轉(zhuǎn)移斷面。

    k值可依據(jù)式(2)確定。式中,Pi表示第i次最短路徑,Pi.len為該路徑長度,由路徑組成支路的電抗相加得到,M為路徑長度之比的閾值,取整數(shù)3。

    zndw3-gs2.gif

    根據(jù)該式可知:當Pk+1與P1沒有交集時,路徑Pk+1上的轉(zhuǎn)移電流最大,趨近于電流源的1/4(潮流轉(zhuǎn)移因子趨近于0.25)。k值的選取公式可以確保搜索到所有λ>λo的相關(guān)支路并盡量減少額外支路,但為了防止時間過長一般限制k最大取4。

    前k次最短路徑搜索是圖論中的著名問題,本文依據(jù)文獻[13]提出了利用割點和鄰近節(jié)點的前k次最短路徑搜索方法。

1.4 潮流轉(zhuǎn)移斷面的特點

    根據(jù)潮流轉(zhuǎn)移因子定義的潮流轉(zhuǎn)移斷面涉及到以下關(guān)鍵技術(shù)問題或特點:

    (1)過載切除支路的選取:理論上全網(wǎng)任一支路都存在過載切除的可能,但對它們進行逐一選取與分析是不現(xiàn)實的。實際中的處理方法是在系統(tǒng)運行中發(fā)生線路過載時,才快速啟動潮流轉(zhuǎn)移斷面搜索,為該支路切除后進行緊急控制作出預防指導。

    (2)潮流轉(zhuǎn)移斷面易出現(xiàn)局部支路集聚現(xiàn)象:由于此類輸電斷面的劃分方法存在原理性缺陷,其組成支路通常相互毗鄰,進而出現(xiàn)了斷面微觀現(xiàn)象,無法表征潮流的整體流動趨勢。

    (3)網(wǎng)絡(luò)規(guī)模的適應(yīng)性問題:耗時長短是評價潮流轉(zhuǎn)移斷面搜索方法優(yōu)劣的最重要指標。但無論何種搜索策略其本質(zhì)都是圖論問題,當網(wǎng)絡(luò)規(guī)模較大時搜索時長會成指數(shù)增長。雖然很多學者做出了各種嘗試但效果并不令人滿意,實際使用中往往通過設(shè)定搜索節(jié)點上限來縮短搜索時長[3-5]

2 輸電通道

2.1 輸電通道的識別思路

    實際輸電網(wǎng)架一般帶有明顯的地域特征,表現(xiàn)為區(qū)域間的電網(wǎng)聯(lián)系較區(qū)域內(nèi)更為薄弱,容易遭受破壞。因此各區(qū)域間的聯(lián)絡(luò)線集合成為電網(wǎng)運行控制和監(jiān)視的重點。這種聯(lián)絡(luò)線集合稱作輸電通道,一般帶有電壓等級高、輸電容量大、潮流流向相同的特點。

    聯(lián)絡(luò)線兩端節(jié)點間的電氣距離相較于區(qū)域內(nèi)各鄰近節(jié)點間要大,據(jù)此可以通過比較各鄰近節(jié)點間電氣距離大小來判斷選定聯(lián)絡(luò)線。再依據(jù)鄰近關(guān)系或電氣距離將非聯(lián)絡(luò)線節(jié)點向聯(lián)絡(luò)線節(jié)點收縮分區(qū),經(jīng)過以上操作,整個電網(wǎng)已經(jīng)劃分為若干區(qū)域。最終根據(jù)聯(lián)絡(luò)線搜索區(qū)域間聯(lián)接割集并判斷潮流同向性得到輸電通道。

2.2 確定區(qū)域聯(lián)絡(luò)線

    電氣距離用以度量兩個節(jié)點間的耦合程度,電氣距離越大,表明節(jié)點聯(lián)系越薄弱。根據(jù)文獻[14],定義節(jié)點i和j間的電氣距離為:

    zndw3-gs3.gif

式中,aij表示系統(tǒng)在某一擾動作用下,節(jié)點i和j的電壓變化量比值,即:

    zndw3-gs4.gif

    ai-j可以由全雅克比矩陣計算求得,具體過程可參見文獻[15]。

    定義聯(lián)絡(luò)線為系統(tǒng)中聯(lián)系薄弱的支路,就可以通過比較Di-j來選定聯(lián)絡(luò)線。詳細流程如下:

    (1)依次計算網(wǎng)絡(luò)中任意相鄰節(jié)點間的Di-j,得到列表D;

    (2)比較得出列表D中最大值Dmax,并計算剩余Di-j的平均值Dmean

    (3)若最大距離與平均距離的差值小于給定的閾值,即Dmax-Dmean<Δε,則結(jié)束,否則轉(zhuǎn)(4);

    (4)確定Dmax對應(yīng)的兩相鄰節(jié)點間支路為一條聯(lián)絡(luò)線,在列表D中剔除Dmax后轉(zhuǎn)(2)。

    其中,最大電氣距離與平均電氣距離差值的閾值Δε應(yīng)視具體網(wǎng)絡(luò)給定。

2.3 聯(lián)絡(luò)線節(jié)點收縮分區(qū)與輸電通道搜索

    根據(jù)聯(lián)絡(luò)線選定結(jié)果進行節(jié)點收縮分區(qū)。以圖2為例,具體流程說明如下:

    (1)通過第2.2節(jié)的操作步驟得到聯(lián)絡(luò)線節(jié)點集合SL以及普通節(jié)點集合SP

    (2)對于任一節(jié)點vi∈SP,按如下規(guī)則予以收縮(圖中以虛線支路和空心圓標識該過程):

    ①若vi屬于邊界節(jié)點,則向直接相連的環(huán)上收縮,如v10和v73

    ②若vj∈SL,vi與vj相鄰,且vj唯一,則vi向vj收縮,如v15和v25

    ③若vj∈SL,vi與vj相鄰,但vj不唯一,則vi向電氣距離較近的vj收縮,如v54和v93

    (3)經(jīng)過以上步驟已經(jīng)完成了網(wǎng)絡(luò)初步分區(qū),對于SP中仍未收縮的節(jié)點vk,計算vk到各個初步分區(qū)的平均電氣距離,向較近區(qū)域進行收縮(圖中以陰影區(qū)域標識該過程)。

    經(jīng)過以上步驟完成了基于聯(lián)絡(luò)線的電網(wǎng)節(jié)點收縮分區(qū),之后可結(jié)合圖論算法依據(jù)聯(lián)絡(luò)線搜索各區(qū)域間的網(wǎng)絡(luò)割集[12],通過判斷潮流同向性得到輸電通道。

2.4 電磁環(huán)網(wǎng)

    電磁環(huán)網(wǎng)是指不同電壓等級的輸電線路通過變壓器電磁回路聯(lián)接而構(gòu)成的環(huán)路[16]。電磁環(huán)網(wǎng)是輸電通道的一個特例,可以采用第2.3節(jié)的方法進行搜索,但應(yīng)注意不同電壓等級間的網(wǎng)絡(luò)節(jié)點不能收縮合并。同時應(yīng)綜合考慮潮流轉(zhuǎn)移因子的概念,以滿足高電壓線路斷開對低電壓回路潮流有較強的影響能力。

2.5 輸電通道的特點

    基于聯(lián)絡(luò)線節(jié)點收縮分區(qū)形成的輸電通道涉及到以下關(guān)鍵技術(shù)或特點:

    (1)聯(lián)絡(luò)線的選取:利用電氣距離概念判定相鄰節(jié)點間的聯(lián)系薄弱程度,通過設(shè)定電氣距離閾值選定聯(lián)絡(luò)線并開展節(jié)點收縮分區(qū)。

    (2)與潮流轉(zhuǎn)移斷面不同,輸電通道可以在宏觀上表征潮流的整體流動趨勢,與某一具體支路無關(guān)。

    (3)電磁環(huán)網(wǎng)屬于輸電通道的一個特例。

3 仿真分析

3.1 潮流轉(zhuǎn)移斷面

    新英格蘭39節(jié)點系統(tǒng)存在割點v16和v26,把全網(wǎng)分成了三塊,考慮支路l5-6切除的情況開展潮流轉(zhuǎn)移斷面搜索。根據(jù)第1.3節(jié)內(nèi)容,只需在l5-6所在塊中進行搜索,如圖1所示,搜索結(jié)果列于表1。

zndw3-t1.gif

zndw3-b1.gif

    根據(jù)第1.3節(jié)式(2)關(guān)于k值的選取公式,由于P3.len/P1.len=4.9>3,得到k=2。

    若參照文獻[5]的方法會包含這兩條最短路徑以外的大量額外支路;參照[7]的方法會漏選3條支路且需要分成17個較多區(qū)域。

    表1中的潮流轉(zhuǎn)移因子是根據(jù)文獻[17]的理論分析計算而來。由該表可以發(fā)現(xiàn)搜索到的2條路徑P1和P2所有組成支路的潮流轉(zhuǎn)移因子絕對值都大于λo(0.2~0.3),計算還發(fā)現(xiàn)除去這些支路以外,最大的潮流轉(zhuǎn)移因子絕對值僅為0.054,出現(xiàn)在路徑P3的支路l4-3上。搜索結(jié)果表明支路l5-6過載切除引起的潮流轉(zhuǎn)移因子較大的支路全部包含在潮流轉(zhuǎn)移斷面S中且沒有產(chǎn)生額外支路。

3.2 輸電通道

    輸電通道的搜索方法通過IEEE118節(jié)點系統(tǒng)[18]進行驗證。設(shè)置聯(lián)絡(luò)線選取時電氣距離的差值閾值Δε=1.5,搜索共計發(fā)現(xiàn)13條聯(lián)絡(luò)線,占線路總數(shù)的7%,在圖2中以加粗線條表示。搜索結(jié)果列于表2。

zndw3-t2.gif

zndw3-b2.gif

    由圖2和表2可知IEEE118節(jié)點共搜索到輸電通道6個,每個輸電通道都有相應(yīng)的聯(lián)絡(luò)線,各通道內(nèi)潮流流向一致。搜索結(jié)果表明基與電氣距離選定聯(lián)絡(luò)線并在此基礎(chǔ)上進行節(jié)點收縮分區(qū)開展輸電通道的搜索方法正確有效。

    本文采用的電氣距離概念能夠表征各節(jié)點間電氣聯(lián)系緊密程度,以此選定的聯(lián)絡(luò)線具有實際運行意義。設(shè)計的搜索方法可以消除先分區(qū)再選聯(lián)絡(luò)線產(chǎn)生的關(guān)鍵支路收縮問題[9]

4 結(jié)論

    本文從潮流轉(zhuǎn)移斷面、輸電通道以及電磁環(huán)網(wǎng)三個不同角度進行輸電斷面的全方位認識分析工作,并在此基礎(chǔ)上進行搜索方法研究與實例仿真。

    潮流轉(zhuǎn)移斷面源于支路連續(xù)過載跳閘現(xiàn)象,能夠描述某一支路過載切除時潮流變化的影響范圍和程度,可以根據(jù)割點和鄰近節(jié)點的前k次最短路徑搜索確定;輸電通道根據(jù)聯(lián)絡(luò)線和網(wǎng)絡(luò)結(jié)構(gòu)劃分確定,能夠描述不同區(qū)域間的潮流輸送特點,可以根據(jù)電氣距離確定聯(lián)絡(luò)線再進行節(jié)點收縮分區(qū),根據(jù)圖論搜索區(qū)域間割集以確定;電磁環(huán)網(wǎng)作為輸電通道的一個特例,能夠描述多電壓等級間的耦合特點。

參考文獻

[1] 尚亞男.輸電斷面認識與搜索方法研究[D].武漢:華中科技大學,2013.

[2] U.S.-Canada Power System Outage Task Force. Interim report: causes of the august 14th blackout in the United States and Canada[R].2003,11.

[3] 倪宏坤,徐玉琴.基于動態(tài)規(guī)劃原理分支界限算法的關(guān)鍵輸電斷面搜索方法[J].華北電力大學學報,2009,36(4):11-15.

[4] 程臨燕,張保會,郝治國,等.基于最小基本回路集合的潮流轉(zhuǎn)移快速搜索[J].電力系統(tǒng)自動化,2010,34(18):21-26.

[5] 王增平,李剛,任建文.基于前K最短路徑的輸電斷面搜索新算法[J].電工技術(shù)學報,2012,27(4):193-201.

[6] 張保會,姚峰,周德才,等.輸電斷面安全性保護及其關(guān)鍵技術(shù)研究[J].中國電機工程學報,2006,26(21):1-7.

[7] 周德才,張保會,姚峰,等.基于圖論的輸電斷面快速搜索[J].中國電機工程學報,2006,26(12):34-41.

[8] 趙峰,孫宏斌,張伯明.基于電氣分區(qū)的輸電斷面及其自動發(fā)現(xiàn)[J].電力系統(tǒng)自動化,2010,24(5):1-6.

[9] 楊燕,文福拴,王潔,等.計及關(guān)鍵支路的輸電斷面自動確定方法[J].電力系統(tǒng)自動化,2012,36(14):6-10.

[10] 苗世洪,馬帥,尚亞男,等.基于割點和路徑搜索的輸電斷面快速識別方法[J].電力系統(tǒng)自動化,2014,38(2):39-45.

[11] 劉廣利,李響.基于輸電斷面劃分原則及方法[J].東北電力技術(shù),2005(11):13-16.

[12] 龔劬.圖論與網(wǎng)絡(luò)最優(yōu)化算法[M].重慶:重慶大學出版社,2009.

[13] 白軼多,胡鵬,夏蘭芳,等.關(guān)于k次短路徑問題的分析與求解[J].武漢大學學報,2009,39(4):492-494.

[14] Lagonotte P,Sabonnadiè J C,Léost J Y,et al.Structural analysis of the electrical system: application to secondary voltage control in France[J].IEEE Trans on Power Systems,1989,4(2):479-486.

[15] 盧志剛,秦四娟,常磊,等.基于Dempster-Shafer證據(jù)理論的無功/電壓控制分區(qū)[J].電網(wǎng)技術(shù),2010,34(10):99-104.

[16] 秦吉,黃方能.電磁環(huán)網(wǎng)運行風險分析[J].電力安全技術(shù),2011,13(12):26-28.

[17] 吳際舜.電力系統(tǒng)靜態(tài)安全分析[M].上海:上海交通大學出版社,1985.

[18] Power Systems Test Case Archive[EB/OL].http://www.ee.washington.edu/research/pstca/.

此內(nèi)容為AET網(wǎng)站原創(chuàng),未經(jīng)授權(quán)禁止轉(zhuǎn)載。
欧美激情办公室aⅴ_国产欧美综合一区二区三区_欧美午夜精品久久久久免费视_福利视频欧美一区二区三区

          久久久之久亚州精品露出| 久久九九热re6这里有精品| 一区二区三区免费看| 亚洲视频一区| 国产欧美91| 亚洲人成在线观看网站高清| 欧美国产日本在线| 在线综合视频| 国产精品网站在线| 久久人人97超碰人人澡爱香蕉| 国产目拍亚洲精品99久久精品| 久久久国产精品亚洲一区| 亚洲国产日韩美| 国产精品hd| 久久综合久久88| 中文欧美字幕免费| 99精品视频一区| 欧美午夜大胆人体| 久久成人资源| 99精品国产在热久久下载| 国产精品视频一区二区高潮| 久久一二三国产| 亚洲性感美女99在线| 黄色在线一区| 欧美天天视频| 噜噜噜噜噜久久久久久91 | 一本色道**综合亚洲精品蜜桃冫 | 久久久久久999| 一本久道久久综合中文字幕| 国产精品日韩专区| 国产精品久久久久久久久免费樱桃| 亚洲午夜一区二区| 国产伦精品一区二区三区| 免费视频亚洲| 久久爱91午夜羞羞| 亚洲天堂网在线观看| 最新日韩在线| 国内精品久久久久久| 久久久一区二区| 亚洲视频在线观看三级| 亚洲一区二区三区影院| 最新成人av网站| 国产午夜精品久久久| 欧美日韩国产综合视频在线观看| 亚洲视频综合在线| 精品不卡视频| 国产精品中文在线| 狠狠色综合色区| 欧美成人自拍视频| 香蕉国产精品偷在线观看不卡| 欧美日韩综合视频| 免费成人激情视频| 欧美一级在线视频| 国产农村妇女毛片精品久久莱园子| 久久久久国产精品厨房| 另类天堂av| 一区二区三区欧美在线观看| 亚洲视频在线观看视频| 亚洲人成网站精品片在线观看| 国产一区日韩欧美| 国产精品高潮在线| 国产日韩欧美自拍| 国产精品成人av性教育| 国产精品日韩电影| 久久精品亚洲一区二区三区浴池 | 国产一区二区三区四区hd| 蜜臀91精品一区二区三区| 欧美久久久久久久久| 欧美激情亚洲激情| 免费国产一区二区| 久久精品国产在热久久| 亚洲一区三区视频在线观看| 久久精品99| 欧美三级网址| 国产精品女人毛片| 在线观看亚洲精品视频| 欧美日韩精品二区第二页| 久久伊人免费视频| 国产精品国产精品| 久久婷婷蜜乳一本欲蜜臀| 欧美高清视频一区二区| 欧美精品97| 国产美女扒开尿口久久久| 国产精品黄视频| 亚洲国产精品va在线看黑人| 影音先锋日韩资源| 影音先锋久久久| 亚洲欧美日韩在线观看a三区| 亚洲国产二区| 亚洲欧美日韩一区二区三区在线 | 国产色婷婷国产综合在线理论片a| 国产伦精品一区二区| 亚洲片在线观看| 久久九九99视频| 国产精品久久久久久久久免费| 久久久久免费观看| 欧美视频在线观看 亚洲欧| 在线观看中文字幕亚洲| 国产精一区二区三区| 国产精品一区二区三区乱码| 国产午夜精品全部视频播放 | 欧美日韩亚洲激情| 美女精品国产| 国产欧美日韩精品一区| 一级成人国产| 欧美黄色成人网| 在线看无码的免费网站| 久久久久久穴| 亚洲永久免费| 久久天天综合| 国产亚洲一区二区三区在线观看 | 国产精品日韩欧美一区二区| av成人免费在线观看| 久久精品亚洲| 国产精品久久久久久久久免费 | 1024亚洲| 久久久久成人网| 国产亚洲一级高清| 狠狠色综合播放一区二区| 性欧美暴力猛交69hd| 亚洲在线国产日韩欧美| 欧美伦理91i| 国产精品素人视频| 在线观看欧美一区| aⅴ色国产欧美| 欧美激情综合五月色丁香小说| 欧美成人dvd在线视频| 欧美日韩久久精品| 亚洲激情av在线| 欧美成人性网| 91久久极品少妇xxxxⅹ软件| 亚洲国产精品高清久久久| 亚洲电影下载| 久久这里有精品视频| 激情综合在线| 亚洲高清不卡在线观看| 亚洲一区在线看| 欧美国产日本高清在线| 亚洲开发第一视频在线播放| 久久精品伊人| 激情亚洲网站| 欧美成人自拍视频| 国产综合久久久久久鬼色| 亚洲免费av网站| 久久午夜激情| 亚洲二区在线| 欧美一区二区国产| 欧美日韩一区二区在线观看视频 | 国产精品区一区| 午夜一区二区三区不卡视频| 欧美日韩国产小视频| 国产一区二区按摩在线观看| 久久精品视频在线看| 一区二区在线观看视频| 欧美成人中文| 亚洲视频在线观看网站| 国产毛片久久| 久久婷婷久久| 国产一区二区三区久久 | 中文高清一区| 国产日产精品一区二区三区四区的观看方式 | 久久久久国产精品厨房| 亚洲福利视频一区二区| 久久久久久久综合狠狠综合| 亚洲国产精品一区在线观看不卡| 久久一区中文字幕| 99伊人成综合| 国产亚洲视频在线| 欧美a一区二区| 亚洲欧美国内爽妇网| 欧美日韩精品免费看| 午夜日韩在线观看| 国产精品国产成人国产三级| 久久激情五月丁香伊人| 亚洲人成人99网站| 久久国产精彩视频| 国产欧美日韩精品丝袜高跟鞋| 亚洲丝袜av一区| 精品av久久707| 久久性天堂网| 一本色道久久88精品综合| 国产亚洲毛片| 欧美日韩综合另类| 久久免费视频观看| 一区二区欧美亚洲| 欧美色网一区二区| 久久久久se| 亚洲字幕一区二区| 国产精品qvod| 在线亚洲一区| 狠狠色丁香婷婷综合久久片| 99精品视频免费观看视频| 欧美老女人xx| 99re在线精品| 欧美午夜精品久久久久久人妖 | 欧美日韩一区二区三区| 久久久久91| 亚洲色图综合久久| 91久久精品国产91久久性色tv| 欧美一区不卡|