999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

基于排列熵算法的電力系統(tǒng)故障信號分析

2015-10-14 07:09:20李從善劉天琪李興源曹喜民劉利兵
電子科技大學(xué)學(xué)報 2015年2期
關(guān)鍵詞:故障診斷故障信號

李從善,劉天琪,李興源,曹喜民,劉利兵

?

基于排列熵算法的電力系統(tǒng)故障信號分析

李從善,劉天琪,李興源,曹喜民,劉利兵

(四川大學(xué)電氣信息學(xué)院 成都 610065)

電網(wǎng)故障的快速準確診斷對加快事故處理和系統(tǒng)恢復(fù)進程、保證電力系統(tǒng)的安全運行具有至關(guān)重要的作用。分析了目前應(yīng)用于電力系統(tǒng)故障診斷的數(shù)據(jù)來源,提出采用廣域量測的電氣量信息進行電力系統(tǒng)的快速診斷方法。首先將廣域量測系統(tǒng)采集的數(shù)據(jù)進行小波變換,通過相關(guān)系數(shù)法濾噪處理,以此來重構(gòu)信號,采用排列熵算法對重構(gòu)的信號進行排列熵分析。由于排列熵反映了一維時間序列的復(fù)雜度,對信號變化具有較高的敏感性,可以應(yīng)用于電力系統(tǒng)故障信號分析方面。相比傳統(tǒng)的小波分析方法,該方法不存在選取小波基問題,且算法簡單,編程易實現(xiàn),能夠滿足在線故障診斷條件。通過對IEEE10機39節(jié)點標準測試系統(tǒng)的故障仿真分析,結(jié)果表明了該方法的有效性和實用性。

相關(guān)系數(shù); 故障診斷; 排列熵算法; 電力系統(tǒng); 小波變換

大規(guī)模電網(wǎng)互聯(lián)在提高系統(tǒng)運行經(jīng)濟性的同時,也使得局部電網(wǎng)的故障對整個系統(tǒng)的影響范圍更廣。盡管現(xiàn)代電力系統(tǒng)采用了大量先進的設(shè)備、保護和控制手段,但仍無法避免電網(wǎng)故障的發(fā)生[1]。故障快速準確診斷是建設(shè)現(xiàn)代智能電網(wǎng)的必然要求,故障診斷歷來是專家和學(xué)者的研究熱點,并且已取得了一定的研究成果[2-17]。從已有的研究成果中可以看出,故障診斷主要涉及兩個方面的內(nèi)容:1) 故障診斷的信息來源;2) 故障診斷方法。目前應(yīng)用于故障診斷的信息來源主要有數(shù)據(jù)監(jiān)控及采集系統(tǒng)(SCADA)、繼電保護故障信息管理系統(tǒng)(RPMS)、廣域測量系統(tǒng)(WAMS)。已有的故障診斷方法主要有人工智能技術(shù)、解析模型方法以及信息融合技術(shù)3類方法,人工智能技術(shù)主要包括專家系統(tǒng)[2-3]、遺傳算法[4]、Petri網(wǎng)絡(luò)[5-6]、粗糙集理論[7-8]、貝葉斯網(wǎng)[9]、多Agent技術(shù)[10]、優(yōu)化技術(shù)[11]。文獻[12-15]通過建立故障解析模型并求解實現(xiàn)對故障的診斷。隨著信息理論的發(fā)展,研究者將信息融合技術(shù)引入到故障診斷中,文獻[16-17]提出了一種綜合SCADA開關(guān)量、故障錄波器電氣量及WAMS系統(tǒng)電氣量的多源信息融合電網(wǎng)故障診斷方法。文獻[18]基于振動信號對高壓斷路器機械故障診斷進行了研究。

當電力系統(tǒng)發(fā)生故障時,首先是故障元件的電氣量信息發(fā)生變化,其次是保護動作,最后才是相應(yīng)的開關(guān)[19]。采用電氣量信息可以最快實現(xiàn)系統(tǒng)故障診斷。隨著近年來WAMS量測系統(tǒng)在電力系統(tǒng)中的應(yīng)用,以相量測量單元(PMU)為基礎(chǔ)的WAMS量測數(shù)據(jù)具有高精度、延時小、更新周期短等優(yōu)點,這為電力系統(tǒng)故障在線診斷提供了數(shù)據(jù)平臺。本文在此數(shù)據(jù)平臺基礎(chǔ)上,提出采用排列熵理論實現(xiàn)電力系統(tǒng)故障診斷。首先將WAMS數(shù)據(jù)進行小波變換,通過相關(guān)系數(shù)濾噪處理,以此來重構(gòu)信號,采用排列熵算法對重構(gòu)的信號進行排列熵分析。與傳統(tǒng)的小波分析方法對比,該方法具有算法簡單,編程易實現(xiàn),不存在小波基選取困難問題,能夠滿足在線故障診斷的要求,這為電力系統(tǒng)故障信號分析提供了新的思路。通過對IEEE10機39節(jié)點標準測試系統(tǒng)的故障仿真分析,結(jié)果表明了該方法的有效性和實用性。

1 信號降噪處理

1.1 小波包分析

小波分析優(yōu)于傅里葉變換,它在時域和頻域同時具有良好的局部化性質(zhì),但每次只對信號的低頻部分進行分解,保留高頻部分。小波包分析對此進行了改進,可同時在低頻和高頻部分進行分解,自適應(yīng)地確定信號在不同頻段的分辨率,使之與信號頻譜相匹配,從而提高時頻分辨率[20]。以5層小波分解為例,兩者分解結(jié)構(gòu)如圖1所示。

由圖1可以看出,小波包分析是對小波分析的擴張,具有多維多分辨率分析的特點,能為信號提供一種更加精細的分析方法,它將頻帶進行多層次劃分,對多分辨率分析沒有細分的高頻部分進一步分解,提高了時頻分辨率。

a. 小波分析

b. 小波包分析

圖1 5層分解結(jié)構(gòu)圖

1.2 小波包分解與重構(gòu)算法

1.3 相關(guān)系數(shù)法濾波

一般情況下,包含噪聲的信號經(jīng)小波分解后,高頻分量即為噪聲。以往基于小波變換閾值降噪法的理論依據(jù)是[20],含噪聲的信號經(jīng)小波分解后有用信號的小波變換系數(shù)大于噪聲的小波變換系數(shù)。對分解得到的小波系數(shù)設(shè)定一個閾值,小波系數(shù)低于該閾值置為0,高于該閾值的小波系數(shù)完整保留或做相應(yīng)的“收縮”處理。最后將處理后獲得的小波系數(shù)利用逆小波變換進行重構(gòu),得出有效信號,實現(xiàn)降噪的目的。由于電網(wǎng)發(fā)生故障時,電氣量信號包含了大量的暫態(tài)高頻分量,若按照上述閾值降噪法,可能會將有用的故障特征成分一并剔除。

基于上述原因,本文提出相關(guān)系數(shù)法實現(xiàn)噪聲分量的篩選,其基本思想是噪聲信號與原信號的相關(guān)性很小,因此可以通過相關(guān)系數(shù)來辨別噪聲信號分量,相關(guān)系數(shù)定義為:

2 應(yīng)用于電力系統(tǒng)故障信號分析的排列熵算法

文獻[21-24]提出了一種衡量一維時間序列復(fù)雜度的平均熵參數(shù)——排列熵,它具有較高的魯棒性,能夠較好地反映時間序列數(shù)據(jù)的微小變化,目前主要應(yīng)用于氣象、醫(yī)學(xué)方面以及機械故障預(yù)測領(lǐng)域,在電力系統(tǒng)方面還未得到應(yīng)用。該方法與Kolmogorow熵和Lyapuunov指數(shù)相比具有概念簡單、計算簡單的優(yōu)點,可以更好地檢測出復(fù)雜系統(tǒng)的動力學(xué)突變,下面介紹其基本原理。

(6)

(8)

基于信號降噪與排列熵分析的電力系統(tǒng)故障診斷步驟為:

1) 首先將含噪信號經(jīng)小波分解,得到各層小波系數(shù);2) 設(shè)定一個閾值,根據(jù)式(1)計算各分量與原信號的相關(guān)性,若相關(guān)系數(shù)小于,則直接將其剔除;3) 將經(jīng)過篩選后的小波分量重構(gòu)得到濾波后的信號;4) 計算濾波后信號的排列熵,通過熵值診斷系統(tǒng)故障。

總體流程如圖2所示。

圖2 排列熵分析故障診斷流程

3 仿真分析

為驗證本文方法的有效性,以圖3所示的10機39節(jié)點標準測試系統(tǒng)為例,在電力系統(tǒng)仿真軟件PSCAD中,對系統(tǒng)發(fā)生概率最大的單相故障進行了仿真分析,并測得了系統(tǒng)故障后、、三相電流信號,將其做為診斷信號。

圖3 10機39節(jié)點系統(tǒng)

3.1 提取故障信號

通過在PSCAD里搭建10機39節(jié)點系統(tǒng),模擬單相瞬時性故障。不失一般性,本文選擇了節(jié)點3和節(jié)點18之間的線路發(fā)生相短路故障,故障發(fā)生在100 ms,持續(xù)50 ms消失。通過在真實信號(不含噪聲的信號)基礎(chǔ)上疊加信噪比為5的隨機噪聲,形成含量測噪聲的電氣量信號,如圖4a所示。

3.2 相關(guān)系數(shù)法去噪

采用小波變換對含噪聲的、、相電流信號進行5層小波分解,共得到32個小波分量,并用式(3)計算其相關(guān)系數(shù),以表示各個小波分量的相關(guān)系數(shù),以相電流為例,計算結(jié)果如表1所示。

表1 A相電流各分量相關(guān)系數(shù)

剔除小波相關(guān)系數(shù)小于0.09的小波分量,采用剩余的小波分量重構(gòu)信號,如圖4b所示。對比可知,經(jīng)過相關(guān)系數(shù)濾波去噪后能夠大致看出信號的變化,但不夠明顯。

3.3 排列熵分析

對采集的整個時間序列數(shù)據(jù),分別計算濾波后信號的、、三相的排列熵。計算結(jié)果如表2所示。

表2 排列熵值

從表2可以看出,相電流的排列熵明顯大于相和相。而、兩相電流的排列熵幾乎相等。從而由熵值可以判斷出相發(fā)生了故障。為了詳細分析整個時間序列排列熵的變化,采用文獻[24]的方法。將經(jīng)過濾噪后的整個時間序列數(shù)據(jù)分為若干個長度為的子序列,這些子序列相互之間使用最大重疊情形,即將每個子序列向后移動一個數(shù)據(jù)點得到下一個子序列。然后分別對、、三相電流信號計算其排列熵,仿真參數(shù),,。計算結(jié)果如圖5所示。

由圖5a可以明顯看出,相的排列熵在100 ms和150 ms發(fā)生了明顯突變,從而證明在這兩個時刻,電流發(fā)生了突變。而相和相排列熵值很平穩(wěn)沒有較大變化。因此從排列熵值的變化不僅可以判斷出系統(tǒng)故障,并且可以明顯看出故障發(fā)生的時刻。從而證明了該方法的有效性。

圖5 排列熵值變化

圖6 d1變化圖

為了做對比,本文又采用小波分析方法對信號做了仿真分析。采用db5小波,對信號進行6層分解,提取重構(gòu)信號的高頻部分1如圖6所示。通過對比,可得該方法與小波分析具有相同的結(jié)果。

4 結(jié) 論

本文提出將排列熵理論應(yīng)用于電力系統(tǒng)故障診斷中。通過小波包分解與相關(guān)系數(shù)法去噪后能夠較好地提取出故障的有用信息,并通過排列熵分析,不僅能夠診斷出系統(tǒng)故障,通過熵值的變化還可以得出故障發(fā)生的時刻,仿真結(jié)果驗證了該方法的有效性。本文是對排列熵在電力系統(tǒng)故障診斷中的有益嘗試,并取得了滿意的效果,其算法的簡單性和有效性可以滿足電力系統(tǒng)在線故障診斷的要求,為電力系統(tǒng)故障信號分析提供了新的思路,這對其在電力系統(tǒng)實際故障診斷中具有一定的參考意義。

[1] 劉道兵. 電網(wǎng)故障診斷的解析化建模與求解[D]. 北京: 華北電力大學(xué), 2012.

LIU Dao-bing. Analytical modeling and solving of power grid fault diagnosis[D]. Beijing: North China Electric Power University, 2012.

[2] 趙偉, 白曉民, 丁劍, 等. 基于協(xié)同式專家系統(tǒng)及多智能體技術(shù)的電網(wǎng)故障診斷方法[J]. 中國電機工程學(xué)報, 2006 26(20): 1-8.

ZHAO Wei, BAI Xiao-min, DING Jian, et al. A new fault diagnosis approach of power grid based on cooperative expert system and multi-agent technology[J]. Proceedings of the CSEE, 2006, 26(20): 1-8.

[3] LEE H J, AHN B S, PARK Y M. A fault diagnosis expert system for distribution substations[J]. IEEE Transactions on Power Delivery, 2000, 15(1): 92-97.

[4] LIU Yi, LI Ying, CAO Yi-jia, et al. Forward and backward models for fault diagnosis based on parallel genetic algorithms[J]. Journal of Zhejiang University Science A, 2008, 9(10): 1420-1425.

[5] 楊健維, 何正友, 臧天磊. 基于方向性加權(quán)模糊Petri網(wǎng)的電網(wǎng)故障診斷方法[J]. 中國電機工程學(xué)報, 2010, 30(34): 42-49.

YANG Jian-wei, HE Zheng-you, ZANG Tian-lei. Power system fault-diagnosis method based on directional weighted Fuzzy Petri nets[J]. Proceedings of the CSEE, 2010, 30(34): 42-49.

[6] SUN Jing, QIN Shi-yin, SONG Yong-hua. Fault diagnosis of electric power systems based on fuzzy Petri nets[J]. IEEE Transactions on Power Systems, 2004, 19(4): 2053-2059.

[7] 王英英, 羅毅, 涂光瑜. 基于粗糙集與決策樹的配電網(wǎng)故障診斷方法[J]. 高電壓技術(shù), 2008, 34(4): 794-798.

WANG Ying-ying, LUO Yi, TU Guang yu. Fault diagnosis method for distribution networks based on the rough sets and decision tree theory[J]. HighVoltage Engineering, 2008, 34(4): 794-798.

[8] 張耀天, 何正友, 趙靜, 等. 基于粗糙集理論和樸素貝葉斯網(wǎng)絡(luò)的電網(wǎng)故障診斷方法[J]. 電網(wǎng)技術(shù), 2007, 31(1): 37-43.

ZHANG Yao-tian, HE Zheng-you, ZHAO Jing, et al. A power network fault diagnosis method based on rough set theory and naive Bayesian networks[J]. Power System Technology, 2007, 31(1): 37-43.

[9] 周曙. 基于貝葉斯網(wǎng)的電力系統(tǒng)故障診斷方法研究[D]. 成都: 西南交通大學(xué), 2010.

ZHOU Shu. Power system fault diagnosis method based on Bayesian network research[D]. Chengdu: Southwest Jiaotong University, 2010.

[10] LEKKAS G P, AVOURIS N M, PAPAKONSTANTINOU G K. Development of distributed problem solving system for dynamic environments[J]. IEEE Transactions on System, Man and Cybemetics, 1995, 25(3): 400-414.

[11] 張曉玲, 陳先祿. 優(yōu)化技術(shù)在發(fā), 變電所接地網(wǎng)故障診斷中的應(yīng)用[J]. 高電壓技術(shù), 2000, 26(4): 64-66.

ZHANG Xiao-ling, CHEN Xian-lu.The technique of optimization applied in the grounding grid’s failure diagnosis of the power plant and substation[J]. High Voltage Engineering, 2000, 26(4): 64-66.

[12] 高湛軍, 陳青, 程婷婷, 等. 應(yīng)用于電網(wǎng)故障診斷的統(tǒng)一信息模型[J]. 電力系統(tǒng)自動化, 2010, 34(9):49-52.

GAO Zhan-jun, CHEN Qing, CHENG Ting-ting, et al. A Uniform information model for power grid fault diagnosis[J]. Auto matio n of Elect ric Power Systems, 2010, 34(9): 49-52.

[13] 翁漢琍, 毛鵬, 林湘寧. 一種改進的電網(wǎng)故障診斷優(yōu)化模型[J]. 電力系統(tǒng)自動化, 2007, 31(7): 66-70.

WENG Han-li, MAO Peng, LIN Xiang-ning. An improved model for optimizing power system fault diagnosis[J]. Auto matio n of Elect ric Power Systems, 2007, 31(7): 66-70.

[14] 郭文鑫, 廖志偉, 文福拴, 等. 計及警報信息時序特性的電網(wǎng)故障診斷解析模型[J].電力系統(tǒng)自動化, 2008, 32(22): 26-31.

GUO Wen-xin, LIAO Zhi-wei, WEN Fu-shuan, et al. An analytic model for power network fault diagnosis with the temporal informationof alarm messages taken in to account[J]. Auto Mation of Elect ric Power Systems, 2008, 32(22): 26-31.

[15] GUO Wen-xin, WEN Fu-shuan, LEDWICH G, et al. An analytic model for fault diagnosis in power systems considering malfunctions of protective relays and circuit breakers[J]. IEEE Transactions on Power Delivery, 2010, 25(3): 1393-1401.

[16] 劉毅, 高振興, 郭創(chuàng)新, 等. 一種考慮多層信息融合的電網(wǎng)故障診斷輔助決策方法[J]. 電力系統(tǒng)保護與控制, 2010, 38(24): 14-18.

LIU Yi, GAO Zhen-xing, GUO Chuan-xin, et al. An assistant decision-making for fault diagnosis based on multilayer information fusion[J]. Power System Protection and Control, 2010, 38(24): 14-18.

[17] 高振興, 郭創(chuàng)新, 俞斌, 等. 基于多源信息融合的電網(wǎng)故障診斷方法研究[J]. 電力系統(tǒng)保護與控制, 2011, 39(6): 17-23.

GAO Zhen-xing, GUO Chuan-xin, YU Bin, et al. Study of a fault diagnosis approach for powergrid with information fusion based on multi-data resources[J]. Power System Protection and Control, 2011, 39(6): 17-23.

[18] 關(guān)永剛, 黃瑜瓏. 基于振動信號的高壓斷路器機械故障診斷[J]. 高電壓技術(shù), 2000, 26(3): 66-68.

GUAN Yong-gang, HUANG Yu-long, QIAN Jia-li. An overview on mechanical failure diagnosis technique for high voltage circuit breaker based on vibration signal[J]. High Voltage Engineering, 2000, 26(3): 66-68.

[20] 符玲, 何正友, 麥瑞坤, 等. 小波熵證據(jù)的信息融合在電力系統(tǒng)故障診斷中的應(yīng)用[J]. 中國電機工程學(xué)報, 2008, 28(13): 64-69.

FU Ling, HE Zheng-you, MAI Rui-kun, et al. Information fusion method of entropy evidences and its application to fault diagnosis in power system[J]. Proceedings of the CSEE, 2008, 28(13): 64-69.

[21] 張翠芳. 基于小波變換的信號降噪技術(shù)的研究[D]. 南京: 南京郵電大學(xué), 2007.

ZHANG Cui-fang. Based on wavelet transform signal noise reduction technology research[D]. Nanjing: Nanjing University of Posts and Telecommunications, 2007.

[22] 侯威, 封國林, 董文杰, 等. 利用排列熵檢測近40年華北地區(qū)氣溫突變的研究[J]. 物理學(xué)報, 2006, 55(5): 2663- 2668.

HOU Wei, FENG Guo-lin, DONG Wen-jie, et al. A technique for distinguishing dynamical species in the temperature time series of north China[J]. Acta Phys Sin, 2006, 55(5): 2663-2668.

[23] 姚文坡, 劉鐵兵, 戴加飛, 等. 腦電信號的多尺度排列熵分析[J]. 物理學(xué)報, 2014, 63(7): 1-7.

YAO Wen-po, LIU Tie-bing, DAI Jia-fei, et al. Multiscale permutation entropy analysis of electroencephalogram[J]. Acta Phys Sin, 2014, 63(7): 1-7.

[24] 馮輔周, 饒國強, 司愛威. 排列熵算法研究及其在振動信號突變檢測中的應(yīng)用[J]. 振動工程學(xué)報, 2012, 25(2): 221-224.

FENG Fu-zhou, RAO Guo-qiang, SI Ai-wei, et al. Researchand application of the arithmetic of PE in testing the sudden change of vibration signal[J]. Journal of Vibration Engineering, 2012, 25(2): 221-224.

編 輯 漆 蓉

Power System Fault Signal Analysis Based on Permutation Entropy Algorithm

LI Cong-shan, LIU Tian-qi, LI Xing-yuan, CAO Xi-min, and LIU Li-bing

(School of Electrical Engineering and Information, Sichuan University Chengdu 610065)

Rapid and accurate diagnosis of fault has a crucial role for speeding up the recovery process and ensuring the safe operation of the power system. This paper analyzes the data sources currently used in power system fault diagnosis, proposes the use of wide-area information for rapid diagnosis of power systems. First, wavelet transform is used to process the wide area measurement data, through correlation coefficient method to eliminate the noise and reconstruct the signal, then permutation entropy calculation use the of the reconstructed signal. Since the complexity of the arrangement of entropy reflects a one-dimensional time series, the signal changes with high sensitivity, therefore, it can be applied to power system fault diagnosis. Compared with the traditional wavelet analysis method, this method avoids the problem of selecting wavelet base and meets the conditions for online fault diagnosis. In addition, it is simple and easy to program. The approach is applied in IEEE10 machine 39-bus system fault diagnosis. The diagnostic results show the applicability and effectiveness of the method.

correlation coefficient; fault diagnosis; permutation entropy (PE) algorithm; power system; wavelet transform

TM711

A

10.3969/j.issn.1001-0548.2015.02.013

2014-07-01;

2014-10-15

國家973項目(2013CB228204)

李從善(1985-),男,博士生,主要從事電力系統(tǒng)穩(wěn)定與控制方面的研究.

猜你喜歡
故障診斷故障信號
信號
鴨綠江(2021年35期)2021-04-19 12:24:18
完形填空二則
故障一點通
基于FPGA的多功能信號發(fā)生器的設(shè)計
電子制作(2018年11期)2018-08-04 03:25:42
奔馳R320車ABS、ESP故障燈異常點亮
基于LabVIEW的力加載信號采集與PID控制
因果圖定性分析法及其在故障診斷中的應(yīng)用
故障一點通
江淮車故障3例
基于LCD和排列熵的滾動軸承故障診斷
主站蜘蛛池模板: www成人国产在线观看网站| 国产成人精品亚洲77美色| 国产无吗一区二区三区在线欢| 性色生活片在线观看| 亚洲精品国产精品乱码不卞| 免费国产小视频在线观看| 精品国产www| 国产导航在线| 少妇精品网站| 亚洲第一成人在线| 四虎永久在线| 中文字幕亚洲精品2页| 激情亚洲天堂| 亚洲精品爱草草视频在线| 精品伊人久久大香线蕉网站| 久久黄色视频影| 国产福利在线观看精品| 欧美一区二区三区不卡免费| 996免费视频国产在线播放| 中文字幕人成人乱码亚洲电影| 国产在线视频福利资源站| 国产9191精品免费观看| 国产幂在线无码精品| 久久91精品牛牛| 国产日韩丝袜一二三区| 97超爽成人免费视频在线播放| 这里只有精品国产| 人妻21p大胆| 亚洲成a人片| 国产精品lululu在线观看| 成人在线亚洲| 中文字幕在线一区二区在线| 欧美精品一二三区| 日本欧美中文字幕精品亚洲| 午夜精品久久久久久久2023| 国产精品熟女亚洲AV麻豆| 国产精品对白刺激| 色噜噜综合网| 国产一在线观看| 亚洲V日韩V无码一区二区| 波多野结衣在线se| 色亚洲成人| 另类欧美日韩| 色精品视频| 丰满人妻一区二区三区视频| 国产精品成人啪精品视频| 国产精品香蕉在线观看不卡| av无码久久精品| 欧美在线中文字幕| 国产精品成人观看视频国产| 四虎影视8848永久精品| 国产精品尤物在线| 成人午夜免费视频| 永久免费无码日韩视频| 久久精品只有这里有| 亚洲国产av无码综合原创国产| 亚洲黄色成人| 日韩专区欧美| 国产一区二区三区在线精品专区 | 亚洲一区免费看| 亚洲欧洲日产无码AV| 最新加勒比隔壁人妻| 亚洲乱伦视频| 亚洲色图综合在线| 国产美女91呻吟求| 毛片在线区| 日韩精品成人在线| www亚洲精品| 国产欧美另类| 国产成熟女人性满足视频| 全部免费毛片免费播放| 国产微拍精品| 国产白浆一区二区三区视频在线| 亚洲免费三区| 国产黄在线观看| 麻豆精品在线视频| 国产毛片片精品天天看视频| 精品福利一区二区免费视频| 国产a v无码专区亚洲av| 国产理论最新国产精品视频| 四虎亚洲国产成人久久精品| 国产精选自拍|