亚洲免费av电影一区二区三区,日韩爱爱视频,51精品视频一区二区三区,91视频爱爱,日韩欧美在线播放视频,中文字幕少妇AV,亚洲电影中文字幕,久久久久亚洲av成人网址,久久综合视频网站,国产在线不卡免费播放

        ?

        基于阻抗模型故障特征匹配技術(shù)的配電網(wǎng)單相接地故障測距研究

        2015-12-30 09:19:45康忠健劉睿瑩
        電網(wǎng)與清潔能源 2015年8期
        關(guān)鍵詞:測量點零序測距

        康忠健,劉?,?/p>

        (中國石油大學(xué)(華東) 信息與控制工程學(xué)院,山東 青島 266580)

        基于阻抗模型故障特征匹配技術(shù)的配電網(wǎng)單相接地故障測距研究

        康忠健,劉?,?/p>

        (中國石油大學(xué)(華東) 信息與控制工程學(xué)院,山東 青島 266580)

        單相接地故障是配電網(wǎng)最常發(fā)生的故障,準(zhǔn)確而快速地確定配電網(wǎng)單相接地故障位置對排除系統(tǒng)故障,減少事故損失具有非常重要的作用。針對當(dāng)前配電網(wǎng)單相接地故障精確定位問題,通過采用網(wǎng)絡(luò)故障時測量點處的電壓預(yù)存信息,同時利用故障中饋線出口處測得的零序電壓和零序電流信號,提出了一種采用差分進化算法的基于阻抗模型故障特征匹配的配電網(wǎng)單相接地故障測距的新方法。該方法基于系統(tǒng)阻抗模型定義了系統(tǒng)故障特征適應(yīng)度指標(biāo),利用枚舉法確定故障分支線路,并采用差分進化算法精確定位故障位置。中石化西北油田艾丁變10 kV線路配電網(wǎng)仿真結(jié)果表明了所提出的配電網(wǎng)單相接地故障測距方法的有效性。

        配電網(wǎng);單相接地;故障測距;差分進化算法;故障特征匹配

        在我國,中壓配電網(wǎng)的中性點大部分采用的是非直接接地方式,包括中性點不接地方式[1]和中性點經(jīng)消弧線圈接地方式[2],屬于小電流接地配電系統(tǒng)[3]。在中性點非直接接地的電力系統(tǒng)中,發(fā)生概率最高的故障是單相接地故障。而配電系統(tǒng)采用小電流接地運行方式的好處是當(dāng)配電系統(tǒng)出現(xiàn)單相接地故障時,故障線路流過的電流僅為流過線路對地電容的電流[4],數(shù)值不大,可在電力系統(tǒng)安全規(guī)定[5]內(nèi)繼續(xù)維持供電運行一至兩個小時,這樣可以提高電力系統(tǒng)的供電可靠性。但是,此時非故障相的相電壓提高為線電壓,對電力設(shè)備的絕緣造成損壞,影響電力系統(tǒng)的安全運行。

        現(xiàn)有的配電網(wǎng)故障測距方法可以分為阻抗分析法[6]、行波 法[7]、注入法[8]、故障區(qū)間定位方法和故障特征匹配測距[9]、智能法[10]等。各種方法均有其一定的優(yōu)缺點,均不能在配電網(wǎng)中完全適用,大都無法在只有單個測量點的單電源配電網(wǎng)中應(yīng)用。

        為解決單電源配電網(wǎng)單相接地單測量點故障測距問題,本文提出了一種新的定位方法。該方法在將各節(jié)點故障時測量點電壓信息預(yù)存完畢的前提下,建立了配電網(wǎng)零序阻抗模型,對配網(wǎng)的零序故障特征進行了提取,利用枚舉法和差分進化算法進行區(qū)間定位和故障精確定位,是一種基于配網(wǎng)區(qū)饋線出口處電壓和電流信息的配電網(wǎng)單相接地故障測距方法,并應(yīng)用實際配電網(wǎng)模型對該方法進行了全面的仿真驗證。

        1 配電網(wǎng)阻抗模型的建立

        本文提出的測距原理全部基于零序網(wǎng)絡(luò)進行分析,故采用單相零序阻抗模型進行建模。

        中性點經(jīng)消弧線圈接地的配電網(wǎng)阻抗模型包括主電源、變壓器、線路以及負(fù)荷的阻抗模型。各個模型中所需各項參數(shù)均為已知量。配電網(wǎng)的各部分零序阻抗模型建立過程如下。

        1.1 主電源及所帶消弧線圈的零序阻抗模型建立

        由戴維南定理可知,系統(tǒng)電源阻抗模型可等值為一個含內(nèi)阻的理想電壓源。采用文獻(xiàn)[11]中電源等效模型的建模方法。其中,配電系統(tǒng)在主饋線處(電源出線位置)的零序電流和電壓皆可通過測量獲得。從電源點看進去,電源內(nèi)阻和消弧線圈可看作一個整體對待。 主電源及其所帶消弧線圈整體的零序阻抗如式(1)所示。

        式中:Us0為配網(wǎng)饋線處測得的零序電壓,規(guī)定由高電位指向低電位為電壓的正方向;IF0為同一時刻配網(wǎng)饋線處測得的零序電流,規(guī)定從母線流向線路為電流的正方向。

        1.2 配電線路的阻抗模型建立

        本文的配電線模型采用文獻(xiàn)[3]提出的π型等效模型,如圖1所示。

        根據(jù)參考文獻(xiàn)[3],線路的阻抗參數(shù)(r0、x0、b0)可由已知的線路型號查表求得,線路長度l也為已知量。

        圖1 π型等效模型Fig.1 The PI equivalent model of line

        1.3 配電網(wǎng)零序節(jié)點導(dǎo)納矩陣的形成

        由于引入變壓器和負(fù)載的零序模型會增加導(dǎo)納矩陣的形成難度,并且變壓器和負(fù)載的引入與否并不影響該方法的正確性。故模型中假設(shè)負(fù)荷變壓器采用Y/Δ接線方式,負(fù)載均采用Δ連接,則在建立系統(tǒng)的零序等效阻抗時不考慮負(fù)荷和變壓器的零序阻抗。

        根據(jù)配電網(wǎng)的網(wǎng)絡(luò)結(jié)構(gòu)和各部分模型參數(shù),形成n個節(jié)點的單電源供電系統(tǒng)的零序節(jié)點導(dǎo)納矩陣[1]如式(4)所示。

        式中:Y0nn代表n節(jié)點的零序自導(dǎo)納,Y0nm則代表n節(jié)點和m節(jié)點之間的互導(dǎo)納。

        至此,配電系統(tǒng)的零序阻抗模型建立完畢。

        2 基于阻抗模型故障特征匹配的配電網(wǎng)單相接地故障測距方法

        2.1 單電源配電網(wǎng)的故障特征提取

        同理,可以得到配電網(wǎng)系統(tǒng)的零序節(jié)點阻抗方程

        系統(tǒng)的節(jié)點阻抗矩陣為

        式中:Ui,Ii分別是節(jié)點i的零序電壓相量和零序注入電流相量。Y0n×n與Z0n×n互為逆矩陣,即

        由零序節(jié)點電壓方程可知,當(dāng)節(jié)點處注入的短路零序電流為I˙fj0時在各節(jié)點產(chǎn)生的故障零序電壓分量ΔV為

        由式(8)可知,當(dāng)節(jié)點j處注入的短路零序電流為I˙fj0時,在測量點處產(chǎn)生的零序故障電壓分量為

        由于電源測量點κ處的故障電壓分量可由該節(jié)點κ所測量的故障中和故障前零序電壓測量值相減得到,即

        由式(9)和式(10)可得由電源測量點處的故障電壓分量計算故障點j節(jié)點的零序故障電流為

        圖2為本文提出方法分析所采用的簡圖,該配網(wǎng)為單電源配網(wǎng),圖中節(jié)點1是電源節(jié)點(即本文測距所用的測量點),節(jié)點2是與電源節(jié)點直接相連的出線節(jié)點。因本文指對配網(wǎng)中的單相接地故障進行討論,而故障前零序電壓電 流可近 似為 0,故無需對故障前的狀態(tài)進行測量。

        各節(jié)點故障時測量點處的零序電壓可在配網(wǎng)模型中仿真得到并預(yù)存到數(shù)據(jù)庫中。其中當(dāng)故障點發(fā)生在 j處,測量點(即圖2中的節(jié)點1)處的預(yù)存零序電壓信息定義為U0(j)。測量點故障中的零序電壓和零序電流,分別記為U1H、I1H。

        圖2 分析原理簡圖Fig.2 Principle diagram

        通過預(yù)存信息可求得當(dāng)故障節(jié)點為 j時,故障處的對地零序電流:

        可求出與電源節(jié)點相連的節(jié)點(即圖2中的節(jié)點2)的零序電壓變化量:

        由1、2節(jié)點的電壓可求出故障后節(jié)點1、2之間的零序電流:

        定義配電網(wǎng)中某一節(jié)點 j發(fā)生故障時計算出的電源出線上的零序電流與實際測得的電源出線上的零序電流的絕對偏差為本方法所設(shè)定的配電網(wǎng)的故障特征值。即在j節(jié)點的故障特征值為:

        式中:I0(j)為通過計算得出的 j點為故障點時配網(wǎng)主饋線首端流通的零序電流;I1H為實際測量獲得的配網(wǎng)主饋線首端流通的零序電流。

        由電網(wǎng)絡(luò)理論可知:當(dāng) j越靠近配電網(wǎng)故障點時系統(tǒng)的故障特征值越小,當(dāng) j為配電網(wǎng)故障點時系統(tǒng)的故障特征值為0。通過比較判斷E的大小,就可以判斷故障區(qū)間和故障位置。

        2.2 配電網(wǎng)單相接地故障測距方法

        本文所提出的配電網(wǎng)單相故障測距方法分為兩步。首先,采用枚舉法識別出故障的潛在區(qū)域;其次采用差分進化算法實現(xiàn)故障的精確定位。

        2.2.1 故障潛在區(qū)域的識別

        假設(shè)故障發(fā)生在節(jié)點j計算此時系統(tǒng)故障特征量E(j),通過枚舉法假設(shè)故障依次發(fā)生在配網(wǎng)各個節(jié)點,根據(jù)局故障特征值的定義各個計算故障特征值的大小,形成故障特征量矩陣E。由上文分析可知,故障點處的故障特征值最小,又為了消除計算誤差等因素。故選擇故障特征量E最小的3個節(jié)點為故障區(qū)間。因此,與故障能量最小的3個節(jié)點相關(guān)聯(lián)的線路即視為可能的故障線路,從而確定出故障潛在區(qū)域,再逐一進行判斷,得出最可能的2條線路。具體實現(xiàn)方法如下[12-14]:

        假設(shè)單相接地故障點依次發(fā)生在可能故障線路的中點處。從網(wǎng)絡(luò)結(jié)構(gòu)上看,該中點相當(dāng)于在配電網(wǎng)中新增了一個節(jié)點,記為第n+1節(jié)點,根據(jù)此時的系統(tǒng)結(jié)構(gòu)對系統(tǒng)線路的節(jié)點編號和線路長度等參數(shù)進行修正。如圖3所示,假設(shè)判斷出的關(guān)聯(lián)線路為第i條線(首末節(jié)點分別為節(jié)點2和節(jié)點j),則新增節(jié)點即為第i條線的中點。重新修正時將第i條線的末節(jié)點改為n+1節(jié)點,長度改為原有的一半。增設(shè)第n+1條線首節(jié)點即為節(jié)點n+1,末節(jié)點為節(jié)點j,長度為第i條線原長的一半。根據(jù)修正后的結(jié)構(gòu)修正零序?qū)Ъ{矩陣,并結(jié)合前文的分析可知,若故障發(fā)生在第i條線中點附近,則在增設(shè)的n+1節(jié)點處配電網(wǎng)的故障特征值理論上應(yīng)為最小。

        圖3 增加節(jié)點后的原理分析圖Fig.3 Principle diagram after adding a new node

        通過這種方法得到各故障特征值,其中數(shù)值最小的2個值所在的線路即為最可能發(fā)生故障的線路。

        故障區(qū)間識別的流程圖如圖4所示。

        2.2.2 故障精確定位

        在從故障區(qū)域中判斷出2條最可能故障線路后,仍需進一步對故障進行精確定位。假設(shè)單相接地故障點依次發(fā)生在判斷所得的2條最可能的故障線路上,且故障發(fā)生點距該線路首節(jié)點x處,同理,等效為在故障點處新增一個節(jié)點n+1節(jié)點,按照上節(jié)介紹的方法對系統(tǒng)線路參數(shù)進行修正。一個節(jié)點的配電網(wǎng)的零序阻抗矩陣模型維數(shù)將變?yōu)椋╪+1)×(n+1)。根據(jù)分析可知,在實際故障發(fā)生點所處的n+1節(jié)點發(fā)生故障時,配電網(wǎng)的故障特征值理論上應(yīng)為最小。

        利用差分進化算法[13]對故障條件進行自動搜索,差分進化算法的基礎(chǔ)是選擇、交叉和變異操作,而構(gòu)成差分進化算法的關(guān)鍵則是個體適應(yīng)度函數(shù)和變量的選取。

        圖4 故障區(qū)間定位流程圖Fig.4 Flow chart of the fault area location

        1)適應(yīng)度函數(shù)。

        適應(yīng)度函數(shù)(fitness function)是評價個體優(yōu)劣的主要判據(jù),也是個體優(yōu)勝劣汰的準(zhǔn)則。差分進化算法的3種差分操作也是以適應(yīng)度函數(shù)為基礎(chǔ)的。根據(jù)研究問題的不同,適應(yīng)度函數(shù)的具體定義也隨之變化,本課題的適應(yīng)度函數(shù)是以故障點n+1的故障特征值來定義的,即故障點的故障特征值最小,適應(yīng)度函數(shù)定義如式(16)。

        2)變量。

        本文選取潛在故障線路上的故障距離x為變量。

        當(dāng)搜索到的故障點越接近真正的故障點的時候,適應(yīng)度函數(shù)的值將會越小,當(dāng)適應(yīng)度函數(shù)值滿足一定的精度要求時,輸出測距結(jié)果。因此,精確故障測距結(jié)果就為故障特征值最小時所對應(yīng)的故障距離,具體流程圖如圖5所示。

        3 數(shù)值仿真驗證

        3.1 仿真模型

        在一個實際的中石化西北油田變電站的10 kV線路TA-2S上對本文提出的方法進行測試。系統(tǒng)拓?fù)鋱D如圖6所示,線路信息如表1所示。其中,配網(wǎng)線路全長為58.2 km,節(jié)點1~18為主干線節(jié)點編號,節(jié)點19~64為分支線節(jié)點編號。已知主干線導(dǎo)線規(guī)格為LGJ-120/20 mm2,分支線規(guī)格為LGJ-50/8 mm2。

        圖5 故障精確定位流程圖Fig.5 Flow chart of the precise fault location

        圖6 仿真系統(tǒng)的拓?fù)浣Y(jié)構(gòu)Fig.6 The topology of the simulation model

        測量點配置在配網(wǎng)饋線處,即等效電源的出線點。通過simulink仿真中的測量及數(shù)據(jù)導(dǎo)出模塊可以獲得饋線處零序電壓和流通的零序電流。通過對測得的數(shù)據(jù)進行全波傅里葉處理,可以得到相應(yīng)的有效值,按照圖3的流程圖進行matlab編程,就可求出不同節(jié)點的故障特征值,作為測距依據(jù)。

        3.2 判斷故障區(qū)域仿真結(jié)果

        故障測距的第一步是判斷故障區(qū)域,系統(tǒng)線路數(shù)為63條,將A相單相接地故障設(shè)定在線路的中點位置。在判斷出的故障關(guān)聯(lián)節(jié)點中進行排查,通過對故障特征值大小的比較,從中判別出2條最可能的故障線路。當(dāng)節(jié)點42發(fā)生單相接地短路故障時,求得n個故障特征值如圖7所示,從圖中可知故障特征值最小的節(jié)點在42節(jié)點附近。潛在故障線路判斷部分仿真結(jié)果如表2所示。

        表1 線路信息Tab.1 The line parameters

        圖7 節(jié)點42故障時各節(jié)點故障特征值曲線Fig.7 Fault characteristic value of each node when No.42 node occurs single-phase ground fault

        3.3 故障精確定位結(jié)果

        選取故障線路、故障距離和接地電阻為變量,配電網(wǎng)的系統(tǒng)零序阻抗模型下的故障特征值作為適應(yīng)度函數(shù),利用差分進化算法對故障條件進行自動搜索,精確故障測距的結(jié)果就為故障特征值最小時對應(yīng)的故障距離,精確測距仿真結(jié)果如表3所示。

        表2 潛在故障線路仿真結(jié)果Tab.2 The results of potential fault branches

        表3 精確測距仿真結(jié)果Tab.3 The simulation results of precise fault location

        從表2和表3可得出,故障分支能夠被準(zhǔn)確地判別出來,精確故障測距與實際故障點距離最大誤差為128 m,測距誤差為2個桿距左右。本文所提出的配電網(wǎng)單相接地故障測距方法可以準(zhǔn)確地確定出故障分支和故障點位置,具有很高的測距精度。

        4 結(jié)論

        本文針對配電網(wǎng)單相接地故障測距這一難題,提出了一種通過預(yù)存網(wǎng)絡(luò)故障時測量點處的電壓信息的基于阻抗模型故障特征匹配的單電源配電網(wǎng)單相接地短路單端故障測距的方法。仿真結(jié)果表明:本文所提出的配電網(wǎng)單相接地故障測距方法可以快速準(zhǔn)確地確定出故障分支和故障點位置,具有很高的測距精度。該方法只需饋線出口處故障中零序電壓和零序電流信號,無需故障前的電壓電流信息,現(xiàn)場應(yīng)用方便,具有廣闊的應(yīng)用前景。

        [1]MIAOYouzhong,SUNYaming,YANGHua.Anew principle of transient current grounded relay for feeder in ungrounded distribution systems[J].Proceedings of the CSEE,2004,24(3):62-66.

        [2]陳怡,蔣平,萬秋蘭,等.電力系統(tǒng)分析[M].北京:中國電力出版社,2005.

        [3]康忠健,田愛娜,馮艷艷.含DG配電網(wǎng)故障測距技術(shù)綜述[J].電網(wǎng)與清潔能源,2013,29(4):6-12. KANG Zhongjian,TIAN Aina,F(xiàn)ENG Yanyan.A review of fault 1ocation technology in distribution network with DG[J]. Power System and clean Energy,2013,29(4):6-12(in Chinese).

        [4]康忠健,田愛娜,馮艷艷.基于HHT配電網(wǎng)不對稱接地故障測距新方法[J].電網(wǎng)與清潔能源,2012,28(11):1-6. KANG Zhongjian,TIAN Aina,F(xiàn)ENG Yan-yan.A new method of asymmetric ground fault 1ocation based on the Hilbert-Huang Transform in distribution network[J].Power System and clean Energy,2012,28(11):1-6(in Chinese).

        [5]KRISHNATHEVAR R,NGU E E.Generalized impedance-based fault 1ocation for distribution systems[J].IEEE Transactions on Power Delivery,2012,27(1):449-451.

        [6]CAI Yumei,HE Zhengyou,WANG Zhibing,et al. Application of travelling wave based fault 1ocation in 10 kV railway automatic blocking and continuous power transmission 1ines[J].Power System Technology,2005,29(1):15-19.

        [7]SANG Zaizhong,PAN Zhencun,LI Lei.“S’s signal injection”and fault 1ine identifying and distance measuring and fault point 1ocating technique[J].Proceedings of electric Power System and Automation,1998(4):35-39.

        [8]吳剛.配電網(wǎng)單相接地故障的區(qū)間定位和測距[D].宜昌:三峽大學(xué),2012.

        [9]張鈞.配電網(wǎng)智能故障診斷與諧波源定位研究[D].成都:西南交通大學(xué),2012.

        [10]BRAHMA S M.Faultlocation in power distribution system with penetration of distributed generation[J].IEEE Transactions on Power Delivery,2011,26(3):1545-1553.

        [11]FENG Yanyan,KANG Zhongjian.Research on fault 1ocation in distributed network with DG based on complex correlation thevenin equivalent and strong tracking fil-ter[J].Sensors&Transducers Journal,2014,162(1):152-160.

        [12]陳根,唐焱,王新橋.基于ATP的高壓電纜金屬護套多點接地故障仿真[J].高壓電器,2014,50(4):49-51. CHEN Gen,TANG Yan,WANG Xinqiao.Simulation of multipoint grounding faults in metal sheath of high-voltage cable based on ATP[J].High Voltage Apparatus,2014,50(4):49-51(in Chinese).

        [13]甄昊涵,郭潔,李江濤,等.UHVDC閥閉鎖移相參數(shù)對站內(nèi)接地故障過電壓的影響[J].高壓電器,2014,50(5):42-45. ZHEN Haohan,GUO Jie,LI Jiangtao,et al.Impact of phase-shift parameters on overvoltage caused by ground fault in convertor station of UHVDC [J].High Voltage Apparatus,2014,50(5):42-45(in Chinese).

        [14]YANG Qiwen,CAI Liang,XUE Yuncan.A survey of differential evolution algorithms[J].PR&AI,2008,21(4):506-513.

        (編輯 馮露)

        A Fault Location Method for Single-Phase Grounding Fault in Distribution Network

        KANG Zhongjian,LIU Ruiying
        (College of Information and Control Engineering,China University of Petroleum,Qingdao 266580,Shandong,China)

        The single-phase grounding fault is the most common fault among all types of faults in the distribution network.It is of great importance to 1ocate the single phase grounding fault accurately and quickly reduce economic 1osses. Aiming at the 1ocation of single-phase grounding fault in the distribution network,this paper puts forward a new fault 1ocation method which uses a differential evolution algorithm and the Impedance Model Fault Characteristics Matching Techno-1ogy.The proposed method uses stored fault voltage information,the measured zero sequence voltage and current signals of the feeder in the fault duration time.And the fitness index of the system fault characteristics is defined.After identifying the potential fault branches by using the enumeration method,the fault can be 1ocated accurately by using the differential evolution algorithm.The detailed simulation results in a 10 kV distribution network model of the SINPOPEC West-Northern oil field is given in this paper.The simulation results suggest the validity of the proposed single-phase grounding fault 1ocation method.

        2015-04-11。

        康忠?。?971),男,博士,教授,博士生導(dǎo)師,研究方向為電力系統(tǒng)分析與控制、電力系統(tǒng)故障檢測與診斷和油田自動化設(shè)備開發(fā)及相關(guān)教學(xué)與科研工作;

        劉?,摚?990),女,碩士研究生,研究方向為電力系統(tǒng)故障檢測與診斷;

        曹 揚(1991),男,碩士研究生,研究方向為電力系統(tǒng)故障檢測與診斷。

        國家自然科學(xué)基金項目(61271001),中國石油大學(xué)(華東)自主創(chuàng)新科研計劃項目(14CX05039A)。

        Project Supported by National Natural Science Foundation of China(61271001);the Fundamental Research Funds for the Central Universities of China under Grant(14CX05039A).

        1674-3814(2015)08-0040-07

        TM726

        A

        KEY W0RDS:distribution network;single-phase grounding fault;fault 1ocation;the differential evolution algorithm;fault characteristics matching

        猜你喜歡
        測量點零序測距
        飛機部件數(shù)字化調(diào)姿定位測量點的優(yōu)選與構(gòu)造算法
        淺析沖壓件測量點的規(guī)劃
        模具制造(2019年10期)2020-01-06 09:13:08
        類星體的精準(zhǔn)測距
        科學(xué)(2020年3期)2020-01-06 04:02:51
        基于CAD模型的三坐標(biāo)測量機測量點分布規(guī)劃
        PM2.5空中探測器的設(shè)計
        淺談超聲波測距
        電子制作(2017年7期)2017-06-05 09:36:13
        6kV供電系統(tǒng)零序保護誤動作處理
        大電流接地系統(tǒng)接地故障時零序電壓分析
        基于PSOC超聲測距系統(tǒng)設(shè)計
        相對差分單項測距△DOR
        太空探索(2014年1期)2014-07-10 13:41:50
        亚洲精品综合久久中文字幕| 69精品人人人人| 亚洲国产精品久久电影欧美 | 欧美成人精品福利在线视频| 亚洲黄片高清在线观看| 国产精品亚洲综合久久| 成人国产一区二区三区| 国产 精品 自在 线免费| 国产午夜成人av在线播放| 曰批免费视频播放免费直播| 中文字幕影片免费在线观看| 久久AV中文综合一区二区| 国产自产在线视频一区| 久久午夜av一区二区| 国产精品无码一区二区三区电影 | 亚洲专区一区二区在线观看| 亚洲一区二区高清在线| 亚洲精品视频一区二区三区四区| 老女老肥熟女一区二区| 久久国产免费观看精品3| a级毛片100部免费看| 成年女人永久免费看片| 人妻少妇中文字幕久久69堂| 精品国产一区二区三区九一色| 精品人妻系列无码人妻漫画| 无码国产69精品久久久久孕妇| 国产无遮挡又黄又爽在线视频| 国产成人久久综合热| 亚洲高清中文字幕精品不卡| 国产av一区二区内射| 久久亚洲道色综合久久| 内射欧美老妇wbb| 免费看久久妇女高潮a| 福利一区二区三区视频午夜观看| 国产精品一区二区午夜久久| 日本高清成人一区二区三区| 小雪好紧好滑好湿好爽视频| 爱情岛论坛亚洲品质自拍hd| 欧美丝袜秘书在线一区| 女人天堂av人禽交在线观看| 男人激烈吮乳吃奶视频免费|