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

    基于坐標(biāo)變換的雙饋風(fēng)力發(fā)電機(jī)組葉片質(zhì)量不平衡故障診斷

    2016-05-10 08:40:26繩曉玲萬書亭李永剛成立峰
    電工技術(shù)學(xué)報(bào) 2016年7期
    關(guān)鍵詞:雙饋發(fā)電機(jī)組風(fēng)力

    繩曉玲 萬書亭 李永剛 成立峰

    (1.華北電力大學(xué)電氣與電子工程學(xué)院 保定 071003

    2.華北電力大學(xué)能源動(dòng)力與機(jī)械工程學(xué)院 保定 071003)

    ?

    基于坐標(biāo)變換的雙饋風(fēng)力發(fā)電機(jī)組葉片質(zhì)量不平衡故障診斷

    繩曉玲1,2萬書亭2李永剛1成立峰2

    (1.華北電力大學(xué)電氣與電子工程學(xué)院保定071003

    2.華北電力大學(xué)能源動(dòng)力與機(jī)械工程學(xué)院保定071003)

    摘要首先分析了葉片質(zhì)量不平衡故障時(shí)雙饋風(fēng)力發(fā)電機(jī)組的電氣特性,提出了基于dq坐標(biāo)變換規(guī)則的定、轉(zhuǎn)子電流特性分析法;然后采用希爾伯特解調(diào)法對(duì)電流信號(hào)進(jìn)行處理,突出了電流的故障頻率;最后在Matlab/Simulink環(huán)境下進(jìn)行了仿真,并在雙饋風(fēng)力發(fā)電機(jī)組實(shí)驗(yàn)平臺(tái)上進(jìn)行了驗(yàn)證。結(jié)果表明,葉片質(zhì)量不平衡故障時(shí),發(fā)電機(jī)定、轉(zhuǎn)子電流中都包含故障頻率——葉片轉(zhuǎn)頻,該方法在葉片不平衡故障的識(shí)別中是有效的,具有一定的工程應(yīng)用價(jià)值。

    關(guān)鍵詞:雙饋風(fēng)力發(fā)電機(jī)組葉片質(zhì)量不平衡定、轉(zhuǎn)子電流坐標(biāo)變換

    0引言

    風(fēng)能作為一種清潔、綠色的可再生能源,已成為世界各國新能源發(fā)展的重要方向[1]。隨著風(fēng)力發(fā)電機(jī)組裝機(jī)容量的持續(xù)增加,如何降低機(jī)組維修成本和停機(jī)時(shí)間以及提高機(jī)組出力成為制約風(fēng)力發(fā)電發(fā)展的關(guān)鍵問題[2,3]。

    風(fēng)電機(jī)組大多運(yùn)行在環(huán)境比較惡劣的戶外或海上,風(fēng)沙的磨損或結(jié)冰沉積等會(huì)造成葉片質(zhì)量不平衡;另外疲勞應(yīng)力的作用會(huì)使葉片產(chǎn)生裂紋,灰塵、雜物、雨水等進(jìn)入裂紋也會(huì)引起葉片質(zhì)量的不平衡,并進(jìn)一步加劇葉片和其他部件的疲勞和振動(dòng),若繼續(xù)運(yùn)行會(huì)對(duì)機(jī)組產(chǎn)生非常大的危害,嚴(yán)重時(shí)葉片斷裂造成重大損失[4]。針對(duì)風(fēng)機(jī)葉片質(zhì)量不平衡問題,文獻(xiàn)[5]建立了風(fēng)機(jī)模型,利用有限元方法分析了葉片質(zhì)量不平衡引起的機(jī)組振動(dòng)特點(diǎn)。文獻(xiàn)[6,7]在模擬實(shí)驗(yàn)平臺(tái)上采集了機(jī)組振動(dòng)信號(hào)進(jìn)行頻譜分析,結(jié)果顯示質(zhì)量不平衡情況下機(jī)組在水平方向的振動(dòng)會(huì)增加。但振動(dòng)信號(hào)的采集需要安裝振動(dòng)傳感器,這首先增加了監(jiān)測(cè)成本,其次對(duì)于安裝在幾十米高空的機(jī)艙很難進(jìn)行實(shí)時(shí)監(jiān)測(cè)。文獻(xiàn)[8-10]則分析了葉片質(zhì)量不平衡與發(fā)電功率的關(guān)系。由于功率的計(jì)算需要三相電流和電壓,采集和計(jì)算都較麻煩,因此這些文獻(xiàn)只分析了功率有效值或單相功率的變化,故分析結(jié)果難免偏頗。文獻(xiàn)[11]分析了直驅(qū)風(fēng)力發(fā)電機(jī)組葉片質(zhì)量不平衡情況下發(fā)電機(jī)定子電流與軸轉(zhuǎn)頻的關(guān)系,利用定子電流信號(hào)間接分析轉(zhuǎn)頻變化特性來診斷葉片質(zhì)量不平衡。文獻(xiàn)[12]同樣也是針對(duì)直驅(qū)風(fēng)力發(fā)電機(jī)組,提出了利用定子電流信號(hào)模平方的方法來診斷葉片質(zhì)量不平衡。

    本文針對(duì)雙饋風(fēng)力發(fā)電機(jī)組葉片質(zhì)量不平衡故障時(shí)發(fā)電機(jī)電氣特性以及故障診斷等問題展開研究,找出了故障診斷依據(jù),實(shí)現(xiàn)了故障的早發(fā)現(xiàn)早診斷,避免了重大事故發(fā)生,提高了機(jī)組運(yùn)行效率。

    1葉片質(zhì)量不平衡機(jī)理分析

    葉片質(zhì)量不平衡故障可等效為在葉片上存在一個(gè)距離輪轂為R的不平衡虛擬質(zhì)量塊m[6,8]。等效不平衡質(zhì)量塊產(chǎn)生的重力矩會(huì)對(duì)風(fēng)機(jī)的輸出轉(zhuǎn)矩產(chǎn)生影響[13],風(fēng)機(jī)輸出的機(jī)械轉(zhuǎn)矩Tw可表示為

    Tw=Tw0+mgRsin(ωwt+φw)

    (1)

    式中,Tw0為氣動(dòng)轉(zhuǎn)矩;ωw為風(fēng)機(jī)葉片旋轉(zhuǎn)角速度;φw為葉片初始位置角度。

    圖1為風(fēng)機(jī)某一葉片質(zhì)量不平衡時(shí)的示意圖。等效質(zhì)量塊m隨著葉片以ωw的速度旋轉(zhuǎn),它產(chǎn)生的重力矩會(huì)造成風(fēng)機(jī)主軸(也稱葉片轉(zhuǎn)子軸)轉(zhuǎn)速的波動(dòng)。如圖1a所示,不平衡質(zhì)量塊由頂部向下旋轉(zhuǎn)到底部的過程中會(huì)給主軸加速;而當(dāng)該葉片由底部向上旋轉(zhuǎn)到頂部的過程中,等效質(zhì)量塊的存在會(huì)給主軸減速,如圖1b所示。這樣的周期變化頻率等于葉片轉(zhuǎn)頻[8]。另外文獻(xiàn)[14,15]中也論證了當(dāng)轉(zhuǎn)矩周期性波動(dòng)時(shí),軸轉(zhuǎn)速也會(huì)出現(xiàn)相同頻率的波動(dòng)。雙饋風(fēng)力發(fā)電機(jī)組通過增速齒輪箱將風(fēng)機(jī)和發(fā)電機(jī)連接起來,發(fā)電機(jī)轉(zhuǎn)子旋轉(zhuǎn)角速度等于Nωw(N為齒輪箱的增速比)。風(fēng)機(jī)主軸增速、減速的變化會(huì)反映到發(fā)電機(jī)轉(zhuǎn)子速度上,使轉(zhuǎn)子發(fā)生增速、減速的周期波動(dòng),波動(dòng)頻率等于葉片轉(zhuǎn)頻。

    根據(jù)雙饋風(fēng)力發(fā)電機(jī)組的運(yùn)動(dòng)方程有[16]

    (2)

    式中,Te為風(fēng)力發(fā)電機(jī)組的電磁轉(zhuǎn)矩;ωr為發(fā)電機(jī)轉(zhuǎn)子電角速度;ρ為極對(duì)數(shù);J為機(jī)組的轉(zhuǎn)動(dòng)慣量。

    將式(1)代入式(2),并令初相角φw=0,可求出轉(zhuǎn)子電角速度ωr為

    ωr=ωr0+ωrgcosωwt

    (3)

    式中,ωr0為氣動(dòng)轉(zhuǎn)矩對(duì)應(yīng)的轉(zhuǎn)子電角速度;ωrg為轉(zhuǎn)子速度周期波動(dòng)幅值,ωrg=ρmgR/Jωw。

    圖1 葉片質(zhì)量不平衡等效模型Fig.1 Equivalent model of blade mass imbalance

    2雙饋風(fēng)力發(fā)電機(jī)組特性分析

    2.1雙饋風(fēng)力發(fā)電機(jī)組工作原理

    雙饋風(fēng)力發(fā)電機(jī)組的結(jié)構(gòu)類似于繞線式感應(yīng)電機(jī),定、轉(zhuǎn)子均為三相對(duì)稱繞組,所不同的是轉(zhuǎn)子繞組與一個(gè)三相勵(lì)磁變換器相連,并在轉(zhuǎn)子繞組中施以轉(zhuǎn)差頻率為sω1的三相對(duì)稱電源,其中ω1為同步速。為了實(shí)現(xiàn)穩(wěn)定的機(jī)電能量轉(zhuǎn)換,定、轉(zhuǎn)子旋轉(zhuǎn)磁場必須均以同步速旋轉(zhuǎn),保持相對(duì)靜止,因此穩(wěn)定時(shí)各頻率應(yīng)滿足如下關(guān)系[16]

    ω1=ωr+sω1=ωr+ωz

    (4)

    式中,ωz為轉(zhuǎn)子電流角頻率;s為轉(zhuǎn)差率。

    當(dāng)風(fēng)力發(fā)電機(jī)組穩(wěn)態(tài)運(yùn)行時(shí),ω1應(yīng)等于電網(wǎng)電壓角頻率。在電網(wǎng)電壓恒定的情況下,ω1可認(rèn)為是一定值[16]。再結(jié)合式(3)和式(4)可得出質(zhì)量不平衡時(shí)轉(zhuǎn)子電流角頻率ωz以及其瞬時(shí)相位θr分別為

    ωz=ωz0-ωrgcos(ωwt)

    (5)

    (6)

    式中,ωz0=ω1-ωr0。

    2.2雙饋風(fēng)力發(fā)電機(jī)組轉(zhuǎn)、定子電流分析

    設(shè)ira為轉(zhuǎn)子a相電流,并設(shè)其初始相角為零,則ira可表示為(推導(dǎo)過程見附錄)

    (7)

    式中,Ir為轉(zhuǎn)子a相電流的幅值。

    由式(7)可知,在某一葉片存在質(zhì)量不平衡時(shí),轉(zhuǎn)子勵(lì)磁電流中除了頻率為fz0的正弦分量外,還會(huì)出現(xiàn)頻率為fz0±fw(ωz0=2πfz0,ωw=2πfw)的調(diào)制諧波分量。除此之外,調(diào)制電流的幅值為基波電流幅值的ωrg/(2ωw)倍,由式(3)可知,ωrg=ρmgR/(Jωw),對(duì)于大型風(fēng)電機(jī)組,不平衡重力矩mgR遠(yuǎn)小于機(jī)組轉(zhuǎn)動(dòng)慣量J,而ωw一般大于1,故調(diào)制電流的幅值相比基波電流幅值是很小的。

    為了進(jìn)一步分析雙饋發(fā)電機(jī)定子電流的特性,需首先將三相靜止坐標(biāo)系下的轉(zhuǎn)子電流ira、irb、irc變換到d-q旋轉(zhuǎn)坐標(biāo)系下。根據(jù)abc-dq變換規(guī)則[17],轉(zhuǎn)子d、q軸電流ird、irq計(jì)算過程見式(8)和式(9)。

    =ird0-irdw++irdw-

    (8)

    =irq0-irqw++irqw-

    (9)

    由式(8)和式(9)可知,轉(zhuǎn)子勵(lì)磁電流可以分解到3個(gè)d-q坐標(biāo)系中,3個(gè)d-q坐標(biāo)系的旋轉(zhuǎn)速度分別是ωz0、ωz0+ωw和ωz0-ωw。其中以ωz0速度旋轉(zhuǎn)的d、q軸電流分別為ird0、irq0;以ωz0+ωw速度旋轉(zhuǎn)的d、q軸電流分別為-irdw+、-irqw+;以ωz0-ωw速度旋轉(zhuǎn)的d、q軸電流分別為irdw-、irqw-。

    根據(jù)雙饋風(fēng)力發(fā)電機(jī)組數(shù)學(xué)模型以及定子磁鏈定向原則有[16]

    (10)

    式中,ψsd和ψsq分別為定子磁鏈的d、q軸分量;ψs為定子磁鏈的幅值;isd和isq分別為定子電流的d、q軸分量;Ls為定子自感;Lm為定、轉(zhuǎn)子間互感。

    將式(8)和式(9)帶入式(10)可得

    (11)

    也即定子的d、q軸電流可表達(dá)為

    (12)

    由式(12)可知,定子電流也能分解到3個(gè)旋轉(zhuǎn)的d、q軸上,設(shè)3個(gè)d、q軸的轉(zhuǎn)速分別為ω0、ω+和ω-,為了使定子磁場與轉(zhuǎn)子磁場保持相對(duì)靜止,應(yīng)有ω0=ωz0+ωr0=ω1、ω+=(ωz0+ωw)+ωr0=ω1+ωw、ω-=(ωz0-ωw)+ωr0=ω1-ωw。故定子電流分解后的3個(gè)d、q軸轉(zhuǎn)速分別為ω1、ω1+ωw和ω1-ωw。以ω1速度旋轉(zhuǎn)的d、q軸電流分別為isd0、isq0;以ω1+ωw速度旋轉(zhuǎn)的d、q軸電流分別為-isdw+、-isqw+;以ω1-ωw速度旋轉(zhuǎn)的d、q軸電流分別為isdw-和isqw-。

    根據(jù)dq-abc坐標(biāo)變換規(guī)則以及式(8)和式(9)的推導(dǎo)過程,可得出發(fā)電機(jī)定子a相電流的表達(dá)式為

    isa=Iscos(ω1t+φ0)-Is1cos(ω1t+ωwt+φ1)+

    Is2cos(ω1t-ωwt+φ2)

    (13)

    式中

    由式(13)可知,當(dāng)某一葉片存在質(zhì)量不平衡時(shí),定子電流中除了頻率為f1的基波外,還會(huì)出現(xiàn)頻率為f1±fw(ω1=2πf1,ωw=2πfw)的調(diào)制諧波分量。同時(shí)在式(7)中也分析出來,轉(zhuǎn)子勵(lì)磁電流中也存在頻率為fz0±fw的調(diào)制諧波分量,其中fw通常被稱作故障頻率或調(diào)制頻率。由此便得出了利用定、轉(zhuǎn)子電流診斷風(fēng)力發(fā)電機(jī)組葉片質(zhì)量不平衡故障的依據(jù):定、轉(zhuǎn)子電流中除了各自的基頻外,還存在調(diào)制頻率fw,該頻率為風(fēng)機(jī)葉片轉(zhuǎn)頻。

    通過前述分析可知,定、轉(zhuǎn)子電流的諧波幅值與基波電流相比非常小,因此直接分析電流的FFT(快速傅里葉變換)很難看出調(diào)制頻率。而信號(hào)處理中的希爾伯特解調(diào)方法能夠從信號(hào)中解調(diào)出低頻的調(diào)制分量,因此在本文中通過對(duì)電流信號(hào)進(jìn)行希爾伯特解調(diào)分析來提取故障頻率。

    3希爾伯特解調(diào)法

    希爾伯特解調(diào)法通過希爾伯特變換的方法從信號(hào)中提取調(diào)制信號(hào),分析調(diào)制函數(shù)的變化來提取故障特征,在機(jī)械振動(dòng)信號(hào)及其故障診斷方面得到廣泛應(yīng)用,但在電氣信號(hào)分析方面應(yīng)用還較少。文獻(xiàn)[18]曾將該方法應(yīng)用于電機(jī)故障診斷。

    若一連續(xù)的時(shí)間信號(hào)為x(t),其希爾伯特變換可看成是x(t)通過一濾波器的輸出[18]

    (14)

    x(t)與其希爾伯特變換構(gòu)成x(t)的解析信號(hào)Zt

    (15)

    由式(15)可得

    (16)

    式中,A(t)為幅值包絡(luò);θ(t)為相位包絡(luò);μ(t)為頻率包絡(luò)。

    A(t)、θ(t)和μ(t)分別包含x(t)的幅值、相位和頻率調(diào)制信息,對(duì)這些包絡(luò)進(jìn)行頻譜分析便可得到調(diào)制頻率,實(shí)現(xiàn)解調(diào)。

    接下來以轉(zhuǎn)子電流為例,按照上述步驟分析希爾伯特頻率包絡(luò)解調(diào)的過程。轉(zhuǎn)子電流可表達(dá)為

    (17)

    求出ira的希爾伯特變換后,可得其解析信號(hào)為

    (18)

    式中,相位包絡(luò)θ(t)可表示為

    (19)

    對(duì)θ(t)求導(dǎo)可得

    (20)

    將式(20)中的μ(t)去直流并進(jìn)行FFT變換便可得到調(diào)制頻率ωw,這就是希爾伯特變換的解調(diào)功能。同理也可求出定子電流的調(diào)制頻率ωw。

    基于希爾伯特變換的解調(diào)法能夠分離出原信號(hào)中的調(diào)制信息進(jìn)行解調(diào),消除了原始信號(hào)的基頻,只剩下調(diào)制頻率,得到的結(jié)果比較清晰直觀。

    4仿真分析

    4.11.5 MW機(jī)組仿真

    為了驗(yàn)證前文所述理論分析的正確性,利用Matlab/Simulink模塊,搭建了一個(gè)1.5 MW的雙饋風(fēng)力發(fā)電機(jī)組平臺(tái)。平臺(tái)模型結(jié)構(gòu)如圖2所示。仿真模型參數(shù)見表1,雙饋風(fēng)力發(fā)電機(jī)組轉(zhuǎn)子側(cè)采用定子磁鏈定向矢量控制[16]。

    圖2 仿真平臺(tái)Fig.2 Simulation platform

    參數(shù)數(shù)值參數(shù)數(shù)值額定功率/MW1.5額定電壓/V690額定風(fēng)速/(m/s)12定子電阻(pu)0.0059葉片額定轉(zhuǎn)速/(r·min-1)13.2定子漏感(pu)0.1425最佳風(fēng)能利用系數(shù)0.48轉(zhuǎn)子電阻(pu)0.0042最佳葉尖速比5.8轉(zhuǎn)子漏感(pu)0.13齒輪箱增速比90極對(duì)數(shù)3慣性常數(shù)5s電網(wǎng)頻率/Hz50

    仿真前首先定義了不平衡度系數(shù)的概念,不平衡度系數(shù)b可表示為

    (21)

    式中,Twg為不平衡重力矩。

    平臺(tái)搭建好之后,首先進(jìn)行了正常工況的仿真。然后在保證機(jī)組穩(wěn)態(tài)運(yùn)行的條件下,故障程度由低到高逐漸增加,分別進(jìn)行了1%、3%、5%、8%和10%五種不平衡工況的仿真。為了便于對(duì)比,這五種工況下風(fēng)機(jī)的運(yùn)行參數(shù)相同,均為恒定風(fēng)速10 m/s,葉片轉(zhuǎn)速12.1 r/min。

    圖3為正常以及10%不平衡工況下葉片轉(zhuǎn)速的對(duì)比圖。從時(shí)域圖3a可看出,在不平衡工況下葉片轉(zhuǎn)速有明顯的周期波動(dòng);經(jīng)FFT變換后其頻率為0.2 Hz,也即葉片轉(zhuǎn)頻(12.1/60≈0.2),如圖3b所示。

    圖3 正常及10%不平衡工況下葉片轉(zhuǎn)速Fig.3 Blade rotate-speed of normal and 10% imbalance

    圖4為正常以及10%不平衡工況下定子電流仿真分析的結(jié)果。由圖4a可知,正常以及不平衡工況下定子電流在時(shí)域圖形上基本重合。在圖4b中,正常及不平衡定子電流在50 Hz處幾乎完全重合,將49.5~50.5 Hz頻段的縱坐標(biāo)放大幾百倍后,才能看到二者的區(qū)別。圖4c是對(duì)定子電流進(jìn)行解調(diào)的結(jié)果:不平衡故障下的定子電流在葉片轉(zhuǎn)頻(即調(diào)制頻率)0.2 Hz處有明顯的峰值,而正常情況下的定子電流則沒有。由此可見,希爾伯特包絡(luò)解調(diào)可凸顯故障頻率,結(jié)果清晰直觀。

    圖4 正常及10%不平衡工況下定子電流Fig.4 Srator current of normal and 10% imbalance

    圖5為正常以及10%不平衡工況下轉(zhuǎn)子電流仿真分析的結(jié)果。由圖5a和圖5b可知,正常以及不平衡時(shí)轉(zhuǎn)子電流的基頻均約為5 Hz,但與圖4a中定子電流時(shí)域波形不同,在圖5a中,正常及故障下的轉(zhuǎn)子電流有明顯相位差,這與前面式(6)的分析一致。由于轉(zhuǎn)子電流基頻較小,因此在故障前后轉(zhuǎn)子電流的相位變化明顯。而定子電流基頻相對(duì)較大,因此在定子電流的時(shí)域波形圖上看不出明顯差別。不平衡轉(zhuǎn)子電流中解調(diào)出葉片轉(zhuǎn)頻0.2 Hz,如圖5c所示。

    除此之外,對(duì)于其他幾種不平衡工況下的定、轉(zhuǎn)子電流也進(jìn)行了相同的分析,結(jié)果與上述分析一致。圖6為在1%、3%、5%、8%和10%五種不平衡工況下解調(diào)出的葉片轉(zhuǎn)頻0.2 Hz處所對(duì)應(yīng)的幅值。從圖6

    中可看出,隨著不平衡程度的增加,故障頻率所對(duì)應(yīng)的幅值也隨著增加。

    圖5 正常及10%不平衡工況下轉(zhuǎn)子電流Fig.5 Rotor current of normal and 10% imbalance

    圖6 各不平衡工況對(duì)應(yīng)的故障頻率幅值Fig.6 Amplitude of fault frequency for different degree of imbalance

    4.2實(shí)驗(yàn)機(jī)組仿真

    目前商業(yè)化的風(fēng)電機(jī)組都越來越大型化,兆瓦級(jí)風(fēng)機(jī)是當(dāng)前主流,但對(duì)于實(shí)驗(yàn)室類型的風(fēng)機(jī)來說,一般都是小型的。筆者所在課題組有一套11 kW雙饋風(fēng)力發(fā)電機(jī)組,它是某兆瓦級(jí)風(fēng)機(jī)的縮比模型。為了將仿真和實(shí)驗(yàn)結(jié)果進(jìn)行對(duì)比,利用Matlab搭建了一個(gè)11 kW實(shí)驗(yàn)機(jī)組的仿真模型。該仿真平臺(tái)參數(shù)與實(shí)驗(yàn)室11 kW機(jī)組參數(shù)一致,具體見表2。

    表2 11 kW DFIG仿真參數(shù)

    平臺(tái)搭建好之后,進(jìn)行了正常工況和2.4%不平衡工況下的仿真,其中2.4%不平衡度工況對(duì)應(yīng)1 kg不平衡質(zhì)量塊的實(shí)驗(yàn)(見后文實(shí)驗(yàn)部分)。仿真風(fēng)況為恒定風(fēng)速6.4 m/s,對(duì)應(yīng)葉片轉(zhuǎn)速75 r/min。仿真結(jié)果如圖7和圖8所示。

    圖7 正常及不平衡工況下定子電流Fig.7 Srator current of normal and imbalance

    圖8 正常及不平衡工況下轉(zhuǎn)子電流Fig.8 Rotor current of normal and imbalance

    由圖7a可知,正常及不平衡情況下的定子電流在時(shí)域波形圖上幾乎完全重合;圖7b是將48~52 Hz頻段的縱坐標(biāo)放大了近103倍后才看到二者的區(qū)別。由圖7c可知,從不平衡的定子電流中解調(diào)出故障頻率1.25 Hz(75/60=1.25),而正常情況下則沒有,區(qū)別較明顯。

    與1.5 MW仿真類似,在圖8a中,正常及不平衡的轉(zhuǎn)子電流在時(shí)域圖形上有相位差,但在圖8b FFT頻譜圖中二者區(qū)別不大。在圖8c中,不平衡轉(zhuǎn)子電流中在1.25 Hz處有明顯峰值。

    5實(shí)驗(yàn)分析

    為了進(jìn)一步驗(yàn)證本文所提理論方法的正確性,在一個(gè)11 kW雙饋風(fēng)力發(fā)電平臺(tái)上進(jìn)行了實(shí)驗(yàn)。實(shí)驗(yàn)平臺(tái)如圖9所示。

    圖9 11 kW雙饋風(fēng)力發(fā)電實(shí)驗(yàn)平臺(tái)Fig.9 Test platform of 11 kW doubly-fed wind turbine

    該實(shí)驗(yàn)平臺(tái)的工作原理是:三葉片風(fēng)力機(jī)由調(diào)速電機(jī)(安裝在風(fēng)力機(jī)內(nèi)部)拖動(dòng),模擬風(fēng)速變化時(shí)葉片的旋轉(zhuǎn),從而帶動(dòng)主軸一起旋轉(zhuǎn)。主軸通過增速齒輪箱與發(fā)電機(jī)相連。調(diào)速電機(jī)控制、變速恒頻控制等都集成在主控柜中,如圖9b所示。

    發(fā)電機(jī)為繞線式異步發(fā)電機(jī),參數(shù)見表2。為了模擬質(zhì)量不平衡故障,在葉片上附加了一個(gè)重約1 kg(不平衡度b=2.4%)的質(zhì)量塊。為了防止對(duì)葉片表面造成損傷以及出于安全的考慮,將質(zhì)量塊用強(qiáng)力AB膠和若干條膠帶固定在葉片內(nèi)腔壁較平坦的地方,如圖9a所示位置的內(nèi)側(cè)。

    采集儀選用Yokogawa公司生產(chǎn)的DL850型示波記錄儀,并采用型號(hào)為96001的鉗式電流探頭采集定子電流,而轉(zhuǎn)子電流則選用HIA-C01型電流傳感器采集,該傳感器基于霍爾效應(yīng),適于低頻信號(hào)。

    實(shí)驗(yàn)參數(shù):風(fēng)力機(jī)設(shè)定為恒定風(fēng)模式,葉片轉(zhuǎn)速為75 r/min,轉(zhuǎn)差率s為0.11,發(fā)電機(jī)并網(wǎng)運(yùn)行。數(shù)據(jù)采集頻率1 024 Hz,采樣時(shí)間3 min,分別采集了機(jī)組正常工況和1 kg不平衡工況下的定、轉(zhuǎn)子電流,并進(jìn)行了分析。

    圖10為正常工況下采集的定、轉(zhuǎn)子電流,圖11為不平衡工況下的定、轉(zhuǎn)子電流。首先對(duì)比圖10與圖11中的定子電流ia,基本看不出二者區(qū)別,這與仿真圖7a是一致的。對(duì)這兩種工況下的定子電流進(jìn)行FFT分析,電流基頻均為50 Hz,與圖7b類似。再對(duì)比圖10與圖11中的轉(zhuǎn)子電流,正常及不平衡工況下轉(zhuǎn)子電流仍有相位差。在FFT頻譜圖中兩者的基頻均約為5.5 Hz,與圖8b類似。

    圖10 正常定、轉(zhuǎn)子電流Fig.10 Stator and rotor current of normal

    圖11 不平衡定、轉(zhuǎn)子電流Fig.11 Stator and rotor current of imbalance

    圖12為對(duì)正常和不平衡情況下采集的定、轉(zhuǎn)子電流進(jìn)行希爾伯特解調(diào)分析。從圖12a中可看出,在不平衡的定子電流中解調(diào)出故障頻率1.25 Hz;而在圖12b中,則從不平衡情況下的轉(zhuǎn)子電流中解調(diào)出故障頻率1.25 Hz。

    圖12 定、轉(zhuǎn)子電流解調(diào)Fig.12 Demodulation of stator and rotor current

    對(duì)比11 kW機(jī)組仿真和實(shí)驗(yàn)結(jié)果可發(fā)現(xiàn),定、轉(zhuǎn)子電流的時(shí)域波形、FFT頻域以及希爾伯特解調(diào)的結(jié)果都基本一致,不同之處主要是兩種情況下電流的幅值有差別,這是因?yàn)榉抡娣治龅碾娏魇窍嚯娏?,而?shí)驗(yàn)時(shí)電流傳感器采集的是線電流,這對(duì)分析結(jié)果影響不大。通過上述仿真和實(shí)驗(yàn)結(jié)果的分析驗(yàn)證了本文所提理論分析的正確性。

    6結(jié)論

    本文針對(duì)雙饋風(fēng)力發(fā)電機(jī)組葉片質(zhì)量不平衡故障,從雙饋發(fā)電機(jī)的數(shù)學(xué)模型出發(fā),提出了基于dq坐標(biāo)變換的定、轉(zhuǎn)子電流特性分析法來診斷葉片不平衡故障,并采用希爾伯特解調(diào)法從電流信號(hào)中解調(diào)出故障頻率,結(jié)果比較清晰直觀。通過理論分析、仿真和實(shí)驗(yàn)分析得出以下結(jié)論:

    1)葉片質(zhì)量不平衡故障會(huì)導(dǎo)致雙饋發(fā)電機(jī)定、轉(zhuǎn)子電流中都包含以葉片一倍轉(zhuǎn)頻為調(diào)制頻率的故障頻率,并隨著不平衡度系數(shù)的增加,故障頻率對(duì)應(yīng)幅值也增加,在一定程度上能夠表征故障的嚴(yán)重程度。

    2)除此之外,不平衡故障情況下轉(zhuǎn)子電流與正常情況下相比其相位發(fā)生變化,這是轉(zhuǎn)子電流頻率存在波動(dòng)所導(dǎo)致的。

    基于雙饋發(fā)電機(jī)定、轉(zhuǎn)子電流特性來監(jiān)測(cè)風(fēng)機(jī)葉片的運(yùn)行狀態(tài),及時(shí)消除引起不平衡的因素,防微杜漸,對(duì)于整個(gè)機(jī)組的安全、穩(wěn)定運(yùn)行具有重要的實(shí)際意義。

    附錄

    式(7)的推導(dǎo)過程

    ira=Ircosθr

    式(8)和式(9)中sinθr的推導(dǎo)過程

    注:在σ?1時(shí),cos(σ)≈1,sin(σ)≈σ。

    參考文獻(xiàn)

    [1]GWEC.Global wind report annual market update 2011[R].Bonn:Global Wind Energy Council,2011.

    [2]Hameed Z,Hong Y S,Cho Y M,et al.Condition monitoring and fault detetion of wind turbines and related algorithms:a review[J].Renewable and Sustainable Energy Reviews,2009,13(1):1-39.

    [3]杭俊,張建忠,程明,等.風(fēng)力發(fā)電系統(tǒng)狀態(tài)監(jiān)測(cè)與故障診斷技術(shù)綜述[J].電工技術(shù)學(xué)報(bào),2013,28(4):261-271.

    Hang Jun,Zhang Jianzhong,Cheng Ming,et al.An overview of condition monitoring and fault diagnostic for wind energy conversion system[J].Transactions of China Electrotechnical Society,2013,28(4):261-271.

    [4]蔡新,潘盼,朱杰,等.風(fēng)力發(fā)電機(jī)葉片[M].北京:中國水利水電出版社,2014.

    [5]Ramlau R,Niebsch J.Imbalance estimation without test masses for wind turbines[J].Journal of Solar Energy Engineering,2009,131(1):0110101-0110107.

    [6]Jiang Dongxiang,Huang Qian,Hong Liangyou.Theoretical and experimental study on wind wheel imbalance for a wind turbine[C]//World Non-Grid-Connected Wind Power and Energy Conference,Nanjing,China,2009:351-355.

    [7]Zhao Minghao,Jiang Dongxiang,Li Shaohua.Research on fault mechanism of icing of wind turbine blades[C]//World Non-Grid-Connected Wind Power and Energy Conference,Nanjing,China,2009:356-359.

    [8]Xiang Gong,Wei Qiao.Simulation Investigation of wind turbine imbalance faults[C]//International Conference on Power System Technology,Hangzhou,China,2010:1-7.

    [9]Caselitz P,Giebhardt J.Rotor condition monitoring for improved operational safety of offshore wind energy converters[J].Journal of Solar Energy Engineering,2005,127(2):253-261.

    [10]楊濤,任永,劉霞,等.風(fēng)力機(jī)葉輪質(zhì)量不平衡故障建模及仿真研究[J].機(jī)械工程學(xué)報(bào),2012,48(6):130-135.

    Yang Tao,Ren Yong,Liu Xia,et al.Research on the modeling and simulation of wind turbine rotor Imbalance fault[J].Journal of Mechanical Engineering,2012,48(6):130-135.

    [11]Xiang Gong,Wei Qiao.Imbalance Fault detection of direct-drive wind turbines using generator current signals[J].IEEE Transactions on Energy Conversion,2012,27(2):468-476.

    [12]杭俊,張建忠,程明,等.直驅(qū)永磁同步風(fēng)電機(jī)組葉輪不平衡和繞組不對(duì)稱的故障診斷[J].中國電機(jī)工程學(xué)報(bào),2014,34((9):1384-1391.

    Hang Jun,Zhang Jianzhong,Cheng Ming,et al.Fault diagnosis of blade imbalance and winding asymmetry ofdirect-driven wind turbine with permanent magnet synchronous generator[J].Proceedings of the CSEE,2014,34(9):1384-1391.

    [13]Niebsch J,Ramlau R,Nguyen T T.Mass and aerodynamic imbalance estimates of wind turbines[J].Energies,2010,3(4):696-710.

    [14]Martin B,Pierre G,Bertrand R,et al.Models for bearing damage detection induction motors using stator current monitoring[J].IEEE Transactions on Industrial Electronics,2008,55(4):1813-1822.

    [15]Yang Wenxian,Tavner P J,Wilkinson M.Wind turbine condition monitoring and fault diagnosis using both mechanical and electrical signatures[C]//International Conference on Advanced Intelligent Mechatronics,Xi’an,2008:1296-1301.

    [16]賀益康,胡家兵,徐烈.并網(wǎng)雙饋異步風(fēng)力發(fā)電機(jī)運(yùn)行控制[M].北京:中國電力出版社,2011.

    [17]張波,顏湘武,王樹岐.基于dq變換的三相電壓暫降生成方法[J].電工技術(shù)學(xué)報(bào),2011,26(10):55-61.Zhang Bo,Yan Xiangwu,Wang Shuqi.Method of three-phase voltage sag forming based on dq transformation[J].Transactions of China Electrotechnical Society,2011,26(10):55-61.

    [18]趙妍,李志民,李天云.一種基于譜峭度的異步電機(jī)故障診斷方法[J].電工技術(shù)學(xué)報(bào),2014,29(5):189-196.

    Zhao Yan,Li Zhimin,Li Tianyun.A method for fault diagnosis of induction motor based on spectral kurtosis[J].Transactions of China Electrotechnical Society,2014,29(5):189-196.

    繩曉玲女,1983年生,博士研究生,講師,研究方向?yàn)轱L(fēng)力發(fā)電機(jī)組的狀態(tài)監(jiān)測(cè)和故障診斷。

    E-mail:sssxlyc@163.com(通信作者)

    萬書亭男,1970年生,教授,博士生導(dǎo)師,研究方向?yàn)榇笮突痣姍C(jī)組、風(fēng)力發(fā)電機(jī)組在線監(jiān)測(cè)與故障診斷。

    E-mail:13582996591@139.com

    Fault Diagnosis for Blade Mass Imbalance of Wind Turbines with DFIG Based on Coordinate Transformation

    ShengXiaoling1,2WanShuting2LiYonggang1ChengLifeng2

    (1.School of Electrical and Electronic EngineeringNorth China Electric Power University Baoding071003China 2.School of Energy Power and Mechanical EngineeringNorth China Electric Power University Baoding071003China)

    AbstractFirstly, the electrical characteristics of the doubly-fed induction generator (DFIG) under the blade imbalance fault are analyzed.Furthermore, the stator and rotor currents can be obtained according to the dq coordinate transformation.Then the Hilbert demodulation method is used to process the current signal so that the fault frequency can be highlighted.Finally, the simulation work under Matlab/Simulink and the experiment on the DFIG testing platform are taken for verification.It is shown that, the faulty stator and rotor currents contain the same harmonic component as that of the blade’s rotating frequency.The proposed method is effective to recognize the blade mass imbalance fault and therefore has an engineering application potential.

    Keywords:Doubly-fed induction generator, blade, mass imbalance, stator and rotor current, coordinate transformation

    作者簡介

    中圖分類號(hào):TM315

    收稿日期2015-03-29改稿日期2015-07-02

    國家自然科學(xué)基金(51177046)和中央高?;究蒲袠I(yè)務(wù)費(fèi)專項(xiàng)資金(12MS101,2014XS82)資助。

    猜你喜歡
    雙饋發(fā)電機(jī)組風(fēng)力
    走進(jìn)風(fēng)力發(fā)電
    煤氣發(fā)電機(jī)組DEH控制系統(tǒng)的優(yōu)化
    山東冶金(2022年4期)2022-09-14 08:59:30
    雙饋式可變速抽水蓄能機(jī)組運(yùn)行控制
    基于Motor CAD的雙饋發(fā)電機(jī)溫升速算方法研究
    大型風(fēng)力發(fā)電設(shè)備潤滑概要
    基于PLC控制柴油發(fā)電機(jī)組3D 模型
    受平均斜率控制的Crowbar雙饋異步電機(jī)低電壓穿越
    基于PI控制的雙饋風(fēng)電場的無功電壓建模與仿真
    八鋼歐冶爐TRT發(fā)電機(jī)組成功并網(wǎng)發(fā)電
    新疆鋼鐵(2015年3期)2015-02-20 14:13:56
    國內(nèi)首臺(tái)65MW超高壓中間再熱發(fā)電機(jī)組并網(wǎng)發(fā)電
    欧美一级毛片孕妇| 亚洲精品久久成人aⅴ小说| 国产成人av激情在线播放| 国产极品粉嫩免费观看在线| 视频区欧美日本亚洲| 午夜精品国产一区二区电影| 欧美日韩成人在线一区二区| 搡老岳熟女国产| 看免费av毛片| 亚洲综合色网址| 国产精品亚洲av一区麻豆| 日本黄色日本黄色录像| av免费在线观看网站| www.精华液| 国产色视频综合| 亚洲精品美女久久久久99蜜臀| 侵犯人妻中文字幕一二三四区| 国产激情久久老熟女| 老鸭窝网址在线观看| 成年动漫av网址| 妹子高潮喷水视频| 国产精品电影一区二区三区 | 午夜两性在线视频| 色尼玛亚洲综合影院| 色94色欧美一区二区| 免费黄频网站在线观看国产| 亚洲午夜理论影院| 免费观看a级毛片全部| 80岁老熟妇乱子伦牲交| 女性被躁到高潮视频| 国产黄频视频在线观看| 欧美激情 高清一区二区三区| 亚洲成人免费电影在线观看| 欧美日韩亚洲国产一区二区在线观看 | 日韩成人在线观看一区二区三区| 欧美乱码精品一区二区三区| 制服诱惑二区| 国产精品.久久久| 亚洲成人免费电影在线观看| 丰满人妻熟妇乱又伦精品不卡| 国产精品 国内视频| 少妇精品久久久久久久| 色播在线永久视频| 日韩中文字幕视频在线看片| 亚洲精华国产精华精| 热99久久久久精品小说推荐| 久久久精品区二区三区| 日日摸夜夜添夜夜添小说| 亚洲国产欧美在线一区| 久久久久国内视频| 国产极品粉嫩免费观看在线| 亚洲免费av在线视频| 国产黄频视频在线观看| 免费高清在线观看日韩| 一个人免费看片子| 亚洲欧美激情在线| 免费一级毛片在线播放高清视频 | 99国产精品一区二区蜜桃av | 纯流量卡能插随身wifi吗| 美女福利国产在线| 日本撒尿小便嘘嘘汇集6| a级毛片在线看网站| 亚洲人成77777在线视频| 两个人免费观看高清视频| 国产欧美日韩一区二区三区在线| 美女国产高潮福利片在线看| 久久青草综合色| 亚洲成国产人片在线观看| 日本a在线网址| 捣出白浆h1v1| 乱人伦中国视频| 搡老熟女国产l中国老女人| 别揉我奶头~嗯~啊~动态视频| 色94色欧美一区二区| avwww免费| 热re99久久精品国产66热6| 国产99久久九九免费精品| 成人国产av品久久久| 1024视频免费在线观看| 国产欧美日韩精品亚洲av| 中文字幕制服av| 国产精品久久久久成人av| 免费女性裸体啪啪无遮挡网站| 欧美亚洲日本最大视频资源| 欧美日韩国产mv在线观看视频| 精品一区二区三区四区五区乱码| 电影成人av| 每晚都被弄得嗷嗷叫到高潮| 久久婷婷成人综合色麻豆| 丰满少妇做爰视频| 女性被躁到高潮视频| 国内毛片毛片毛片毛片毛片| 精品国产一区二区三区四区第35| 在线观看66精品国产| 黄色怎么调成土黄色| 2018国产大陆天天弄谢| 午夜精品久久久久久毛片777| 日韩三级视频一区二区三区| 精品一区二区三卡| 国产精品久久久久久精品古装| 99国产精品一区二区三区| xxxhd国产人妻xxx| 午夜久久久在线观看| 十分钟在线观看高清视频www| 久久精品亚洲av国产电影网| 免费观看人在逋| www.999成人在线观看| 男女午夜视频在线观看| 日韩中文字幕欧美一区二区| 午夜福利视频精品| 成人18禁在线播放| 欧美黄色淫秽网站| 午夜老司机福利片| 亚洲精品自拍成人| 免费观看人在逋| 男人操女人黄网站| 大片电影免费在线观看免费| 亚洲第一欧美日韩一区二区三区 | 日本撒尿小便嘘嘘汇集6| 在线观看免费日韩欧美大片| 母亲3免费完整高清在线观看| 天堂中文最新版在线下载| 精品久久久久久久毛片微露脸| 亚洲成av片中文字幕在线观看| 另类亚洲欧美激情| 在线观看一区二区三区激情| 自线自在国产av| 欧美 日韩 精品 国产| 女同久久另类99精品国产91| 欧美国产精品va在线观看不卡| 欧美乱码精品一区二区三区| 精品少妇内射三级| 美女国产高潮福利片在线看| 成人精品一区二区免费| 久久午夜综合久久蜜桃| 久久久久视频综合| 一区在线观看完整版| 中亚洲国语对白在线视频| 久久99一区二区三区| 精品国产一区二区三区久久久樱花| 欧美黑人欧美精品刺激| 欧美激情 高清一区二区三区| 午夜福利一区二区在线看| 天天躁狠狠躁夜夜躁狠狠躁| 国产在线一区二区三区精| 国产视频一区二区在线看| 757午夜福利合集在线观看| 热re99久久国产66热| 精品国内亚洲2022精品成人 | 国产三级黄色录像| 欧美日本中文国产一区发布| 一二三四社区在线视频社区8| 天堂动漫精品| 99国产精品免费福利视频| 日本精品一区二区三区蜜桃| 深夜精品福利| 国产精品二区激情视频| 人妻一区二区av| 宅男免费午夜| 欧美老熟妇乱子伦牲交| 天堂动漫精品| 日本精品一区二区三区蜜桃| 国产欧美日韩精品亚洲av| 啦啦啦 在线观看视频| 国产三级黄色录像| 露出奶头的视频| 99国产综合亚洲精品| 国产欧美亚洲国产| 亚洲精品美女久久久久99蜜臀| 性高湖久久久久久久久免费观看| 久久av网站| 久久人人97超碰香蕉20202| 欧美日韩国产mv在线观看视频| 美女扒开内裤让男人捅视频| 国产精品美女特级片免费视频播放器 | 午夜福利免费观看在线| 国产欧美日韩一区二区三| 国产日韩欧美在线精品| 91字幕亚洲| 午夜视频精品福利| 飞空精品影院首页| 国产精品亚洲一级av第二区| 激情视频va一区二区三区| 精品少妇一区二区三区视频日本电影| 岛国在线观看网站| 国产成+人综合+亚洲专区| 亚洲中文日韩欧美视频| 亚洲国产中文字幕在线视频| 精品欧美一区二区三区在线| 欧美另类亚洲清纯唯美| 正在播放国产对白刺激| 国产精品电影一区二区三区 | 成人国产av品久久久| 精品久久久精品久久久| 亚洲中文日韩欧美视频| 男男h啪啪无遮挡| 亚洲中文日韩欧美视频| 一本综合久久免费| 岛国在线观看网站| 国产精品免费视频内射| 欧美成人午夜精品| 欧美日韩亚洲高清精品| 日韩制服丝袜自拍偷拍| av在线播放免费不卡| 国产av一区二区精品久久| 人妻一区二区av| 捣出白浆h1v1| 国产精品自产拍在线观看55亚洲 | 国产亚洲精品第一综合不卡| 97在线人人人人妻| 久久精品国产a三级三级三级| 国产亚洲精品第一综合不卡| av在线播放免费不卡| 1024香蕉在线观看| 大陆偷拍与自拍| 视频区欧美日本亚洲| 成在线人永久免费视频| 99国产综合亚洲精品| 精品国产乱码久久久久久男人| xxxhd国产人妻xxx| 国产精品 国内视频| 午夜福利在线免费观看网站| 久久久水蜜桃国产精品网| 老汉色av国产亚洲站长工具| 国产一区二区三区综合在线观看| 高清视频免费观看一区二区| 亚洲av美国av| 国产成人精品在线电影| 一夜夜www| 久久精品91无色码中文字幕| 一区在线观看完整版| 久久久精品国产亚洲av高清涩受| 久久久精品国产亚洲av高清涩受| 精品久久蜜臀av无| 国产高清videossex| 日日夜夜操网爽| 日日摸夜夜添夜夜添小说| 啦啦啦在线免费观看视频4| 热re99久久国产66热| 精品福利观看| 亚洲全国av大片| 国产又爽黄色视频| 亚洲精品一二三| 精品久久久久久久毛片微露脸| 欧美日韩福利视频一区二区| 啦啦啦视频在线资源免费观看| 侵犯人妻中文字幕一二三四区| 亚洲精品国产区一区二| 欧美在线黄色| 国产一区二区激情短视频| 老鸭窝网址在线观看| 国产精品偷伦视频观看了| av片东京热男人的天堂| 搡老乐熟女国产| a级毛片黄视频| 热99久久久久精品小说推荐| 狠狠婷婷综合久久久久久88av| 国产区一区二久久| 99国产精品一区二区三区| 成人特级黄色片久久久久久久 | 免费黄频网站在线观看国产| a级片在线免费高清观看视频| 国产成人精品在线电影| 日本vs欧美在线观看视频| 一级片免费观看大全| 黄色片一级片一级黄色片| 法律面前人人平等表现在哪些方面| 一级毛片精品| 一边摸一边抽搐一进一出视频| 国产在线一区二区三区精| 热99久久久久精品小说推荐| 免费观看av网站的网址| 91国产中文字幕| 成人国产av品久久久| 女同久久另类99精品国产91| 国产野战对白在线观看| 99在线人妻在线中文字幕 | 国产欧美亚洲国产| 久久久久久久久久久久大奶| 99精品久久久久人妻精品| 国产亚洲午夜精品一区二区久久| 色播在线永久视频| 亚洲中文字幕日韩| 国产精品香港三级国产av潘金莲| 日韩视频一区二区在线观看| 亚洲人成电影观看| 美女扒开内裤让男人捅视频| 精品少妇内射三级| 久久久久久免费高清国产稀缺| 波多野结衣av一区二区av| 少妇被粗大的猛进出69影院| 老司机影院毛片| 成年人黄色毛片网站| 亚洲国产中文字幕在线视频| 一区二区三区乱码不卡18| 日韩中文字幕欧美一区二区| 国产淫语在线视频| 一个人免费在线观看的高清视频| 亚洲五月婷婷丁香| 亚洲欧洲日产国产| 国产亚洲精品一区二区www | av不卡在线播放| 午夜精品国产一区二区电影| 一进一出好大好爽视频| 91成年电影在线观看| 汤姆久久久久久久影院中文字幕| 成年女人毛片免费观看观看9 | 丝袜美腿诱惑在线| 视频在线观看一区二区三区| 日韩大片免费观看网站| 国产在线精品亚洲第一网站| 三级毛片av免费| 国产黄色免费在线视频| 十八禁网站网址无遮挡| 亚洲专区国产一区二区| 一本色道久久久久久精品综合| 91国产中文字幕| 99精品久久久久人妻精品| 久久人妻av系列| 日韩中文字幕欧美一区二区| 欧美精品高潮呻吟av久久| 色94色欧美一区二区| tocl精华| 欧美久久黑人一区二区| 精品人妻1区二区| 免费观看人在逋| 69精品国产乱码久久久| 国产亚洲欧美精品永久| 欧美 亚洲 国产 日韩一| 欧美日韩亚洲高清精品| 午夜福利,免费看| 高清毛片免费观看视频网站 | 久久 成人 亚洲| 午夜两性在线视频| 国产高清激情床上av| 高清黄色对白视频在线免费看| 女人精品久久久久毛片| 色综合婷婷激情| 老司机靠b影院| 午夜福利,免费看| 欧美精品一区二区大全| 亚洲伊人久久精品综合| 亚洲 欧美一区二区三区| 国产精品99久久99久久久不卡| 少妇精品久久久久久久| 久久国产精品大桥未久av| 最黄视频免费看| 最新美女视频免费是黄的| 成年版毛片免费区| 久久国产精品男人的天堂亚洲| cao死你这个sao货| 精品久久蜜臀av无| www日本在线高清视频| 欧美精品高潮呻吟av久久| 国产精品.久久久| 欧美人与性动交α欧美精品济南到| 多毛熟女@视频| 欧美国产精品va在线观看不卡| 免费看a级黄色片| 黑丝袜美女国产一区| 免费不卡黄色视频| 亚洲欧美一区二区三区久久| 精品人妻1区二区| 日韩欧美免费精品| 久久这里只有精品19| 欧美一级毛片孕妇| 咕卡用的链子| 国产成人系列免费观看| 亚洲九九香蕉| a级毛片在线看网站| 最新在线观看一区二区三区| 久久人人爽av亚洲精品天堂| 大片免费播放器 马上看| av线在线观看网站| 中文字幕另类日韩欧美亚洲嫩草| 色综合欧美亚洲国产小说| 咕卡用的链子| 无遮挡黄片免费观看| 亚洲欧美一区二区三区黑人| 欧美成人免费av一区二区三区 | 日本一区二区免费在线视频| 一级毛片精品| 午夜免费成人在线视频| 亚洲中文av在线| 欧美人与性动交α欧美精品济南到| 亚洲欧美日韩高清在线视频 | 亚洲天堂av无毛| 在线观看免费高清a一片| 丝瓜视频免费看黄片| 91av网站免费观看| 精品欧美一区二区三区在线| 曰老女人黄片| 他把我摸到了高潮在线观看 | 亚洲自偷自拍图片 自拍| 免费黄频网站在线观看国产| 亚洲欧美一区二区三区黑人| 自拍欧美九色日韩亚洲蝌蚪91| 黄色视频在线播放观看不卡| avwww免费| 这个男人来自地球电影免费观看| 正在播放国产对白刺激| 亚洲精品中文字幕在线视频| 日本黄色日本黄色录像| 日本一区二区免费在线视频| 99九九在线精品视频| 巨乳人妻的诱惑在线观看| 久久久久国内视频| 久久99热这里只频精品6学生| 久久人人97超碰香蕉20202| 香蕉丝袜av| av一本久久久久| 国产一区二区三区视频了| av有码第一页| 亚洲成国产人片在线观看| 九色亚洲精品在线播放| 精品少妇黑人巨大在线播放| 中文字幕人妻熟女乱码| 日日摸夜夜添夜夜添小说| 新久久久久国产一级毛片| 91av网站免费观看| 久久久精品94久久精品| 蜜桃国产av成人99| 国产av精品麻豆| 国产精品免费大片| 亚洲三区欧美一区| 亚洲人成77777在线视频| 中文字幕色久视频| 久久久久久免费高清国产稀缺| 亚洲熟妇熟女久久| 一本色道久久久久久精品综合| 在线 av 中文字幕| 亚洲视频免费观看视频| 欧美日韩福利视频一区二区| 国产精品久久久久成人av| 欧美老熟妇乱子伦牲交| 午夜老司机福利片| 成年版毛片免费区| 国产欧美日韩精品亚洲av| 亚洲成a人片在线一区二区| 在线 av 中文字幕| 人妻一区二区av| 国产区一区二久久| 麻豆成人av在线观看| 99在线人妻在线中文字幕 | 老鸭窝网址在线观看| 亚洲 欧美一区二区三区| 免费看a级黄色片| 午夜福利在线观看吧| 日日夜夜操网爽| 欧美在线黄色| 亚洲色图 男人天堂 中文字幕| 黑人操中国人逼视频| a级片在线免费高清观看视频| 人妻久久中文字幕网| 精品亚洲成a人片在线观看| 嫁个100分男人电影在线观看| 50天的宝宝边吃奶边哭怎么回事| 色在线成人网| 免费人妻精品一区二区三区视频| 汤姆久久久久久久影院中文字幕| 狂野欧美激情性xxxx| 日韩一区二区三区影片| 日韩熟女老妇一区二区性免费视频| 精品久久蜜臀av无| 一二三四在线观看免费中文在| 日韩三级视频一区二区三区| 久久久久视频综合| 亚洲精品中文字幕在线视频| 国产男女超爽视频在线观看| 欧美国产精品一级二级三级| 好男人电影高清在线观看| 国产精品久久久久久人妻精品电影 | 91大片在线观看| 嫩草影视91久久| 欧美精品一区二区大全| 精品一区二区三区av网在线观看 | 亚洲第一青青草原| 日本vs欧美在线观看视频| 亚洲专区国产一区二区| 日韩熟女老妇一区二区性免费视频| 国产精品免费大片| 亚洲专区中文字幕在线| cao死你这个sao货| av国产精品久久久久影院| 午夜福利乱码中文字幕| 免费在线观看黄色视频的| 国产精品国产高清国产av | 国产高清激情床上av| 欧美日韩视频精品一区| 精品久久久精品久久久| kizo精华| 日韩大码丰满熟妇| 久久久久久久久免费视频了| 国产精品麻豆人妻色哟哟久久| 免费观看a级毛片全部| 亚洲欧美精品综合一区二区三区| 叶爱在线成人免费视频播放| 精品第一国产精品| 精品一区二区三卡| 亚洲欧美激情在线| 亚洲av日韩在线播放| 精品少妇内射三级| 欧美黑人欧美精品刺激| 一级片'在线观看视频| 好男人电影高清在线观看| a在线观看视频网站| 五月天丁香电影| av有码第一页| 国产日韩一区二区三区精品不卡| 搡老岳熟女国产| 国产成人一区二区三区免费视频网站| 黄网站色视频无遮挡免费观看| 国产福利在线免费观看视频| 国产成人精品久久二区二区免费| 国产91精品成人一区二区三区 | 免费在线观看影片大全网站| 母亲3免费完整高清在线观看| 一本久久精品| 国产精品九九99| 在线播放国产精品三级| 久久久国产欧美日韩av| 亚洲人成电影免费在线| 亚洲欧美日韩另类电影网站| 亚洲欧美一区二区三区黑人| 色综合婷婷激情| 老汉色av国产亚洲站长工具| 国产成人啪精品午夜网站| 精品亚洲成a人片在线观看| 波多野结衣一区麻豆| 国产精品亚洲一级av第二区| 欧美精品一区二区大全| 亚洲专区字幕在线| 久久精品人人爽人人爽视色| 午夜免费鲁丝| 另类亚洲欧美激情| 午夜久久久在线观看| 少妇精品久久久久久久| 日本wwww免费看| 午夜视频精品福利| 国产精品偷伦视频观看了| 手机成人av网站| 国产精品一区二区在线观看99| 国产日韩一区二区三区精品不卡| 男女无遮挡免费网站观看| 亚洲人成伊人成综合网2020| 国产一区二区三区视频了| 十八禁人妻一区二区| 777久久人妻少妇嫩草av网站| 狠狠婷婷综合久久久久久88av| 王馨瑶露胸无遮挡在线观看| videosex国产| 久久99热这里只频精品6学生| 国产精品久久久人人做人人爽| 少妇的丰满在线观看| 新久久久久国产一级毛片| 99久久国产精品久久久| 国产无遮挡羞羞视频在线观看| 亚洲av国产av综合av卡| 久久久久网色| 18禁美女被吸乳视频| 91九色精品人成在线观看| 国产精品九九99| 精品一区二区三区av网在线观看 | 亚洲av电影在线进入| 脱女人内裤的视频| 国产一卡二卡三卡精品| 免费久久久久久久精品成人欧美视频| 国产日韩欧美亚洲二区| 黄色片一级片一级黄色片| 黄色 视频免费看| 91国产中文字幕| 一级黄色大片毛片| 一级,二级,三级黄色视频| 亚洲专区国产一区二区| 精品福利永久在线观看| 国产精品熟女久久久久浪| 中文字幕另类日韩欧美亚洲嫩草| 精品国产国语对白av| 色94色欧美一区二区| 免费一级毛片在线播放高清视频 | 精品视频人人做人人爽| 老司机午夜福利在线观看视频 | 国产亚洲欧美在线一区二区| 亚洲精品久久午夜乱码| 91精品三级在线观看| 淫妇啪啪啪对白视频| 男女床上黄色一级片免费看| 国产一区二区 视频在线| 国产精品免费一区二区三区在线 | 日韩人妻精品一区2区三区| 黄色视频在线播放观看不卡| 人成视频在线观看免费观看| 在线十欧美十亚洲十日本专区| 国产成人精品久久二区二区91| 久久久久网色| 最近最新免费中文字幕在线| 亚洲精品在线美女| 中文字幕人妻熟女乱码| 伊人久久大香线蕉亚洲五| 在线观看免费午夜福利视频| 欧美+亚洲+日韩+国产| 亚洲成a人片在线一区二区| 无遮挡黄片免费观看| av网站在线播放免费| 国产免费av片在线观看野外av| 最新美女视频免费是黄的| 99精品欧美一区二区三区四区| 亚洲 欧美一区二区三区| 国产亚洲欧美在线一区二区| 黑人操中国人逼视频| 男女之事视频高清在线观看| 欧美日韩av久久| 99riav亚洲国产免费| 国产无遮挡羞羞视频在线观看| 久9热在线精品视频|