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

        ?

        無功、負(fù)序及諧波電流檢測中的數(shù)字濾波器研究

        2011-07-25 07:03:18張曉濱鐘彥儒王小艷
        電工技術(shù)學(xué)報(bào) 2011年8期
        關(guān)鍵詞:延時諧波濾波器

        張曉濱 鐘彥儒 王小艷

        (西安理工大學(xué)自動化學(xué)院 西安 710048)

        1 引言

        非線性、沖擊性和不平衡負(fù)載的使用造成電網(wǎng)電壓和電流的無功、負(fù)序及諧波分量顯著增加。這會引起電網(wǎng)中其他設(shè)備損壞,并降低電網(wǎng)運(yùn)行效率。為了提高電能質(zhì)量,STATCOM、APF等補(bǔ)償設(shè)備被應(yīng)用[1-4]。這些設(shè)備基于瞬時無功理論進(jìn)行無功、負(fù)序及諧波電流檢測,并依據(jù)檢測結(jié)果進(jìn)行補(bǔ)償。因此檢測的速度與精度直接影響設(shè)備補(bǔ)償性能。電流檢測方法如圖1所示,其中使用的數(shù)字濾波器為低通濾波器,后面稱其為傳統(tǒng)濾波器[5-9]。傳統(tǒng)濾波器的精度和速度是一對矛盾,一方提高另一方必然降低。為解決該問題,文獻(xiàn)[10-11]提出對傳統(tǒng)濾波器的這兩個指標(biāo)進(jìn)行折中優(yōu)化。文獻(xiàn)[12-14]對傳統(tǒng)濾波器進(jìn)行改造,在保證精度的前提下提高速度。

        本文在各種電壓和電流情況下對電流檢測中的數(shù)字濾波器進(jìn)行分析,得出數(shù)字濾波器分類設(shè)計(jì)的想法,并優(yōu)化設(shè)計(jì)出三種數(shù)字濾波器。仿真和實(shí)驗(yàn)結(jié)果證明所設(shè)計(jì)的數(shù)字濾波器相對于傳統(tǒng)濾波器在保證精度的前提下有效提高了速度,且濾波器本身算法的復(fù)雜性并未顯著增加。

        圖1 基于瞬時無功理論的電流檢測方法Fig.1 The current detection based on instantaneous reactive power theory

        2 電流檢測過程中的數(shù)字濾波器分析

        2.1 電壓正弦情況

        電壓和電流都為正弦時進(jìn)行電流檢測不需要濾波器。電流非正弦時分三種情況對濾波器進(jìn)行討論。

        2.1.1電流畸變

        電流的瞬時值表達(dá)式如下所示:

        對于三相三線系統(tǒng),式中n≠3k+3(k=0~∞)。

        瞬時有功功率為

        式中,n1=3k+1,n2=3k+2(k=0~∞)。

        有功功率的直流量由基波正序電流產(chǎn)生,交流量的最低諧波為 150Hz,因此數(shù)字濾波器應(yīng)濾除所有高于150Hz的諧波。稱該濾波器為F150濾波器。

        2.1.2電流不平衡

        電流的瞬時值表達(dá)式如下所示:

        瞬時有功功率為

        有功功率的直流量由基波正序電流產(chǎn)生,交流量為100Hz諧波,因此數(shù)字濾波器應(yīng)濾除100Hz諧波。稱該濾波器為F100濾波器。

        2.1.3電流畸變+不平衡

        電流的瞬時值表達(dá)式如下所示:

        式中,n=1~∞。

        瞬時有功功率為

        有功功率的直流量由基波正序電流產(chǎn)生,交流量的最低諧波為50Hz,因此數(shù)字濾波器應(yīng)濾除所有高于50Hz的諧波。稱該濾波器為F50濾波器。

        2.2 電壓非正弦情況

        當(dāng)電壓非正弦時,基于瞬時無功理論的電流檢測方法存在誤差[15]。對檢測方法進(jìn)行改進(jìn),如圖 2所示。改進(jìn)檢測方法主要是在原有方法基礎(chǔ)上增加了基波正序電壓的提取。用提取矩陣Txy/αβ和數(shù)字濾波器從電壓中得到基波正序電壓分量,然后以基波正序電壓分量為參考檢測出各個電流分量。

        圖2 改進(jìn)電流檢測方法Fig. 2 The improved current detection

        下面將電壓分三種情況討論提取電壓過程中使用的數(shù)字濾波器。

        2.2.1電壓畸變

        電壓的瞬時值表達(dá)式如下所示:

        式中,m≠3k+3(k=0~∞)。

        提取矩陣Txy/αβ為

        與提取矩陣相乘的運(yùn)算結(jié)果為

        式中,m1=3k+1,m2=3k+2(k=0~∞)。

        運(yùn)算結(jié)果的直流量為基波正序電壓,交流量的最低諧波為150Hz,因此F150濾波器滿足要求。

        2.2.2電壓不平衡

        電壓的瞬時值表達(dá)式如下所示:

        與提取矩陣相乘的運(yùn)算結(jié)果為

        運(yùn)算結(jié)果的直流量為基波正序電壓,交流量為100Hz諧波,因此F100濾波器滿足要求。

        2.2.3電壓畸變+不平衡

        電壓的瞬時值表達(dá)式如下所示:

        式中,m=1~∞。

        與提取矩陣相乘的運(yùn)算結(jié)果為

        運(yùn)算結(jié)果的直流量為基波正序電壓,交流量的最低諧波為50Hz,因此F50濾波器滿足要求。

        根據(jù)上面的分析發(fā)現(xiàn)在各種電壓和電流情況下,電流檢測需要的數(shù)字濾波器歸納起來共有三種:F150、F100和 F50。

        3 數(shù)字濾波器的分類優(yōu)化設(shè)計(jì)

        目前應(yīng)用的傳統(tǒng)濾波器通常設(shè)計(jì)為IIR形式的Butterworth低通濾波器[16-17]。通過分析可知在各種電壓和電流情況下,電流檢測需要的數(shù)字濾波器有三種。為了使這三種濾波器在各自適用的條件下,性能都優(yōu)于傳統(tǒng)濾波器,下面分別進(jìn)行優(yōu)化設(shè)計(jì)。

        3.1 F150濾波器的優(yōu)化設(shè)計(jì)

        F150濾波器要求濾除所有高于150Hz的諧波。與傳統(tǒng)濾波器相比,F(xiàn)150濾波器的截止頻率較高,因此速度較快。綜合考慮精度、速度以及運(yùn)算量,F(xiàn)150濾波器設(shè)計(jì)為直接Ⅱ型結(jié)構(gòu)的二階Butterworth低通濾波器,截止頻率為30Hz,通帶衰減≤2dB,阻帶衰減≤-15dB[18]。其傳遞函數(shù)為

        式中,a0=2.605 832e-4,a1=5.211 664e-4,a2=2.605 832e-4,b1=1.953 824,b2=0.954 866。

        3.2 F100濾波器的優(yōu)化設(shè)計(jì)

        F100濾波器要求濾除100Hz諧波,設(shè)計(jì)為直接Ⅱ型結(jié)構(gòu)的二階Butterworth帶阻濾波器,下阻帶截止頻率為97Hz,上阻帶截止頻率為103Hz,通帶衰減≤3dB,阻帶衰減≤-15dB。其傳遞函數(shù)為

        式中,c0=0.979 02,c1=1.949 18,c2=0.979 02,d1=1.949 18,d2=0.958 05。

        3.3 F50濾波器的優(yōu)化設(shè)計(jì)

        F50濾波器要求濾除所有高于50Hz的諧波,這與傳統(tǒng)濾波器要求相同。F50濾波器采用組合方式,用兩個帶阻濾波器和一個低通濾波器串聯(lián)。濾除50Hz諧波的帶阻濾波器設(shè)計(jì)為直接Ⅱ型結(jié)構(gòu)的二階 Butterworth帶阻濾波器,下阻帶截止頻率為47Hz,上阻帶截止頻率為53Hz,通帶衰減≤3dB,阻帶衰減≤-15dB。濾除100Hz諧波的帶阻濾波器和濾除所有高于150Hz諧波的低通濾波器分別采用F100濾波器和F150濾波器。其傳遞函數(shù)為

        式中,c0=0.974 48,c1=1.946 71,c2=0.974 48,

        d1=1.946 71,d2=0.948 96,e0=0.979 02,

        e1=1.949 18,e2=0.979 02,f1=1.949 18,

        f2=0.958 05,a0=2.605 832e-4,a1=5.211 664e-4,

        a2=2.605 832e-4,b1=1.953 824,b2=0.954 866。

        這三種數(shù)字濾波器的算法復(fù)雜性與傳統(tǒng)濾波器相比并無顯著增加,適合于工程應(yīng)用。

        4 仿真驗(yàn)證

        本文在 Matlab軟件上對所設(shè)計(jì)的數(shù)字濾波器進(jìn)行仿真驗(yàn)證。濾波器的輸入信號為瞬時有功功率,仿真中用一個無量綱的數(shù)值信號進(jìn)行代替。將F150濾波器與傳統(tǒng)濾波器進(jìn)行仿真對比,如圖3所示。輸入疊加直流量的150Hz交流信號。F150濾波器提取直流量的精度與傳統(tǒng)濾波器基本相同。當(dāng)直流量由 1(pu)跳變到 1.65(pu)時,傳統(tǒng)濾波器輸出延時約50ms,F(xiàn)150濾波器輸出延時約20ms,濾波速度明顯提高。

        圖3 F150濾波器與傳統(tǒng)濾波器的仿真對比Fig. 3 The simulation comparison between F150 filter and conventional filter

        將 F100濾波器與傳統(tǒng)濾波器進(jìn)行仿真對比,如圖4所示。輸入疊加直流量的100Hz交流信號。F100濾波器提取直流信號的精度與傳統(tǒng)濾波器基本一致。當(dāng)直流量由1(pu)跳變到1.65(pu)時,傳統(tǒng)濾波器輸出延時約 50ms,F(xiàn)100濾波器輸出延時約20ms,濾波速度明顯提高。

        圖4 F100濾波器與傳統(tǒng)濾波器的仿真對比Fig. 4 The simulation comparison between F100 filter and conventional filter

        將F50濾波器與傳統(tǒng)濾波器進(jìn)行仿真對比,如圖5所示。輸入疊加直流量的50Hz交流信號。F50濾波器提取直流信號的精度與傳統(tǒng)濾波器基本一致。當(dāng)直流量由 1(pu)跳變到 1.65(pu)時,傳統(tǒng)濾波器輸出延時約50ms,F(xiàn)50濾波器輸出延時約30ms,濾波速度明顯提高。

        圖5 F50濾波器與傳統(tǒng)濾波器的仿真對比Fig. 5 The simulation comparison between F50 filter and conventional filter

        圖6 F50濾波器的濾波性能仿真Fig. 6 The filtering performance simulation of F50 filter

        接下來對F50濾波器輸入疊加直流量的100Hz和150Hz交流信號。圖6顯示F50濾波器提取直流信號的精度與傳統(tǒng)濾波器基本一致。當(dāng)直流量發(fā)生跳變時F50濾波器輸出延時約30ms,濾波速度明顯提高。

        當(dāng)電網(wǎng)的電壓和電流為未知的非正弦時,瞬時有功功率所含的高次諧波成分未知,因此應(yīng)采用F50濾波器進(jìn)行濾波。分別應(yīng)用傳統(tǒng)濾波器和F50濾波器對瞬時有功功率進(jìn)行濾波,仿真結(jié)果如圖7所示。由圖中可見系統(tǒng)的瞬時有功功率在9s時有一個突變,傳統(tǒng)濾波器輸出延時約50ms,F(xiàn)50濾波器輸出延時約 30ms,濾波速度明顯提高。兩種濾波器的濾波精度基本一致。

        圖7 電壓和電流為非正弦時的仿真對比Fig. 7 The simulation comparison under non-sinusoidal voltage and current condition

        仿真結(jié)果證明所設(shè)計(jì)的三種數(shù)字濾波器在各自適用的條件下,其濾波精度與傳統(tǒng)濾波器基本相同,但濾波速度明顯提高。

        5 實(shí)驗(yàn)驗(yàn)證

        本文在TMS320F2407A上編程驗(yàn)證所設(shè)計(jì)的數(shù)字濾波器。實(shí)驗(yàn)中用一個無量綱的信號代替瞬時有功功率作為濾波器的輸入。將F150濾波器和傳統(tǒng)濾波器進(jìn)行實(shí)驗(yàn)對比,如圖 8所示。輸入疊加直流量的150Hz交流信號。F150濾波器提取直流信號的精度與傳統(tǒng)濾波器基本一致。當(dāng)直流量由 1(pu)跳變到 1.65(pu)時,傳統(tǒng)濾波器輸出延時約 50ms,F(xiàn)150濾波器輸出延時約 20ms,濾波速度明顯提高。

        圖8 F150濾波器與傳統(tǒng)濾波器的實(shí)驗(yàn)對比Fig. 8 The experiment comparison between F150 filter and conventional filter

        將 F100濾波器與傳統(tǒng)濾波器進(jìn)行實(shí)驗(yàn)對比,如圖9所示。輸入疊加直流量的100Hz交流信號。F100濾波器提取直流信號的精度與傳統(tǒng)濾波器基本一致。當(dāng)直流量由1(pu)跳變到1.65(pu)時,傳統(tǒng)濾波器輸出延時約 50ms,F(xiàn)100濾波器輸出延時約20ms,濾波速度明顯提高。

        圖9 F100濾波器與傳統(tǒng)濾波器的實(shí)驗(yàn)對比Fig. 9 The experiment comparison between F100 filter and conventional filter

        將F50濾波器與傳統(tǒng)濾波器進(jìn)行實(shí)驗(yàn)對比,如圖10所示。輸入疊加直流量的50Hz交流信號。F50濾波器提取直流信號的精度與傳統(tǒng)濾波器基本一致。當(dāng)直流量由 1(pu)跳變到 1.65(pu)時,傳統(tǒng)濾波器輸出延時約50ms,F(xiàn)50濾波器輸出延時約30ms,濾波速度明顯提高。

        圖10 F50濾波器與傳統(tǒng)濾波器的實(shí)驗(yàn)對比Fig. 10 The experiment comparison between F50 filter and conventional filter

        接下來對F50濾波器輸入疊加直流量的100Hz和 150Hz交流信號。圖11顯示F50濾波器提取直流信號的精度與傳統(tǒng)濾波器基本一致。當(dāng)直流量發(fā)生跳變時F50濾波器輸出延時約30ms,濾波速度明顯提高。

        圖11 F50濾波器的濾波性能實(shí)驗(yàn)Fig. 11 The filtering performance experiment of F50 filter

        當(dāng)電網(wǎng)的電壓和電流為未知的非正弦時,分別應(yīng)用傳統(tǒng)濾波器和F50濾波器對瞬時有功功率進(jìn)行濾波,實(shí)驗(yàn)結(jié)果如圖12所示。由圖中可見系統(tǒng)的瞬時有功功率在9s時有一個突變,傳統(tǒng)濾波器輸出延時約 50ms,F(xiàn)50濾波器輸出延時約 30ms,濾波速度明顯提高。兩種濾波器的濾波精度基本一致。

        圖12 電壓和電流為非正弦時的實(shí)驗(yàn)對比Fig. 12 The experiment comparison under non-sinusoidal voltage and current condition

        實(shí)驗(yàn)結(jié)果證明所設(shè)計(jì)的三種數(shù)字濾波器在各自的適用條件下,與傳統(tǒng)濾波器相比濾波精度基本相同,濾波速度明顯提高。

        6 結(jié)論

        應(yīng)用基于瞬時無功理論的檢測方法及其改進(jìn)方法進(jìn)行電流檢測。通過分析發(fā)現(xiàn)在各種電壓和電流情況下電流檢測需要的數(shù)字濾波器共有三種,由此得出數(shù)字濾波器分類設(shè)計(jì)的想法,并優(yōu)化設(shè)計(jì)出這三種濾波器。應(yīng)用仿真和實(shí)驗(yàn)將所設(shè)計(jì)的數(shù)字濾波器與傳統(tǒng)濾波器對比,結(jié)果表明所設(shè)計(jì)的數(shù)字濾波器在各自適用的條件下精度不變,速度明顯提高,且算法復(fù)雜性沒有顯著增加,適合工程實(shí)際應(yīng)用。

        [1] 王兆安, 楊君, 劉進(jìn)軍. 諧波抑制和無功功率補(bǔ)償[M]. 北京: 機(jī)械工業(yè)出版社, 1998.

        [2] 謝小榮, 姜齊榮. 柔性交流輸電系統(tǒng)的原理與應(yīng)用[M]. 北京: 清華大學(xué)出版社, 2006.

        [3] 劉定國, 羅安, 帥智康. 單周控制靜止同步補(bǔ)償器的死區(qū)補(bǔ)償策略[J]. 電工技術(shù)學(xué)報(bào), 2007, 22(12):159-163.

        Liu Dingguo, Luo An, Shuai Zhikang. Dead-time effect compensation method of unified constantfrequency integration static compensator[J]. Transactions of China Electrotechnical Society, 2007,22(12): 159-163.

        [4] 帥智康, 羅安, 劉定國, 等. 串聯(lián)諧振混合型有源電力濾波器及濾波特性分析[J]. 電工技術(shù)學(xué)報(bào),2009, 24(5): 125-134.

        Shuai Zhikang, Luo An, Liu Dingguo, et al. A novel hybrid active power filter with series resonance circuit[J]. Transactions of China Electrotechnical Society, 2009, 24(5): 125-134.

        [5] Akagi H. New trends in active filters for power conditioning[J]. IEEE Transactions on Industrial Application, 1996, 32 (6): 1312-1322.

        [6] 楊君, 王兆安. 三相電路諧波電流兩種檢測方法的對比研究[J]. 電工技術(shù)學(xué)報(bào), 1995, 10(2):43-48.

        Yang Jun, Wang Zhao’an. A study on the comparison of two methods used to detect the harmonic currents of three-phase circuits[J]. Transactions of China Electrotechnical Society, 1995, 10(2): 43-48.

        [7] Akagi H, Kanazawa Y, Nabae A. Instantaneous reactive power compensators comprising switching devices without energy storage components[J]. IEEE Transactions on Industrial Application, 1984,IA-20(3): 625-630.

        [8] 劉進(jìn)軍, 王兆安. 基于旋轉(zhuǎn)空間矢量分析瞬時無功功率理論及其應(yīng)用[J]. 電工技術(shù)學(xué)報(bào), 1999, 14(1):49-54.

        Liu Jinjun, Wang Zhao’an. Analysis of theory and application of instantaneous reactive power based on space vector method[J]. Transactions of China Electrotechnical Society, 1999, 14(1): 49-54.

        [9] 殷波, 陳允平, 鄧恒, 等. α-β坐標(biāo)系下瞬時無功理論與傳統(tǒng)功率理論的統(tǒng)一數(shù)學(xué)描述及物理意義[J].電工技術(shù)學(xué)報(bào), 2003, 18(5): 42-45.

        Yin Bo, Chen Yunping, Deng Heng, et al. Uniform mathematical description of instantaneous reactive power theory and conventional power theory and its physical meaning in α-β coordinates[J]. Transactions of China Electrotechnical Society, 2003, 18(5): 42-45.

        [10] 王群, 姚為正, 王兆安. 高通和低通濾波器對諧波檢測電路檢測效果的影響研究[J]. 電工技術(shù)學(xué)報(bào),1999, 14(5): 22-26.

        Wang Qun, Yao Weizheng, Wang Zhao’an. A study about influence of high and low pass filters on detecting effect of harmonics detection circuits[J].Transactions of China Electrotechnical Society, 1999,14(5): 22-26.

        [11] 王群, 姚為正, 王兆安. 低通濾波器對諧波檢測電路的影響[J]. 西安交通大學(xué)學(xué)報(bào), 1999, 33(4): 5-8.

        Wang Qun, Yao Weizheng, Wang Zhao’an. Effect of low pass filter on harmonics detection circuit[J].Journal of Xi’an Jiaotong University, 1999, 33(4): 5-8.

        [12] 周柯, 羅安, 夏向陽, 等. 一種改進(jìn)的 ip-iq諧波檢測方法及數(shù)字低通濾波器的優(yōu)化設(shè)計(jì)[J]. 中國電機(jī)工程學(xué)報(bào), 2007, 27(34): 96-101.

        Zhou Ke, Luo An, Xia Xiangyang, et al. An improved ip-iq harmonic current detection method and digital low-pass filter’s optimized design[J]. Proceedings of the CSEE, 2007, 27(34): 96-101.

        [13] 袁兆祥, 鄭建勇, 胡敏強(qiáng), 等. 有源濾波器中新型檢測濾波器的設(shè)計(jì)[J]. 中國電機(jī)工程學(xué)報(bào), 2008,28(9): 45-49.

        Yuan Zhaoxiang, Zheng Jianyong, Hu Minqiang, et al.Design of a novel detection filter of active power filter[J]. Proceedings of the CSEE, 2008, 28(9):45-49.

        [14] El-Habrouk M, Darwish M K, Mehta P. Active power filters: A review[J]. IEE Proceedings Electric Power Applications, 2000, 147(5): 403-413.

        [15] 劉翔宇, 楊仁剛, 王文成. 基于Ip-Iq法諧波檢測中數(shù)字低通濾波器的設(shè)計(jì)及其DSP實(shí)現(xiàn)[J]. 電力自動化設(shè)備, 2006, 26(8): 81-84.

        Liu Xiangyu, Yang Rengang, Wang Wencheng.Design of digital low pass filter in harmonic detection based on Ip-Iq transformation and its implementation using DSP[J]. Electric Power Automation Equipment,2006, 26(8): 81-84.

        [16] 田大強(qiáng), 蔣平, 唐國慶. 有源濾波器中數(shù)字低通濾波器設(shè)計(jì)及其實(shí)現(xiàn)[J]. 電工電能新技術(shù), 2003,22(3): 77-80.

        Tian Daqiang, Jiang Ping, Tang Guoqing. Design and implementation of digital filter in active power filter[J]. Advanced Technology of Electrical Engineering and Energy, 2003, 22(3): 77-80.

        [17] Zhang Bo, Yi Shongwen, He Xiaoming. A novel harmonic current detection technique based on a g eneralized dqk coordinate transform for active power filter and fault protection of power system[C].Proceedings of the International Conference on Advances in Power System Control, Operation &Measurement, Guangzhou, 2000.

        猜你喜歡
        延時諧波濾波器
        基于無擾濾波器和AED-ADT的無擾切換控制
        基于級聯(lián)步進(jìn)延時的順序等效采樣方法及實(shí)現(xiàn)
        從濾波器理解卷積
        電子制作(2019年11期)2019-07-04 00:34:38
        開關(guān)電源EMI濾波器的應(yīng)用方法探討
        電子制作(2018年16期)2018-09-26 03:26:50
        基于TMS320C6678的SAR方位向預(yù)濾波器的并行實(shí)現(xiàn)
        Two-dimensional Eulerian-Lagrangian Modeling of Shocks on an Electronic Package Embedded in a Projectile with Ultra-high Acceleration
        虛擬諧波阻抗的并網(wǎng)逆變器諧波抑制方法
        基于ELM的電力系統(tǒng)諧波阻抗估計(jì)
        基于ICA和MI的諧波源識別研究
        電測與儀表(2015年8期)2015-04-09 11:50:00
        桑塔納車發(fā)動機(jī)延時熄火
        少妇太爽高潮在线播放| 国产午夜影视大全免费观看| 无码毛片高潮一级一免费| 日韩精品少妇专区人妻系列| 91九色最新国产在线观看| 亚洲av久久久噜噜噜噜| 国产成人+亚洲欧洲+综合| 久久AⅤ天堂Av无码AV| 一区二区三区亚洲视频| 免费看男女做羞羞的事网站| 欧美变态口味重另类在线视频| 欧美zozo另类人禽交| 91国产自拍精品视频| 无码中文亚洲av影音先锋| 少妇厨房愉情理伦片bd在线观看 | 亚洲一区二区高清在线| 久久中文字幕一区二区| 无码国产69精品久久久久孕妇| 奇米狠狠色| 亚洲精品99久91在线| 亚洲综合极品美女av| 国产精品人妻一码二码尿失禁| 不卡国产视频| 亚洲综合久久精品少妇av | 精品久久久久久无码不卡 | 在线观看一区二区蜜桃| 国产专区一线二线三线码| 婷婷四房播播| 日本岛国一区二区三区| 西川结衣中文字幕在线| 中文字幕av无码一区二区三区| 日韩国产精品一本一区馆/在线 | 日本一区二区三区在线播放| 国产91色综合久久免费| 午夜成人鲁丝片午夜精品| 老汉tv永久视频福利在线观看 | 尤物yw午夜国产精品视频 | 亚洲色无码中文字幕| 国产乱精品女同自线免费| 99久久伊人精品综合观看| 91成人午夜性a一级毛片|