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

    一種基于LCD—Hilbert變換和奇異譜熵的配電網(wǎng)暫時過電壓類型識別方法

    2018-05-14 13:31:20金濤許立彬高偉郭謀發(fā)陳永往
    電機與控制學(xué)報 2018年11期

    金濤 許立彬 高偉 郭謀發(fā) 陳永往

    摘要:針對當(dāng)前電力系統(tǒng)配網(wǎng)領(lǐng)域的過電壓識別問題,提出一種基于時頻分布特征的配電網(wǎng)暫時過電壓分類識別方法。在暫時過電壓中,單相金屬性接地過電壓的三相電壓能量分布均勻度最低,間歇性弧光接地零序電壓的直流分量最高,鐵磁諧振零序電壓信號在幅值和頻率集中頻帶上存在很大差異。計算零序電壓能量貢獻率、零序電壓平均值和三相電壓奇異譜熵,提取過電壓信號的時域能量分布特征;采用局部特征尺度分解(LCD)和Hilbert變換結(jié)合帶通濾波算法計算零序電壓重心頻帶,提取過電壓信號的頻域能量分布特征;結(jié)合閾值判別法實現(xiàn)暫時過電壓類型識別。該方法不需要分類器,算法簡單,計算時間少。仿真和實驗表明該識別方法在不同故障條件下均有更高的識別率。

    關(guān)鍵詞:中性點不接地系統(tǒng);暫時過電壓;能量貢獻率;平均值;奇異譜熵;重心頻帶

    DOI:10.15938/j.emc.2018.11.000

    中圖分類號:TM 73

    文獻標(biāo)志碼:A

    文章編號:1007-449X(2018)11-0000-00

    0引言

    運行經(jīng)驗表明,過電壓是影響配電網(wǎng)安全運行的重要因素之一[1-2]。暫時過電壓持續(xù)時間較長,易引起設(shè)備絕緣損壞,從而引發(fā)各類短路故障,危及配電網(wǎng)供電可靠性。因此,及時檢測配電網(wǎng)出現(xiàn)的暫時過電壓,準(zhǔn)確區(qū)分過電壓類型,對于配電網(wǎng)災(zāi)害預(yù)防和故障分析具有重要意義。

    國內(nèi)外普遍采用的過電壓識別方法是先提取特征量,再結(jié)合支持向量機、極端學(xué)習(xí)機、徑向基核函數(shù)(radial basis function, RBF)神經(jīng)網(wǎng)絡(luò)等分類算法,得到識別結(jié)果。文獻[3]采用Morlet小波方法進行故障特征提取,但對于不同種類的過電壓,小波方法存在最優(yōu)小波基選擇困難的問題。文獻[4]采用S變換進行信號分解,通過核主成分分析法提取鐵磁諧振過電壓的特征量結(jié)合支持向量機進行鐵磁諧振的識別。文獻[5]采用經(jīng)驗?zāi)B(tài)分解提取電弧故障特征,并構(gòu)建極端學(xué)習(xí)機,用于配網(wǎng)不同負(fù)載電弧故障的識別,但經(jīng)驗?zāi)B(tài)分解存在模態(tài)混疊和端點效應(yīng)的問題,對于包含不同時間尺度的過電壓信號,分解后的分量不唯一,該方法的適應(yīng)性有待提高。文獻[6]和文獻[7]采用RBF神經(jīng)網(wǎng)絡(luò)作為分類方法,神經(jīng)網(wǎng)絡(luò)對樣本數(shù)量要求較高,且RBF神經(jīng)網(wǎng)絡(luò)存在局部極小值問題。文獻[8]建立過電壓分層識別系統(tǒng),采用多級支持向量機對特征量逐級選擇,應(yīng)用于暫時過電壓類型識別,但多級支持向量機需要大量訓(xùn)練樣本,且訓(xùn)練周期長。文獻[9]結(jié)合傅里葉變換和小波變換提取三類比值作為特征參量對鐵磁諧振過電壓信號進行分類識別。文獻[10]利用希爾伯特黃變換(HilbertHuang transform, HHT)作為過電壓信號的時頻分析工具,基于瞬時幅值譜、Hilbert邊際譜和Hilbert時頻譜提取特征量,對于某些過電壓信號,分解時同樣會出現(xiàn)端點效應(yīng)、模態(tài)混疊的問題。文獻[11]和文獻[12]采用數(shù)學(xué)形態(tài)學(xué)濾波算子提取零序過電壓信號波形的突變特征,作為識別弧光接地過電壓的判據(jù),所提方法較為新穎,但不同的信號采樣頻率和干擾信號會影響結(jié)構(gòu)元素的選取,進而影響信號處理結(jié)果[13]。文獻[14]基于S變換模矩陣奇異值結(jié)合分類器識別操作及雷電過電壓,在高頻瞬態(tài)過電壓信號的識別上具有一定的參考價值。局部特征尺度分解(local characteristicscale decomposition,LCD)是一種新的非線性、非平穩(wěn)信號的自適應(yīng)時頻分析方法。該方法繼承了經(jīng)驗?zāi)B(tài)分解法(empirical mode decomposition,EMD)的優(yōu)點,并改善了EMD存在的端點效應(yīng)、模態(tài)混疊和分解速度等問題,已廣泛應(yīng)用于機械故障診斷領(lǐng)域[15]。

    本文采用LCDHilbert變換結(jié)合帶通濾波算法提取零序電壓信號的頻域分布特征,將其與三相電壓、零序電壓的時域分布特征有機結(jié)合,采用閾值判別法對配電網(wǎng)暫時過電壓進行類型識別。該方法不需要使用分類器,可有效克服樣本獲取困難、耗時長等問題,且在迭代次數(shù)上優(yōu)于EMD,用于暫時過電壓的實時類型識別具有良好的效果。

    1配電網(wǎng)暫時過電壓特性分析

    配電網(wǎng)暫時過電壓包括諧振過電壓和工頻電壓升高?;诖耍匮芯恐行渣c不接地系統(tǒng)中出現(xiàn)的對絕緣危害較大的不對稱接地故障引起的工頻電壓升高和電磁式電壓互感器飽和引起的鐵磁諧振過電壓。下面分析過渡電阻阻值大小對不對稱接地故障過電壓倍數(shù)的影響。圖1為單相接地故障的復(fù)合序網(wǎng)圖,Rf為過渡電阻,Uf|0|為f點正常時的電壓;Z∑(1)、Z∑(2)、Z∑(0)分別對應(yīng)正序網(wǎng)、負(fù)序網(wǎng)、零序網(wǎng)f點的等值阻抗。

    設(shè)故障相為A相,EA=EA∠0°為故障前A相的正常電壓,則EA=Uf|0|,計算故障后A相電壓UA、B相電壓UB為:

    UA=3RfZ∑(1)+Z∑(2)+Z∑(0)+3RfEA,(1)

    UB=a2Uf(1)+aUf(2)+Uf(0)=

    a2(EA-Z∑(1)If(1))+a(-Z∑(2)If(2))+

    (-Z∑(0)If(0))。(2)

    式中:a為旋轉(zhuǎn)因子,a=1∠120°,a2=1∠240°,a2+a+1=0。X∑(1)、X∑(2)和X∑(0)分別為正序、負(fù)序和零序網(wǎng)的等效電抗。因為在實際系統(tǒng)中,X∑(i)R∑(i)(i=1,2,0),因此忽略系統(tǒng)的三序電阻,則X∑(i)≈Z∑(i)(i=1,2,0),且X∑(1)≈X∑(2),If(1)=If(2)=If(0)=Uf|0|Z∑(1)+Z∑(2)+Z∑(0)+3Rf,則式(2)可變?yōu)?/p>

    UB=a2EA-j(X∑(0)-X∑(1))If(1)=

    a2EA-j(X∑(0)-X∑(1))EAj(2X∑(1)+X∑(0))+3Rf。(3)

    令k0=X∑(0)X∑(1),則式(3)可變?yōu)?/p>

    UB=EB-j(k0-1)EAj(k0+2)+3RfX∑(1)=

    EB-1-j3+3RfX∑(1)j(k0+2)+3RfX∑(1)EA。(4)

    當(dāng)Rf=0時,UA=0,由于X∑(0)X∑(1),因此k0=

    SymboleB@ ,UB=EB-EA=3EB∠-30°;當(dāng)Rf阻值增大至與零序電抗X∑(0)數(shù)量級相同時,kf=RfX∑(1)=∞,則式(4)可變?yōu)?/p>

    UB=EB-1-3jk0kf+3EA。(5)

    Rf越大,kf越大,|UB|越接近|EB|,因此在高阻接地故障時非故障相電壓UB的過電壓倍數(shù)很低,故僅考慮單相金屬性接地的引起的工頻電壓升高。

    圖2和圖3是典型的過電壓三相電壓和零序電壓時域波形,過電壓發(fā)生時刻均為0.01 s,采樣頻率為20 kHz,以故障前半個工頻周波為起始時刻,給出5個工頻周期的過電壓波形示例。

    單相金屬性接地過電壓發(fā)生時,如圖2(a)所示,故障相電壓幅值跌落至0,非故障相電壓幅值升高至原線電壓幅值,因此相電壓之間能量分布相差很大,盡管故障點、故障時刻和過渡電阻會發(fā)生變化,這個特征仍保持不變,因此可作為識別單相金屬性接地過電壓的依據(jù)?;」饨拥剡^電壓的三相波形如圖2(b)所示,當(dāng)C相電壓達(dá)到最大值時(0.01 s)發(fā)生第一次絕緣擊穿,故障點的電弧在工頻電流通過第一個過零點時熄滅,此時,非故障相上的自由電荷將沿三相對地電容重新分布,于是在各相上便產(chǎn)生了同等的位移電壓,熄弧時段三相電壓恢復(fù)為正弦電壓。給出其零序電壓波形如圖2(d)所示,由于位移電壓的存在,零序電壓每隔半個周波出現(xiàn)一次直流分量,可作為識別弧光接地過電壓的時域特征。

    鐵磁諧振過電壓以單相故障消失作為“激發(fā)”條件,0.01 s時刻單相接地故障消失,原故障相電壓振幅逐漸增大,三相電壓的能量分布較為均勻,但不同類型的諧振過電壓振蕩頻率不同(如圖3(a)、圖3(b)、圖3(e)所示):分頻諧振可能包含1/2次、1/3次、1/5次等分次諧波含量,振蕩頻率為工頻的分?jǐn)?shù)倍;高頻諧振的振蕩頻率為工頻的整數(shù)次倍;中性點位移電壓屬于基頻零序電壓時,產(chǎn)生基頻諧振,顯然其基波含量占比最大。鐵磁諧振的振蕩頻率不同必然會反映在零序電壓上,為提高識別速率,可將零序電壓能量集中的頻帶作為區(qū)分鐵磁諧振過電壓的特征量。

    2配電網(wǎng)暫時過電壓識別方法分析

    根據(jù)以上分析,提出一種基于LCDHilbert變換和奇異譜熵的過電壓識別方法,對典型的暫時過電壓進行分層識別。如圖4所示,給出涉及到的幾類過電壓的分層結(jié)構(gòu)圖。間歇性弧光過電壓雖屬操作過電壓,但因電弧間歇性熄滅和重燃,持續(xù)時間較長,故將間歇性弧光接地過電壓歸為暫時過電壓一起分析。

    首先,基于零序電壓在時域上的能量分布差異,采用能量貢獻率作為特征量,排除由斷路器操作引發(fā)的過電壓,從而獲得本文的識別對象,如圖4虛線框內(nèi)所示;第二,基于三相電壓能量分布差異,采用奇異譜熵作為特征量,識別單相金屬性接地;第三,基于零序電壓直流分量的高低,采用平均值的絕對值作為特征量,識別間歇性弧光接地,同時,基于鐵磁諧振零序電壓能量集中頻帶的不同,采用重心頻帶作為特征量,識別高頻諧振、基頻諧振和分頻諧振。下面介紹地方法的原理。

    2.1時域特征量提取

    在分析過程中,操作過電壓是配電網(wǎng)中較為常見的故障,在分析過程中需要先行排除。操作過電壓持續(xù)時間較短,一般小于10 ms,暫時過電壓持續(xù)時間較長,從幾個工頻周期到1.5~2 h不等。在所截取的時段內(nèi),操作過電壓的能量分布多集中于故障后0.5個周波,而暫時過電壓能量分布較均勻,因此,暫時過電壓和操作過電壓的零序電壓信號在故障后的0.5個周波的能量與所截取時段的總能量的比值相差較大。定義能量貢獻率E為

    E=∑N1k=1v0(k)2/∑N2k=1v0(k)2。(6)

    式中:N1為0.5個周波的采樣點數(shù);N2為截取的總時段的采樣點數(shù);v0(k)為零序電壓信號序列。

    “熵”最早出現(xiàn)在熱力學(xué)研究中,是用于衡量一個系統(tǒng)中分子狀態(tài)的均勻狀況,分子狀態(tài)越均勻,熵值越高,當(dāng)系統(tǒng)中分子分布最均勻時,熵值達(dá)到最大[16]?;谙嚯妷褐g的能量分布差異,由于接地瞬間單相金屬性接地三相電壓可能有高頻振蕩,因此以故障發(fā)生后2個工頻周波為起始點,截取5個工頻周波的三相電壓信號,計算其奇異譜熵作為識別單相金屬性接地過電壓的特征量。

    設(shè)截取的信號采樣點數(shù)為n,三相電壓信號矩陣為U=[ua,ub,uc]T,對矩陣U的每一行作奇異值分解,可得奇異譜Λ=diag(μ1,μ2,μ3),三相電壓信號奇異譜熵定義為

    S=-∑3i=1qilogqi。(7)

    式中:qi=μi/∑3i=1μi(i=1,2,3)表征了μi對整個奇異譜的貢獻率[17];S可有效表征三相電壓信號奇異值分布的不均勻度,進而反映三相電壓之間幅值分布的差異。

    間歇性弧光接地的零序電壓包含較高的直流分量,采用零序電壓信號的平均值絕對值表征這個特性。零序電壓信號的平均值絕對值M定義為

    M=|∑nk=1v0(k)|n。(8)

    式中:v0(k)是各種過電壓零序電壓信號;n為故障時刻開始5個周波的零序電壓信號的采樣點數(shù)。

    2.2LCDHILBERT變換

    LCDHilbert變換由兩部分組成:首先,通過LCD將原信號分解成有限個內(nèi)稟尺度分量(intrinsic scale component,ISC);其次,對各個內(nèi)稟尺度分量進行Hilbert變換,進而得到瞬時頻率函數(shù)。

    LCD分解計算得到的ISC應(yīng)滿足以下條件:1)信號數(shù)據(jù)的任意兩個相鄰極值點符號互異;2)信號數(shù)據(jù)的任意兩個相鄰極大(或?。┲迭c確定的直線在二者之間的極?。ɑ虼螅┲迭c所對應(yīng)的橫坐標(biāo)處的函數(shù)值與該極?。ɑ虼螅┲迭c的比值保持不變[14]。

    根據(jù)上述兩個條件,將過電壓信號進行局部特征尺度分解,分解成若干個ISC和一個殘余項之和為

    x(t)=∑ni=1ISCi(t)+r(t)。(9)

    式中:x(t)為原始信號;n為ISC的個數(shù);ISCi(t)為第i個ISC;r(t)為剩余分量?;l諧振的零序電壓信號如圖1所示,其經(jīng)LCD分解后的ISC如圖5所示。

    2.3頻域特征量提取

    利用LCDHilbert變換結(jié)合帶通濾波算法可實現(xiàn)過電壓信號在頻帶上的分離并濾去高頻噪聲干擾,增強信號可比性,提高識別成功率。帶通濾波算法的基本原理是:對故障波形進行LCDHilbert變換,得到ISC分量和瞬時頻率函數(shù)之后,根據(jù)瞬時頻率范圍選擇合適的頻帶寬度Δf,以此進行頻帶劃分,則第i個頻帶的頻率范圍為[(i-1)Δf ,iΔf],將所有ISC分量不在該頻率范圍內(nèi)的瞬時能量點都置為零,其余在該頻率范圍內(nèi)的瞬時能量點保持不變,即可得到所有ISC分量在該頻帶的成分。最后將經(jīng)第i個頻帶濾波后的所有ISC分量疊加,可得到原始信號在該頻帶內(nèi)的分量,同理可得到其他頻帶內(nèi)的分量。

    不同頻帶重構(gòu)波形的能量大小可反映過電壓信號的頻域分布情況。以各頻帶重構(gòu)信號的平方和表征該頻帶信號能量,第i頻帶的能量為

    ei=∑mk=1c2i(k)。(15)

    式中:ei為第i頻帶的信號能量;m為重構(gòu)信號的數(shù)據(jù)長度;ci(k)為第i頻帶的重構(gòu)信號。

    物理學(xué)中用物理重心來反映物體受重力作用所集中到的點,將重心的概念應(yīng)用到頻域分析上,選取各個頻帶中能量最大的頻帶作為頻域的“重心”,可反映波形諧波含量的集中位置。定義重心頻帶為

    Ng=imax。(16)

    式中:Ng為重心頻帶;imax為能量最大的頻帶編號。

    2.4方法的計算流程

    確定了特征量E、S、M和Ng后,針對中性點不接地系統(tǒng),給出配電網(wǎng)暫時過電壓類型的識別流程如圖6所示。暫時過電壓類型識別步驟如下:

    1)獲取配電網(wǎng)母線處的三相電壓和零序電壓信號,對三相電壓序列求取奇異譜熵S,判斷其是否小于閾值,若是,判為單相金屬性接地過電壓;

    2)若奇異譜熵S大于設(shè)定的閾值,求取零序電壓序列平均值的絕對值M,判斷其是否大于閾值,若是,判為間歇性弧光接地過電壓;

    3)若M小于設(shè)定的閾值,采用LCDHilbert變換結(jié)合帶通濾波算法重構(gòu)相應(yīng)的零序電壓波形,求取重心頻帶Ng,按照劃定好的范圍進行分類。

    3仿真實驗分析和討論

    根據(jù)以上分析,采用ATPEMTP電磁暫態(tài)仿真軟件搭建中性點不接地配網(wǎng)模型,如圖7所示。其中:G為110 kV無窮大系統(tǒng)電源;T1為110 kV/10.5 kV主變壓器,T2為10 kV/0.4 kV配電變壓器;負(fù)荷采用感性負(fù)載R+jX=(20+j40) Ω等值替代。饋線共有全架空線路(OL)、全電纜線路(CL)和線-纜混合線路3種類型,線路參數(shù)如表1所示。K1為時控開關(guān),用于模擬接地故障,Rf為過渡電阻。母線PT為電磁式電壓互感器,以“F”開頭的編號均為仿真實驗設(shè)置的故障點。

    根據(jù)工頻熄弧理論分析所得的弧光接地過電壓倍數(shù)較接近配電網(wǎng)出現(xiàn)的弧光接地的過電壓倍數(shù),該理論認(rèn)為故障相每隔一個工頻周期電弧就會熄滅和重燃一次,并且假定電弧在工頻電流過零時熄滅,達(dá)到最大恢復(fù)電壓時重燃[12],因此基于該理論,通過周期性投切時控開關(guān)K1進行弧光接地過電壓的模擬。采用單相接地故障消失作為鐵磁諧振的“激發(fā)”條件,為了實現(xiàn)鐵磁諧振的參數(shù)匹配要求[18],通過投切不同長度的饋線形成不同的匹配電容,達(dá)到改變諧振回路三相導(dǎo)納的目的,以產(chǎn)生基頻諧振、分頻諧振和高頻諧振。圖7中,投入饋線L3、饋線L6進行分頻諧振的仿真,投入饋線L5、饋線L6進行基頻諧振的仿真,投入饋線L6進行高頻諧振的仿真。其中線路的參數(shù)如表1所示。

    將過渡電阻設(shè)置為1、2、5、10、20 Ω,故障初相角選擇±90°(對應(yīng)正負(fù)電壓峰值時刻,若分析對象為鐵磁諧振,則是指單相接地故障消除時刻初始相角選擇±90°),每種過電壓記錄10組波形,計算其奇異譜熵。如圖9所示,每組數(shù)據(jù)的1~5個點對應(yīng)-90°初相角、過渡電阻分別為1~20 Ω時的奇異譜熵,第6~10個點對應(yīng)90°初相角、過渡電阻分別為1~20 Ω時的奇異譜熵??梢钥闯觯簡蜗嘟饘傩越拥剡^電壓的奇異譜熵最低,均小于1.2;其余四類過電壓的奇異譜熵均大于1.5,隨著過渡電阻的增大,單相金屬性接地的奇異譜熵呈上升趨勢。計算在20 Ω過渡電阻情況下,不同故障點、不同故障初相角的100組三相電壓數(shù)據(jù)的奇異譜熵,判別單相金屬性接地過電壓的閾值可取1.4。

    弧光接地熄弧時段零序電壓存在不為零直流分量,計算弧光接地與鐵磁諧振的零序電壓平均值如圖10所示。可以看出,弧光接地的零序電壓平均值的絕對值最大。因此,采用零序電壓平均值的絕對值M作為識別弧光接地過電壓的特征量,計算不同過渡電阻、不同故障點的100組零序電壓數(shù)據(jù)的M值,確定其閾值為2 000。

    以鐵磁諧振的零序電壓作為分析對象,采用帶通濾波算法對LCD分解獲得的各個ISC依據(jù)Hilbert瞬時頻率函數(shù)分頻帶多支重構(gòu):通過式(10)~式(14)計算不同過渡電阻、不同故障初相角等因素影響下的5個工頻周波的零序電壓波形的瞬時頻率向量,共200組數(shù)據(jù)的仿真結(jié)果,發(fā)現(xiàn)信號的瞬時頻率集中在0~600 Hz范圍內(nèi),因此將此范圍作為帶通濾波算法的總頻帶寬度;將總頻帶劃分為6個頻帶:0~10 Hz、10~20 Hz、20~30 Hz、30~40 Hz、40~80 Hz、80~600 Hz,對各頻帶依次編號,0~10 Hz分量為頻帶1,10~20 Hz為頻帶2,以此類推;設(shè)當(dāng)前濾波頻帶為i,若瞬時頻率不在頻帶頻率范圍內(nèi),則將該點對應(yīng)的ISC瞬時能量置零,最后將經(jīng)第i頻帶濾波后的所有ISC疊加,可得原始信號在該頻帶內(nèi)的分量,同理可得到其他頻帶內(nèi)的分量?;l諧振零序電壓重構(gòu)后的波形如圖11所示,幅值最高的重構(gòu)信號出現(xiàn)在基波頻帶,與基頻諧振的原理分析一致。

    計算不同過渡電阻、不同故障初相角影響下鐵磁諧振零序電壓的重心頻帶,如表2所示??梢钥闯觯悍诸l諧振重心頻帶為2或3,這是因為分頻諧振1/2次或1/3次諧波含量最高;基頻諧振重心頻帶為5,這是因為基頻諧振的基波含量最高;高頻諧振重心頻帶為6,這是因為高頻諧振的高次諧波含量最高。因此重心頻帶Ng可作為識別分頻諧振、基頻諧振、高頻諧振的特征量。判別高頻諧振的Ng取為6,判別基頻諧振的Ng取為5,判別分頻諧振的Ng閾值取為4。

    通過改變故障點、故障相別、故障初相角、過渡電阻4個因素,對圖7所示的中性點不接地系統(tǒng)發(fā)生的暫時過電壓進行仿真,以驗證識別方法的有效性,仿真數(shù)據(jù)的實驗條件如表3所示。同時給出其余仿真數(shù)據(jù)的識別統(tǒng)計結(jié)果,如表4所示。將LCD應(yīng)用于所提方法對于暫時過電壓類型識別的正確率可達(dá)到94.98%,將EMD應(yīng)用于所提方法對暫時過電壓類型識別也是有效的,準(zhǔn)確率達(dá)92%以上,但是低于LCD方法,原因在于其混疊效應(yīng)較為嚴(yán)重,而將頻帶范圍為40~80 Hz的基頻諧振誤識別為30~40 Hz的分頻諧振。采用LCD方法的識別速度總體優(yōu)于EMD方法。所提方法不僅可以實現(xiàn)鐵磁諧振的分類,還可獲取諧振頻率范圍,對于鐵磁諧振的抑制具有參考價值。

    為了模擬現(xiàn)場電磁環(huán)境的噪聲干擾,抽取表4中已正確識別的若干組過電壓數(shù)據(jù)加入20 dB的高斯白噪聲,驗證噪聲對識別方法的影響。識別結(jié)果如表5所示,識別正確率為100%,這是由于帶通濾波算法本身具有良好的濾波效果,因此噪聲對于所提方法的影響十分微弱,說明了所提方法具有良好的抗噪性能。

    本文識別方法在表4所示的1 434組數(shù)據(jù)中的平均識別時間是0.13 s。在實驗室進行算法的工程驗證,采用數(shù)據(jù)采集卡作為在線采集終端,在福州大學(xué)配電網(wǎng)物理仿真系統(tǒng)[19]上進行接地性實驗,處理器是電腦主機,數(shù)據(jù)采集卡內(nèi)嵌在電腦主機中,識別算法可以在線識別。

    在實際電力現(xiàn)場應(yīng)用中,已在福建省數(shù)個變電站安裝了10 kV過電壓在線采集終端,其實物圖如圖12所示,其主要功能是過電壓啟動后,將過電壓波形通過GPRS通信上傳到后臺主站,由主站負(fù)責(zé)離線分析,當(dāng)前還無法應(yīng)用于實時控制。

    4結(jié)論

    操作過電壓的持續(xù)時間較暫時過電壓短,單相金屬性接地過電壓的三相能量分布差異較其余四種暫時過電壓大,間歇性弧光接地的零序電壓包含最高的零序分量,鐵磁諧振過電壓的能量集中頻帶均不相同。采用零序電壓能量貢獻率、零序電壓平均值、三相電壓奇異譜熵作為時域特征量,通過LCDHilbert變換和帶通濾波算法重構(gòu)零序電壓波形,以重心頻帶作為頻域特征量,結(jié)合閾值判別法構(gòu)成暫時過電壓的類型識別方法。該方法所提取的特征量,綜合利用了零序電壓和三相電壓信號在頻域和時域的整體分布規(guī)律,完整地刻畫了過電壓信號特征,對不同的故障點、故障相別、故障過渡電阻、故障初相角情況下發(fā)生的暫時過電壓均能正確識別。

    參 考 文 獻:

    [HT6SS]

    [1]黃艷玲,司馬文霞,楊慶,等. 電力系統(tǒng)實測過電壓信號的特征量提取與驗證[J].高電壓技術(shù),2013,39(1):60.

    HUANG Yanling, SIMA Wenxia, YANG Qing,et al. Features extraction and validation from measured overvoltage signals of power system[J]. High Voltage Engineering, 2013, 39(1): 60.

    [2]CERRETTI A, GATTA M F, GERI A, et al. Ground fault temporary overvoltages in MV networks: evaluation and experimental tests[J]. IEEE Transaction Power Delivery, 2012, 27(3): 1592.

    [3]侯新國, 牛超, 楊忠林. 基于最優(yōu)Morlet小波自適應(yīng)包絡(luò)解調(diào)的弱故障特征提取方法[J]. 電機與控制學(xué)報, 2016, 20(10): 88.

    HOU Xinguo, NIU Chao, YANG Zhonglin. Method to extract weak fault feature based on optimal Morlet wavelet adaptive envelope demodulation[J]. Electric Machines and Control, 2016, 20(10): 88.

    [4]MOKRYANI G, SIANO P, PICCOLO A. Identification of ferroresonance based on Stransform and support vector machine[J]. Simulation Modelling Practice and Theory, 2010, 18(9): 1412.

    [5]張麗萍, 繆希仁, 石敦義. 基于EMD和ELM的低壓電弧故障識別方法的研究[J]. 電機與控制學(xué)報, 2016, 20(9): 54.

    ZHANG Liping, MIAO Xiren, SHI Dunyi. Research on low voltage arc fault recognition method based on EMD and ELM [J]. Electric Machines and Control, 2016, 20(9): 54.

    [6]管春,周雒維,盧偉國.基于多標(biāo)簽RBF神經(jīng)網(wǎng)絡(luò)的電能質(zhì)量復(fù)合擾動分類方法[J].電工技術(shù)學(xué)報,2011(8):198.

    GUAN Chun, ZHOU Luowei, LU Weiguo. Recognition of multiple power quality disturbances using multilabel RBF neural networks [J].Transactions of China Electrotechnical Society, 2011(8): 198.

    [7]司馬文霞,王荊,楊慶,等.HilbertHuang變換在電力系統(tǒng)過電壓識別中的應(yīng)用[J].高電壓技術(shù),2010,36(6):1480.

    SIMA Wenxia, WANG Jing, YANG Qing, et al. Application of HilbertHuang transform to power system overvoltage recognition[J]. High Voltage Engineering, 2010, 36(6): 1480.

    [8]杜林,李欣,王麗蓉,等.電力系統(tǒng)暫時過電壓多級支持向量機分層識別[J].電力系統(tǒng)保護與控制,2012,40(4):26.

    DU Lin, LI Xin, WANG Lirong, et al. Temporary overvoltage layered pattern identification based on multistage support vector machine in powersystem[J]. Power System Protection and Control, 2012, 40(4): 26.

    [9]杜林,李欣,吳高林,等. 采用3類特征參量比值法的鐵磁諧振過電壓識別[J].高電壓技術(shù),2011,37(9): 2241.

    DU Lin,LI Xin, WU Gaolin, et al. Ferroresonance overvoltage identification using three feature parameters of ratio method[J]. High Voltage Engineering, 2011, 37(9): 2241.

    [10]WANG Jing, YANG Qing, CHEN Lin, et al. Classification and identification of overvoltage based on HHT and SVM[J]. High Voltage Engineering, 2012, 38(8): 2068.

    [11]司馬文霞,冉銳,袁濤,等.采用數(shù)學(xué)形態(tài)學(xué)的弧光接地過電壓識別方法[J].高電壓技術(shù),2010,36(4):835.

    SIMA Wenxia, RAN Rui, YUAN Tao, et al. Identification of arc grounding overvoltage using mathematical morphology transform[J]. High Voltage Engineering, 2010, 36(4): 835.

    [12]羅寧昭,張曉鋒,劉佳.基于GradientClosing變換的接地故障識別[J].電力系統(tǒng)保護與控制,2011,39(23):74.

    LUO Ningzhao, ZHANG Xiaofeng, LIU Jia. Ground fault identification based on GradientClosing transform[J]. Power System Protection and Control, 2011, 39(23): 74.

    [13]王健民, 張春東. 電網(wǎng)諧波檢測的數(shù)學(xué)形態(tài)學(xué)優(yōu)化算法[J]. 哈爾濱理工大學(xué)學(xué)報, 2013, 18(2):114.

    WANG Jianmin, ZHANG Chundong. The mathematical morphology optimization algorithm of harmonic detection[J]. Journal of Harbin University of Science and Technology, 2013, 18(2):114.

    [14]杜林,李欣,司馬文霞,等. S變換模矩陣和最小二乘SVM在雷電及操作過電壓識別中的應(yīng)用[J].電力自動化設(shè)備,2012,32(8):35.

    DU Lin, LI Xin,SIMA Wenxia, et al. Application of Stransform modular matrix and LSSVM in identification of lightning and switching overvoltages[J]. Electric Power Automation Equipment, 2012,32(8):35.

    [15]程軍圣,鄭近德,楊宇.一種新的非平穩(wěn)信號分析方法——局部特征尺度分解法[J].振動工程學(xué)報,2012,25(2): 215.

    CHENG Junsheng, ZHENG Jinde, YANG Yu. A nonstationary signal analysis approach: the local characteristic scale decomposition method[J]. Journal of Vibration Engineering, 2012, 25(2): 215.

    [16]陳仕龍,曹蕊蕊,畢貴紅,等.利用多分辨奇異譜熵和支持向量機的特高壓直流輸電線路區(qū)內(nèi)外故障識別方法[J].電網(wǎng)技術(shù),2015,39(4): 991.

    CHEN Shilong, CAO Ruirui, BI Guihong, et al. Distinguish faults located inside/outside protection zone of UHVDC transmission line by multiresolution singular spectrum entropy and support vector machine[J]. Power System Technology, 2015, 39(4): 991.

    [17]楊豐源,宋輝,程序,等.基于多分辨高階奇異譜熵分析的局部放電信號特征提取[J].電網(wǎng)技術(shù),2016, 40(10):3266.

    YANG Fengyuan, SONG Hui, CHENG Xu, et al. Feature extraction of partial discharge signal based on multiresolution and highorder singular spectrum entropy analysis[J]. Power System Technology, 2016, 40(10): 3266.

    [18]司馬文霞,郭飛,楊慶,等.變電站鐵磁諧振仿真分析及抑制措施研究[J].電力自動化設(shè)備,2007,27(6): 23.

    SIMA Wenxia, GUO Fei, YANG Qing, et al. Simulative analysis of substation ferroresonance and its suppression measures[J]. Electric Power Automation Equipment, 2007, 27(6): 23.

    [19]郭謀發(fā),游林旭,洪翠,等.基于LCDHilbert譜奇異值和多級支持向量機的配電網(wǎng)故障識別方法[J].高電壓技術(shù),2017,43(4):1246.

    GUO Moufa, YOU Linxu, HONG Cui, et al. Identification method of distribution network faults based on singular value of LCDHilbert spectrums and multilevel SVM[J]. High Voltage Engineering, 2017, 43(4):1246.

    (編輯:邱赫男)

    国产精品熟女久久久久浪| 青春草亚洲视频在线观看| 男女国产视频网站| 国产精品久久久久成人av| 婷婷成人精品国产| 考比视频在线观看| 成人毛片a级毛片在线播放| 亚洲欧美成人综合另类久久久| 男女边吃奶边做爰视频| 91精品伊人久久大香线蕉| 热re99久久国产66热| 一本大道久久a久久精品| 久热久热在线精品观看| 国产精品偷伦视频观看了| 一区二区av电影网| 国产成人a∨麻豆精品| 日本黄色片子视频| 久久亚洲国产成人精品v| 99精国产麻豆久久婷婷| 日日撸夜夜添| 久久韩国三级中文字幕| 成年av动漫网址| 一级a做视频免费观看| 好男人视频免费观看在线| 热re99久久精品国产66热6| 考比视频在线观看| 久久99蜜桃精品久久| 男女边吃奶边做爰视频| 国产乱来视频区| 午夜视频国产福利| 国产极品天堂在线| 色哟哟·www| 日本欧美国产在线视频| 一级毛片aaaaaa免费看小| 免费观看在线日韩| 亚洲综合色网址| 中文字幕精品免费在线观看视频 | 国产探花极品一区二区| 人人妻人人澡人人看| 日韩成人av中文字幕在线观看| 久久久精品区二区三区| 欧美人与善性xxx| 一区二区三区精品91| 一级毛片电影观看| 99热这里只有是精品在线观看| 赤兔流量卡办理| 亚洲激情五月婷婷啪啪| 亚洲国产精品一区三区| 妹子高潮喷水视频| 久久久久国产网址| 亚洲av男天堂| 大香蕉97超碰在线| 国产高清三级在线| 美女大奶头黄色视频| 啦啦啦在线观看免费高清www| 插阴视频在线观看视频| 夜夜爽夜夜爽视频| 纵有疾风起免费观看全集完整版| 一级片'在线观看视频| 免费日韩欧美在线观看| 国产精品国产av在线观看| 交换朋友夫妻互换小说| 多毛熟女@视频| 在线 av 中文字幕| 美女午夜性视频免费| 亚洲av第一区精品v没综合| 国产精品亚洲av一区麻豆| 欧美乱码精品一区二区三区| videos熟女内射| 人人妻人人澡人人爽人人夜夜| 黄色片一级片一级黄色片| xxxhd国产人妻xxx| 久久精品国产a三级三级三级| videosex国产| 少妇的丰满在线观看| 午夜福利影视在线免费观看| 欧美av亚洲av综合av国产av| 91av网站免费观看| 亚洲综合色网址| 曰老女人黄片| 欧美黑人欧美精品刺激| 青草久久国产| 久久人人爽av亚洲精品天堂| 飞空精品影院首页| 国产av精品麻豆| 亚洲成a人片在线一区二区| 国产成+人综合+亚洲专区| 国产男靠女视频免费网站| 中文字幕人妻丝袜制服| av又黄又爽大尺度在线免费看| 亚洲一区中文字幕在线| 久久中文看片网| 久久精品亚洲av国产电影网| 少妇猛男粗大的猛烈进出视频| 伊人久久大香线蕉亚洲五| 国产男女超爽视频在线观看| 午夜激情久久久久久久| 精品国产亚洲在线| 久久中文字幕人妻熟女| 夫妻午夜视频| 免费在线观看影片大全网站| 国产精品九九99| 女警被强在线播放| av不卡在线播放| avwww免费| 亚洲精华国产精华精| 国产日韩一区二区三区精品不卡| 12—13女人毛片做爰片一| 成年人免费黄色播放视频| 免费在线观看黄色视频的| 王馨瑶露胸无遮挡在线观看| avwww免费| 18在线观看网站| 久久人人爽av亚洲精品天堂| 中文字幕制服av| 色94色欧美一区二区| 日韩一卡2卡3卡4卡2021年| 成人永久免费在线观看视频 | 在线播放国产精品三级| 黄色 视频免费看| 国产亚洲欧美在线一区二区| 国产一区二区在线观看av| 一本久久精品| 精品一品国产午夜福利视频| 最黄视频免费看| 丝袜美足系列| 久久毛片免费看一区二区三区| 亚洲五月婷婷丁香| 一本综合久久免费| 精品一区二区三卡| 如日韩欧美国产精品一区二区三区| 国产av国产精品国产| www.自偷自拍.com| 久久亚洲真实| 午夜免费鲁丝| 亚洲国产毛片av蜜桃av| 在线观看免费日韩欧美大片| 亚洲欧美一区二区三区久久| h视频一区二区三区| 黄片大片在线免费观看| 咕卡用的链子| 亚洲成人国产一区在线观看| 亚洲伊人久久精品综合| 五月天丁香电影| 宅男免费午夜| 99热网站在线观看| 亚洲熟女毛片儿| 美女高潮到喷水免费观看| 亚洲成a人片在线一区二区| 欧美日韩av久久| 黄片小视频在线播放| 中文字幕人妻熟女乱码| 激情在线观看视频在线高清 | 成人手机av| 亚洲精品av麻豆狂野| 午夜福利视频在线观看免费| 欧美亚洲 丝袜 人妻 在线| 看免费av毛片| 在线观看免费视频网站a站| 国产黄频视频在线观看| 夜夜骑夜夜射夜夜干| 在线永久观看黄色视频| 国产成人啪精品午夜网站| 999久久久精品免费观看国产| 国产男女超爽视频在线观看| 国产精品欧美亚洲77777| 另类精品久久| 狠狠狠狠99中文字幕| av一本久久久久| 黄色 视频免费看| 中文字幕人妻丝袜制服| 亚洲 欧美一区二区三区| 精品人妻1区二区| 国产成人精品久久二区二区91| 国产日韩一区二区三区精品不卡| 在线观看一区二区三区激情| 亚洲色图 男人天堂 中文字幕| 女人爽到高潮嗷嗷叫在线视频| 午夜福利乱码中文字幕| 后天国语完整版免费观看| 精品高清国产在线一区| av视频免费观看在线观看| 亚洲色图 男人天堂 中文字幕| 一本久久精品| 黄色视频不卡| 欧美大码av| 伊人久久大香线蕉亚洲五| 纵有疾风起免费观看全集完整版| 精品卡一卡二卡四卡免费| 国产又色又爽无遮挡免费看| 国产精品秋霞免费鲁丝片| 色婷婷久久久亚洲欧美| 性高湖久久久久久久久免费观看| 久久久久精品国产欧美久久久| 伊人久久大香线蕉亚洲五| 亚洲黑人精品在线| 色视频在线一区二区三区| 十八禁高潮呻吟视频| 亚洲国产欧美一区二区综合| 啦啦啦视频在线资源免费观看| 日韩视频在线欧美| 欧美成狂野欧美在线观看| 精品久久久精品久久久| 中文字幕人妻丝袜一区二区| 高清av免费在线| 久久久久久免费高清国产稀缺| 久久人妻福利社区极品人妻图片| 国精品久久久久久国模美| 国产精品香港三级国产av潘金莲| 99久久国产精品久久久| 另类精品久久| 精品福利永久在线观看| 亚洲av成人一区二区三| 人人妻人人爽人人添夜夜欢视频| av网站在线播放免费| 国产成人欧美| 丝袜人妻中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 亚洲成av片中文字幕在线观看| 亚洲免费av在线视频| 天堂俺去俺来也www色官网| 香蕉丝袜av| 国产精品一区二区免费欧美| 久久精品熟女亚洲av麻豆精品| 麻豆国产av国片精品| 99国产精品一区二区三区| 久久av网站| 日韩大片免费观看网站| 欧美+亚洲+日韩+国产| 精品福利永久在线观看| 中文字幕精品免费在线观看视频| 亚洲精品美女久久av网站| 大码成人一级视频| 黑人巨大精品欧美一区二区蜜桃| 成人黄色视频免费在线看| www日本在线高清视频| 欧美日本中文国产一区发布| www.精华液| 午夜激情av网站| 99国产精品免费福利视频| av福利片在线| 午夜日韩欧美国产| 一二三四社区在线视频社区8| 欧美日韩亚洲国产一区二区在线观看 | 在线观看免费高清a一片| 欧美黑人精品巨大| 成人国产av品久久久| 亚洲一区二区三区欧美精品| 欧美人与性动交α欧美软件| 国产一区二区三区在线臀色熟女 | 精品欧美一区二区三区在线| 精品亚洲乱码少妇综合久久| 日本一区二区免费在线视频| 精品一区二区三区视频在线观看免费 | 亚洲精品一卡2卡三卡4卡5卡| 黑丝袜美女国产一区| 国产aⅴ精品一区二区三区波| 国产亚洲精品第一综合不卡| 国产麻豆69| 亚洲三区欧美一区| 国产成人精品久久二区二区91| 成人精品一区二区免费| 欧美日本中文国产一区发布| 国产精品一区二区免费欧美| 国产欧美日韩一区二区三| 日韩大片免费观看网站| 精品人妻1区二区| 亚洲熟女精品中文字幕| 亚洲三区欧美一区| cao死你这个sao货| 一本大道久久a久久精品| 国产午夜精品久久久久久| 午夜久久久在线观看| 国产深夜福利视频在线观看| 国产熟女午夜一区二区三区| 黄网站色视频无遮挡免费观看| 亚洲精品中文字幕一二三四区 | 国产成人精品无人区| 在线观看免费高清a一片| 亚洲第一青青草原| 国产高清videossex| 午夜福利视频在线观看免费| 国产麻豆69| 丝袜在线中文字幕| 国产色视频综合| 国产精品久久久久久精品电影小说| 成在线人永久免费视频| 真人做人爱边吃奶动态| 国产在线免费精品| 亚洲精品av麻豆狂野| 亚洲,欧美精品.| 国产精品免费大片| 国产xxxxx性猛交| 久久久精品国产亚洲av高清涩受| 成人av一区二区三区在线看| 两性夫妻黄色片| 一级a爱视频在线免费观看| 久久人妻福利社区极品人妻图片| 丰满少妇做爰视频| 大陆偷拍与自拍| 脱女人内裤的视频| 99精品欧美一区二区三区四区| 国产精品久久电影中文字幕 | 久久天堂一区二区三区四区| 免费av中文字幕在线| 一个人免费看片子| av片东京热男人的天堂| 免费看a级黄色片| 别揉我奶头~嗯~啊~动态视频| 黄网站色视频无遮挡免费观看| av网站免费在线观看视频| 不卡av一区二区三区| 人人妻人人澡人人爽人人夜夜| 久久99热这里只频精品6学生| av视频免费观看在线观看| 欧美成人免费av一区二区三区 | 欧美变态另类bdsm刘玥| 欧美中文综合在线视频| 亚洲国产看品久久| 高清av免费在线| 欧美精品高潮呻吟av久久| aaaaa片日本免费| 国产高清国产精品国产三级| 亚洲精品成人av观看孕妇| 天堂动漫精品| 50天的宝宝边吃奶边哭怎么回事| 精品乱码久久久久久99久播| 欧美大码av| av视频免费观看在线观看| 色综合欧美亚洲国产小说| 涩涩av久久男人的天堂| 首页视频小说图片口味搜索| 欧美av亚洲av综合av国产av| av又黄又爽大尺度在线免费看| 亚洲少妇的诱惑av| 国产成人免费观看mmmm| 日韩成人在线观看一区二区三区| 老司机午夜福利在线观看视频 | 欧美日韩一级在线毛片| 欧美 日韩 精品 国产| 如日韩欧美国产精品一区二区三区| 日韩一区二区三区影片| 中文欧美无线码| 天堂俺去俺来也www色官网| 这个男人来自地球电影免费观看| 亚洲欧美激情在线| 超色免费av| 国产精品98久久久久久宅男小说| 久久精品成人免费网站| 久热爱精品视频在线9| 亚洲一码二码三码区别大吗| 激情在线观看视频在线高清 | 黑人猛操日本美女一级片| 国产xxxxx性猛交| 亚洲国产欧美在线一区| 成年动漫av网址| 美女扒开内裤让男人捅视频| 两个人免费观看高清视频| 国产单亲对白刺激| 一个人免费看片子| 久久久国产一区二区| 国产在线观看jvid| 高清av免费在线| av在线播放免费不卡| 精品久久蜜臀av无| 久久亚洲精品不卡| 国产精品欧美亚洲77777| 国产高清激情床上av| 国产精品 欧美亚洲| 国产精品香港三级国产av潘金莲| 在线播放国产精品三级| 人人妻人人爽人人添夜夜欢视频| 亚洲 欧美一区二区三区| 在线天堂中文资源库| 亚洲色图综合在线观看| 亚洲中文av在线| 成人精品一区二区免费| 欧美日韩亚洲综合一区二区三区_| 天天躁日日躁夜夜躁夜夜| 国产av精品麻豆| 一本久久精品| 纵有疾风起免费观看全集完整版| 在线观看66精品国产| 国产成人影院久久av| 天堂俺去俺来也www色官网| 青草久久国产| 日韩免费av在线播放| 热99re8久久精品国产| av线在线观看网站| 中文字幕色久视频| 日韩欧美一区视频在线观看| 欧美精品av麻豆av| 精品少妇黑人巨大在线播放| 亚洲精品av麻豆狂野| 久久久国产精品麻豆| 黑人操中国人逼视频| 国产精品自产拍在线观看55亚洲 | 在线观看人妻少妇| 精品欧美一区二区三区在线| 久久精品国产亚洲av香蕉五月 | 精品国产超薄肉色丝袜足j| 国产精品偷伦视频观看了| 热re99久久精品国产66热6| 国产欧美日韩一区二区精品| www.自偷自拍.com| 精品久久蜜臀av无| 日韩欧美一区二区三区在线观看 | 精品久久久精品久久久| 亚洲一码二码三码区别大吗| 黑人巨大精品欧美一区二区mp4| 中亚洲国语对白在线视频| 在线 av 中文字幕| 99国产综合亚洲精品| 亚洲精品av麻豆狂野| 久久久水蜜桃国产精品网| 久久青草综合色| 女性生殖器流出的白浆| 高清欧美精品videossex| 国产成人精品无人区| 在线观看66精品国产| 中文欧美无线码| 亚洲国产成人一精品久久久| 国产成人一区二区三区免费视频网站| 涩涩av久久男人的天堂| 午夜久久久在线观看| 蜜桃在线观看..| 午夜免费鲁丝| 亚洲性夜色夜夜综合| 热99国产精品久久久久久7| 精品免费久久久久久久清纯 | 久久毛片免费看一区二区三区| 欧美久久黑人一区二区| 视频区欧美日本亚洲| 亚洲国产欧美网| 欧美日韩亚洲综合一区二区三区_| 欧美精品啪啪一区二区三区| 另类精品久久| 最黄视频免费看| 三上悠亚av全集在线观看| 亚洲一区二区三区欧美精品| a级毛片黄视频| 少妇裸体淫交视频免费看高清 | 国产精品秋霞免费鲁丝片| 丝袜喷水一区| √禁漫天堂资源中文www| 国产不卡av网站在线观看| 免费在线观看日本一区| 久久中文字幕人妻熟女| 国产精品1区2区在线观看. | 成人国语在线视频| 80岁老熟妇乱子伦牲交| 99精品在免费线老司机午夜| 日韩欧美免费精品| 757午夜福利合集在线观看| 视频在线观看一区二区三区| 国产成人欧美| 欧美黄色淫秽网站| 一进一出好大好爽视频| 中文字幕av电影在线播放| 天天添夜夜摸| 9热在线视频观看99| 久久精品亚洲av国产电影网| 丝瓜视频免费看黄片| 91麻豆av在线| 99精品欧美一区二区三区四区| 亚洲精品av麻豆狂野| 香蕉丝袜av| 一级a爱视频在线免费观看| 麻豆国产av国片精品| 国产日韩欧美亚洲二区| 黑人猛操日本美女一级片| 十八禁网站网址无遮挡| 中文字幕精品免费在线观看视频| 午夜福利在线观看吧| 婷婷成人精品国产| 欧美日韩福利视频一区二区| 999精品在线视频| 久久香蕉激情| 亚洲美女黄片视频| 熟女少妇亚洲综合色aaa.| 另类亚洲欧美激情| 亚洲色图综合在线观看| 中文字幕人妻丝袜制服| 国产精品亚洲一级av第二区| 国产极品粉嫩免费观看在线| 新久久久久国产一级毛片| 午夜久久久在线观看| 亚洲成av片中文字幕在线观看| 国产成人精品久久二区二区免费| 国产精品久久久久久人妻精品电影 | 国产精品自产拍在线观看55亚洲 | 一进一出好大好爽视频| 亚洲中文日韩欧美视频| 国产免费福利视频在线观看| 国产日韩欧美视频二区| 日本vs欧美在线观看视频| 久久亚洲精品不卡| 日日摸夜夜添夜夜添小说| 少妇裸体淫交视频免费看高清 | 美女扒开内裤让男人捅视频| 人妻 亚洲 视频| 黄色丝袜av网址大全| 五月开心婷婷网| 色老头精品视频在线观看| 99久久国产精品久久久| 国产高清激情床上av| 国精品久久久久久国模美| 亚洲七黄色美女视频| 欧美成人午夜精品| 精品国产一区二区久久| 在线观看免费视频网站a站| 啦啦啦中文免费视频观看日本| 黄频高清免费视频| 亚洲熟女毛片儿| 色婷婷av一区二区三区视频| 蜜桃在线观看..| 午夜日韩欧美国产| 视频在线观看一区二区三区| av线在线观看网站| 久久久久精品国产欧美久久久| 久久久国产精品麻豆| 黄色视频在线播放观看不卡| 午夜精品久久久久久毛片777| 亚洲第一av免费看| xxxhd国产人妻xxx| 国产福利在线免费观看视频| 亚洲精品国产区一区二| 大码成人一级视频| 国产一区二区三区综合在线观看| 欧美激情久久久久久爽电影 | 每晚都被弄得嗷嗷叫到高潮| 人妻久久中文字幕网| 亚洲久久久国产精品| 久久国产亚洲av麻豆专区| 亚洲黑人精品在线| 丝瓜视频免费看黄片| 久久久久国产一级毛片高清牌| 欧美黑人欧美精品刺激| 国产成+人综合+亚洲专区| 高清毛片免费观看视频网站 | 久久精品国产亚洲av高清一级| 国产日韩欧美视频二区| 日韩大片免费观看网站| 天堂动漫精品| 色老头精品视频在线观看| 日韩有码中文字幕| 国产av一区二区精品久久| 亚洲精品av麻豆狂野| 欧美精品人与动牲交sv欧美| 中文字幕制服av| 午夜成年电影在线免费观看| 一个人免费在线观看的高清视频| 国产精品亚洲一级av第二区| 日韩大片免费观看网站| 十八禁网站免费在线| 国产老妇伦熟女老妇高清| 在线观看免费午夜福利视频| 国产主播在线观看一区二区| 国产成人精品久久二区二区免费| 久久久久国内视频| 色在线成人网| 看免费av毛片| 电影成人av| 午夜免费成人在线视频| 国产精品一区二区免费欧美| 在线看a的网站| 999久久久国产精品视频| 成人18禁在线播放| 免费高清在线观看日韩| 丰满饥渴人妻一区二区三| 首页视频小说图片口味搜索| 国产一区二区三区综合在线观看| 黄色视频在线播放观看不卡| 老司机午夜十八禁免费视频| 人人妻人人爽人人添夜夜欢视频| 岛国在线观看网站| 日本一区二区免费在线视频| 蜜桃在线观看..| 中文字幕精品免费在线观看视频| 国产熟女午夜一区二区三区| 久久久精品免费免费高清| 国产成人精品久久二区二区免费| 桃红色精品国产亚洲av| 啦啦啦免费观看视频1| 国产免费视频播放在线视频| 一级,二级,三级黄色视频| 涩涩av久久男人的天堂| 亚洲精品久久午夜乱码| 日韩人妻精品一区2区三区| 国产男靠女视频免费网站| 亚洲欧美日韩高清在线视频 | 在线亚洲精品国产二区图片欧美| 男女高潮啪啪啪动态图| 国产精品自产拍在线观看55亚洲 | 国产在线免费精品| 又紧又爽又黄一区二区| 狂野欧美激情性xxxx| 欧美日本中文国产一区发布| 一级a爱视频在线免费观看| 成人国产av品久久久| 日日夜夜操网爽| 免费看十八禁软件| 亚洲欧美色中文字幕在线| 国产免费视频播放在线视频| 欧美亚洲 丝袜 人妻 在线| 啦啦啦免费观看视频1| 大片免费播放器 马上看| 两性夫妻黄色片| 国产高清激情床上av| av电影中文网址| 欧美日韩一级在线毛片| 亚洲成人国产一区在线观看| 成年人午夜在线观看视频| 一区二区av电影网| 中文字幕人妻丝袜一区二区|