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

基于噪聲數據驅動CNN和LSTM的罐道故障診斷

2024-07-17 00:00:00馬利芬王偉池耀磊朱宏偉韓磊
無線電工程 2024年4期
關鍵詞:故障診斷振動故障

摘 要:剛性罐道是立井采礦中煤炭運輸以及設備運輸環節的重要組成機械,對罐道做出精確的故障診斷具有保證和提高立井采礦生產效率的重要意義。面對罐道難以實現精確故障診斷的問題,在Improved Convolution Neural Networks withTraining Interference (ITICNN) 的基礎上添加了Inception_v4 與Long Short Term Memory (LSTM),提出了一種高精度卷積神經網絡———Improved Convolution Neural Network Based on ITICNN (IBICNN),并采用了一種方向振動傳感器收集罐道振動信號,用IBICNN 的卷積層提取振動信號數字特征,并用LSTM 提取多方向振動信號之間的相關性信息的剛性罐道故障診斷方法。通過搭建實驗平臺對罐道實驗模型進行了故障診斷研究,取得了99. 4% 的診斷率。為了解決罐道在使用過程中伴有大量噪聲從而難以診斷故障的問題,在Adaptive Batch Normalization (AdaBN) 算法的基礎上進行了改進,并采取自遷移學習的方法,提高了IBICNN 的抗噪聲能力,在噪聲含量為100% 的情況下取得了90. 11% 的診斷率。

關鍵詞:罐道;故障診斷;IBICNN;LSTM;AdaBN

中圖分類號:TP277 文獻標志碼:A 開放科學(資源服務)標志碼(OSID):

文章編號:1003-3106(2024)04-1043-10

0 引言

隨著我國經濟建設發展,煤礦資源得到了廣泛開發,我國井下采礦礦場越來越多,因此立井提升運輸設備得到了大量的應用[1]。罐籠是井下采礦工程中的主要提升運輸設備,其主要功能是將地底的煤、機械和人員提升至地面,因此罐籠的安全穩定運行影響著采煤效率。剛性罐道作為罐籠的運行軌道,出現故障損傷時將影響到罐籠的安全與穩定運行。在采礦作業中由于機械振動、高腐蝕的環境往往引起罐道的損傷,當罐道發生故障時在不影響罐道使用的情況下對罐道提前做出故障診斷,對罐籠的使用具有重大安全意義[2]。

傳統的罐道故障診斷方法有傅里葉變換[3]、小波變換[4]、線性對應和模態分解[5]等,在罐道發生故障時,根據加速度傳感器采集到的振動數據的時域頻域特征,進行有監督的特征提取與分析從而識別出故障類型[6]。雖然傳統方法對罐道進行故障診斷比較直觀便于專業人員進行分析,但面對工況復雜、故障類型較多、故障數據較多以及噪聲干擾較大的問題時難以做出快速高效的故障診斷[7]。

隨著深度學習理論的不斷發展,各種神經網絡被廣泛應用于機械[8]、土木[9]和采礦[10]等各個領域[11]。張仲杭等[12]提出了一種基于紋理圖像與卷積神經網絡的軸承故障診斷方法,該方法將震動數據轉換成紋理圖像從而增加神經網絡的輸入尺度,進而增強了軸承故障診斷的精度。蔡超志等[13]提出了一種基于搭載Inception 模塊的卷積神經網絡對齒輪箱進行了故障診斷研究,研究結果表明在神經網絡中加入Inception 模塊可以提升網絡的抗噪聲能力和診斷精度,但其采用的Inception_v1 較為陳舊,性能提升有限。隨后蔡超志等[14]在TICNN基礎上采用了Inception 替換小卷積層,進而提出了一種高精度的卷積神經網絡———ImprovedConvolution Neural Networks with Training Interference(ITICNN)用于框架結構故障診斷研究并將該模型命名為BICNN,研究結果表明BICNN 不僅具備超高診斷精度而且具備較強的抗噪聲能力。此外,神經網絡具有強大的自適應特征提取能力與分類能力,因此本研究采取數據驅動的方法,在ITICNN 的基礎上加入Inception_v4[15]和Long Short Term Memory(LSTM)[16],提出了一種分類精度較高的卷積神經網絡———Improved Convolution Neural Networks withTraining Interference(IBICNN)。本研究采取IBICNN中的卷積層提取振動信號數字特征,用LSTM 提取加速度傳感器X、Y、Z 三個方向組合信號之間的相關性特征的綜合特征提取方法進行罐道故障診斷研究,對罐道的安全與穩定使用具有重大工程意義。為了對噪聲環境中的罐道做出精確的故障診斷,本研究使用帶噪聲的數據訓練IBICNN,從而提高IBICNN 對噪聲的容納范圍,并使用自遷移學習的Adaptive Batch Normalization(AdaBN)[17]方法,凍結BN 層參數,將帶噪聲振動數據作為源域進行正向傳播訓練,從而消除BN 層參數在源域和目標域形成相互獨立的關系,達到進一步提升抗噪聲能力的目的。

1 數據采集與增強

采取深度學習方法對罐道進行故障診斷時,采集數據位置的便利性尤為重要,選擇合適的位置不僅便于安裝傳感器,而且在罐道發生故障時便于快速地進行故障診斷,因此在實驗過程中將加速度傳感器放置在罐道頂端,由此來對應實際使用過程中罐道的頂部數據采取較為方便的特征,進而更加快速與準確地做出故障診斷。本實驗模擬了罐道在使用過程中的發生率最高的故障,包括傾斜、凸起、凹陷、錯位和間隙5 種故障類型。數據采集方式為每一次在罐籠開始提升時加速度傳感器采集不同方向的振動信號。實驗裝置與故障原理如圖1 所示,圖中①為加速度傳感器,②為提升裝置,③為罐道。

本實驗所采取的加速度傳感器可以收集來自X、Y、Z 三個方向的振動數據,所采集振動數據如圖2所示。

為了方便觀察,本研究采取了1 024 個點位進行了繪制。由圖2 可以看出,在時域內同一種故障情況下加速度傳感器采集到的3 個方向的振動信號雖然不同,但振幅區間相同,不同故障情況下采集到的同一方向振動信號凌亂程度具有較大差異。使用單一方向的震動數據進行故障診斷可以提取的特征數量有限,因此本研究采取對同一種故障的3 個方向的振動數據先進行數組增強保持其在數據增強中具有相同的步長,并且保持每個方向增強后的樣本順序不變,然后將每個方向的數據進行依次組合,最后將組合數據進行歸一化,其原理如圖3 所示。

經過銜接后的數組中具有不同方向的特征,因此使用神經網絡對其訓練可以提取到更加豐富的特征信息。在數據增強過程中,滑動窗口數據增強法因其本身高度的靈活性受到廣大研究學者的青睞,通過設置窗口的長度l 獲得樣本數據的長度,設置滑動步長b 可控制樣本的數量,b 的值越小得到的樣本重復性越高。滑動窗口數據增強原理如下:

式中:N 為數據增強后的樣本數量,L 為源數據點位數量。從式中可以看出,滑動窗口越大則會導致樣本數量越少,滑動步長b 越小則得到的樣本數量越多,b = 1 時可獲取最多的樣本數量為L-l。

數據歸一化有利于提升卷積神經網絡的訓練速度,并可以提升卷積神經網絡的收斂速度與訓練精度。振動數據歸一化常用的方法有線性歸一化以及零均值歸一化,線性歸一化和零值歸一化往往忽略振動信號中的負幅值,因此導致部分特征丟失,本研究將2 種方法結合起來對振動信號進行正負值歸一化,正負值歸一化原理如下:

式中:xi 為樣本x 中的第i 個點位的數值,yi 為樣本第i 個點位歸一化的數值。經過式(2)可以將樣本x 的點位全部歸一化到[-1,1],相當于將振動數據進行了等比例還原,因此保留了振動數據的負值特征。

2 罐道故障診斷流程

罐道往往被運用于噪聲較大的采礦運輸場所[18-19],因此使用卷積神經網絡對罐道進行故障診斷研究要求其具備良好的抗噪聲能力。使用不帶噪聲的振動數據訓練出來的卷積神經網絡具備一定的抗噪聲能力,但泛化能力較低。為了提高卷積神經網絡本身的抗噪聲能力,本研究采取帶一定噪聲的數據訓練神經網絡,從而提高卷積神經網絡的抗噪聲能力上限。本研究罐道故障診斷流程如圖4所示。

首先將數據劃分成訓練集、測試集和驗證集,使用不帶噪聲的訓練集和測試集訓練基于AdaBN 的卷積神經網絡,生成模型并保存。隨之往驗證集中逐步按照百分比添加10% ~ 100% 間隔為10% 的噪聲,并用不帶噪聲數據訓練出來的模型對帶噪聲的驗證集進行分類,當分類結果準確率小于等于90%且大于等于80% 時停止添加噪聲,依次探究出網絡模型最大的抗噪聲能力。將模型可承受的最大噪聲添加至訓練集和測試及并用于訓練卷積神經網絡,從而得出抗噪聲能力更強勢的卷積神經網絡模型。此外,往訓練數據中添加合適的噪聲去訓練卷積神經網絡,相比于傳統使用無噪聲信號訓練卷積神經網絡的方法,可容納噪聲范圍得到了增加,從而有利于高噪聲情況下罐道的故障診斷。

2. 1 卷積神經網絡

卷積神經網絡的性能往往取決于感受視野的大小和參數量,2022 年池耀磊[20]在ITICNN 的基礎上加入了Inception 模塊提出了精確度較高的BICNN,但其輸入樣本數據較短,在卷積神經網絡末端生成的特征參數量過低,因此在網絡訓練過程中損失函數值變化較大。為了對罐道做出精確的故障診斷,本研究在BICNN 的基礎上沿用了第一個卷積層為16 通道、卷積核大小為64、步長為16 的結構,使得卷積神經網絡具備較高的感受視野,并將原有的樣本輸入數據由1 024 增加為3 072,從而提高網絡末端的特征數量,此外將原有的Inception_v1 升級至Inception_v4,使得BICNN 增加了一個1×1 卷積的尺度并降低整體的參數量,從而得到IBICNN 結構。本研究所用的卷積神經網絡如圖5 所示。

IBICNN 首先由卷積核為64 的大卷積層對振動數據進行特征提取,第一層卷積的感受視野即為卷積核的尺寸,因此IBICNN 第一個卷積層可提取到較多的特征信息;其次,IBICNN 的第二個卷積核尺寸為32,用于細化提取第一層卷積輸出的特征參數;隨后,使用了5 組Inception_v4 結構,Inception_v4 中存在4 個尺度的卷積,且卷積核大小均不相同,從而可以提取到更加豐富的特征參數。

2. 2 LSTM

IBICNN 在訓練過程中卷積核僅僅提取當前區域的特征信息,因此忽略了前后卷積區域的相關性信息。此外,本研究采取了數據拼接,卷積層僅僅可以提取到每段信號的獨立信息,難以捕捉到拼接數據與數據之間的特征關系,因此在本研究在全連接層位置加入2 層LSTM 用于提取特征與特征之間的特征關系。LSTM 邏輯如圖6 所示。

圖6 為一個LSTM 單元,經過組合形成整體LSTM 層,xt 為t 時刻的輸入數據,ht 為經過LSTM后的輸出數據。LSTM 的更新公式如下:

it = σ(Wi[ht-1 ,xt]+bi), (3)

ot = σ(Wo [ht-1 ,xt]+bo ), (4)

ft = σ(Wf[ht-1 ,xt]+bf), (5)

Ct = ft ?Ct-1 +it ?tanh(Wc [ht-1 ,xt]+bc ), (6)

ht = ot ?tanh(Ct), (7)

式中:Wi、Wo、Wf 分別表示輸入門、輸出門、遺忘門的權重系數,bi、bo、bf 為偏置矢量。通過這樣的門結構,使得LSTM 神經網絡結構具有了保持長期儲存信息的能力,有效增加了記憶的長度,適用于提取組合信號之間的相關性特征提取。

2. 3 網絡參數

改進后的IBICNN 參數如表1 所示。

3 AdaBN

AdaBN 的原理為使用一個卷積神經網絡訓練形成的BN 參數去替換掉另一個神經網絡的BN 參數,從而防止神經網絡的訓練集和測試集在訓練過程中BN 層的均值和方差在空間領域分布為相互獨立,進而使得模型的泛化能力和魯棒性大大提高。AdaBN 算法的實現需要將IBICNN 在目標領域訓練的BN 層局部參數替換掉IBICNN 在源領域訓練所產生的BN 層的參數。本文在傳統AdaBN 算法的基礎上將BN 層所需要替換的參數提升為4 個,即進行整個BN 層參數替換,從而降低了AdaBN 算法的實現難度。AdaBN 算法如圖7 所示。

AdaBN 的具體實現方式為:首先,將數據分為等比例的訓練集和測試集,先用訓練集訓練IBICNN,訓練過程中再次進行7 ∶ 3 數據劃分,使得IBICNN 在訓練過程中進行正向和反向傳播并保存訓練好的模型;然后,使用訓練集的數據進行遷移學習,在遷移IBICNN 模型過程中凍結所有BN 層,使得BN 層不進行參數學習,同時在遷移學習的過程中只做正向傳播不做反向傳播;最終,保存遷移學習訓練好的模型用于故障診斷。AdaBN 算法流程如圖8 所示。

4 神經網絡訓練

本文本在一臺處理器為Corei78750H,GPU 為NVIDIA GTX1050ti4G,內存16 GB 的筆記本上運行IBICNN,噪聲來源為Tensorflow 中的numpy 函數模擬產生的高斯白噪聲。將不帶噪聲的數據輸入到IBICNN 的訓練結果如圖9 所示。

由圖9 可以看出,在源域和目標域的訓練過程中精確率值最終均穩定在100% 左右,但目標域的目標函數值在訓練170 ~ 200 次的跳躍程度明顯高于源域,由此可以說明源域訓練神經網絡穩定性更高。訓練達到200 次時聚類結果如圖10 所示。

由圖10 可以看出,目標域的訓練結果雖然可以將5 種故障數據分離,并且每種故障均收縮,但每種故障之間的距離接近,當信號中存在噪聲因素時會導致故障誤判。源域訓練結果可以將5 種故障數據完全分開各自歸攏,且故障數據與故障數據之間均存在較大的距離,因此當數據中含有一定的噪聲時,對源域訓練結果影響較低。

使用保存的IBICNN 模型對帶不同程度噪聲的驗證集數據進行故障分類,得到的分類精確度如圖11 所示。

由圖11 可以看出,當驗證集中噪聲含量達到50% 時,診斷精度降低到90% ,因此在原始數據中加入50% 噪聲用于訓練最終IBICNN 模型。訓練結果如圖12 所示。

對比圖12 和圖9 可以看出,往數據中加入50% 噪聲后IBICNN 的目標域訓練曲線更加粗糙、跳躍度更加明顯,但整體仍然趨于收斂,同時可以看出IBICNN 源域訓練受噪聲的影響不大,訓練曲線趨于一致,但加入噪聲后使得目標函數值發生高頻率的跳躍。單方面地分析訓練曲線難以看出模型是否具備分類條件,因此引入混淆矩陣和Tsne 可視化,結果如圖13 所示。

由圖13 可以看出,數據中加入50% 噪聲后,故障數據中存在2% 的情況4 被判斷成情況5,1% 的情況5 被判斷成情況4,但聚類結果中各種故障情況依舊各自歸攏且故障情況之間存在一定的距離,說明生成的模型仍然可以做出有效的故障診斷。通過混淆矩陣可計算出訓練的精確率,計算如下:

式中:TP 為預測成正的樣本數量,TN 為預測成負的樣本數量,TP+TF+FP+FN 為所有樣本數量。

5 消融實驗

當振動數據中存在50% 噪聲時,通過圖13 中的混淆矩陣可計算出本研究所提出的方法精度可以達到99. 4% ,為了探究本研究所改進網絡中各個組成部分的效果,使用TICNN、ITICNN、BICNN 以及IBICNN(AdaBN)對帶有不同噪聲的振動數據進行訓練,經過對比訓練結果曲線的精度分析出各個改進機制的作用。訓練對比結果如表2 所示。

對比表2 中ITICNN 和BICNN 可以看出,BICNN在噪聲含量為10% 時,精度高于ITICNN 1. 26% 左右,同時BICNN 的訓練結果曲線顯示精度均高于ITICNN,由此說明使用Inception 模塊可以增強卷積神經網絡的特征識別能力;對比BICNN 和IBICNN的目標域訓練結果可以看出,IBICNN 目標域訓練結果曲線精度均大于BICNN,由此可以說明本研究采用的Inception_v4 識別精度高于Inception_v1;對比IBICNN 源域和目標域精度,由于IBICNN 的源域訓練僅僅進行正向傳播,從而導致源域的訓練精度遠高于目標與訓練精度,同時可以說明采?。粒洌幔拢?算法可以提高網絡訓練的精度。經過對比分析結果可以得出本研采取Inception_v4 以及AdaBN 算法改進BICNN 可以提高網絡的識別能力的結論。

6 抗噪聲能力測試

使用帶噪聲數據訓練的IBICNN 模型對加入不同程度噪聲的驗證集進行分類,分類精度曲線如圖14 所示。

由圖14 可以看出,當源數據中含有60% 噪聲時分類精確度可以達到100% ,數據中含有100% 噪聲時分類精確度仍然可以達到90% ,滿足使用要求。對比圖9 和圖12 說明本研究所使用的方法可以大大提高故障診斷過程中的診斷精度。通過往驗證集數據中添加不同含量的高斯白噪聲,將SVM、MLP、TICNN、ITICNN、BICNN 以及本研究所提出方法的故障診斷的抗噪聲能力做對比,對比結果如表3 所示。

由表3 可以看出,傳統機器學習方法的SVM 診斷精度較低且精度降低速度較快,當數據中噪聲含量達到一定程度后診斷精度值變為1. 11% ,不再變化,說明網絡已經完全喪失分類能力。對比TICNN、ITICNN、BICNN 以及本研究所提出的方法可以看出,網絡結構越復雜得到的抗噪聲能力越強,當噪聲含量達到60% 時本研究所提方法仍然可以得到100% 的診斷率,當噪聲含量達到100% 時,本方法可以做出精度為90. 11% 的故障診斷,由此說明本研究所提出方法具有較強的抗噪聲能力和較高的診斷精度。

7 結束語

本研究采取了多軸加速度傳感器數據拼接的數據預處理方法,并在ITICNN 基礎上提出了一種高精度的IBICNN 卷積神經網絡,并利用AdaBN 算法,將IBICNN 訓練集與測試集的BN 層參數達到完全相互獨立,從而提升網絡本身的抗噪聲能力。此外,本研究采取了往罐道故障數據中逐步添加噪聲,當模型仍然可以做出90% 診斷率時,將噪聲添加進訓練數據再次訓練IBICNN,從而得到使用噪聲數據訓練出來的模型對噪聲環境下的罐道數據進行故障診斷。研究結果表明,當數據中含有60% 噪聲時使用本研究所提出的方法仍然可以達到100% 的診斷率,當數據中噪聲含量達到100% 時仍然具備90. 11% 的診斷精度。所提出的方法可用于噪聲環境下的罐道故障診斷,對罐道的使用具有重大安全意義。

參考文獻

[1] 吳波. 立井提升剛性罐道系統健康監測研究[D]. 徐州:中國礦業大學,2019.

[2] 季華,肖興明,劉正全. 罐道測試平臺的設計[J]. 煤礦機械,2007(11):25-27.

[3] 郟豪,尚嘯,張強. 航空發動機振動數據稀疏傅里葉變換分析[J]. 現代機械,2023(1):54-59.

[4] 黃姍姍,李志農. 基于高密度小波變換的航空發動機滾動軸承故障診斷方法[J]. 軸承,2023(2):19-25.

[5] 鄒寧,金楊超,郭成,等. 基于EMD 的井中雷達信號預處理[J]. 電子科技大學學報,2022,51(6):875-883.

[6] 馬艷芳,劉雪貞,鄧小飛. 刮板輸送機鏈傳動系統故障檢測與診斷[J]. 煤礦機械,2021,42(4):181-183.

[7] 胡偉康. 基于卷積神經網絡的立井剛性罐道健康監測研究[D]. 淮南:安徽理工大學,2021.

[8] 謝佳琪. 深度學習框架下的旋轉機械故障診斷研究[D]. 蘇州:蘇州大學,2020.

[9] ABDELJABER O,AVCI O,KIRANYAZ S,et al. Inman,Realtime vibrationbased Structural Damage Detection Using Onedimensional Convolutional Neural Networks[J]. Journal of Sound and Vibration,2017,388:154-170.

[10] 顧清華,江松,李學現,等. 人工智能背景下采礦系統工程發展現狀與展望[J]. 金屬礦山,2022(5):10-25.

[11] 邱穎豫,張柯,楊欣毅. 面向旋轉機械故障診斷的深度流形遷移學習[J]. 計算機工程與應用,2022,58(12):289-298.

[12] 張仲杭,蔡超志,徐輝,等. 基于紋理圖像的強抗干擾性軸承故障診斷方法[J]. 軸承,2023(8):79-86.

[13] 蔡超志,池耀磊,郭璐彬. 基于一種高精度卷積神經網絡的框架結構模型故障診斷研究[J]. 制造技術與機床,2022(1):135-140.

[14] 蔡超志,白金鑫,張仲杭等. 基于自適應小波降噪和In-ception 網絡的齒輪箱故障診斷[J]. 制造技術與機床,2022(10):21-28.

[15] 邢科,呂澤均. 基于InceptionV4 模型的通信信號調制方式識別新方法[J]. 現代計算機,2021(12):48-54.

[16] 熊志斌. 基于CEEMDAN 與LSTM 的人民幣匯率分析與預測[J]. 數理統計與管理,2022,41(3):507-525.

[17] ZHANG W,LI C H,PENG G L,et al. A Deep Convolutional Neural Network with New Training Methods for Bearing Fault Diagnosis Under Noisy Environment and Different Working Load [J ]. Mechanical Systems and Signal Processing,2018,100:439-453.

[18] 杜菲,馬天兵,胡偉康,等. 基于小波變換和改進卷積神經網絡的剛性罐道故障診斷[J]. 工礦自動化,2022,48(9):42-48.

[19] 馬天兵,王孝東,杜菲,等. 基于小波包和BP 神經網絡的剛性罐道故障診斷[J]. 工礦自動化,2018,44(8):76-80.

[20] 池耀磊. 基于深度學習理論的框架結構綜合故障診斷研究[D]. 邯鄲:河北工程大學,2022.

作者簡介

馬利芬 女,(1986—),碩士,高級工程師。主要研究方向:基于深度學習的采礦運輸設備智能故障診斷。

王 偉 男,(1974—),博士,高級工程師。主要方向:暖通設備智能化健康檢測與故障診斷。

池耀磊 男,(1996—),碩士,工程師。主要研究方向:基于深度學習理論的機械結構故障診斷。

朱宏偉 男,(1992—),碩士,講師。主要研究方向:深度學習理論在實際工程中的應用。

韓 磊 男,(1982—),碩士,工程師。主要研究方向:基于深度學習理論的電子技術。

基金項目:河北省重大科技成果轉化專項(22293601Z)

猜你喜歡
故障診斷振動故障
振動的思考
科學大眾(2023年17期)2023-10-26 07:39:14
振動與頻率
天天愛科學(2020年6期)2020-09-10 07:22:44
故障一點通
中立型Emden-Fowler微分方程的振動性
奔馳R320車ABS、ESP故障燈異常點亮
因果圖定性分析法及其在故障診斷中的應用
故障一點通
江淮車故障3例
UF6振動激發態分子的振動-振動馳豫
計算物理(2014年2期)2014-03-11 17:01:44
基于LCD和排列熵的滾動軸承故障診斷
主站蜘蛛池模板: 成年看免费观看视频拍拍| 老司国产精品视频| 亚洲国产精品日韩欧美一区| 国产精品亚洲天堂| 欧美一级专区免费大片| 欧美精品啪啪| 国产网友愉拍精品视频| 欧美一级在线看| 亚洲成aⅴ人在线观看| 狠狠色噜噜狠狠狠狠色综合久| 无码专区在线观看| 麻豆精品国产自产在线| 狠狠做深爱婷婷久久一区| 国产精品亚洲а∨天堂免下载| av手机版在线播放| 99九九成人免费视频精品| 国产精品嫩草影院av| 国产麻豆91网在线看| 亚洲精品视频网| 亚洲日本一本dvd高清| 国产一区二区三区在线观看视频 | 美女毛片在线| 四虎精品免费久久| 国产精品久久久久鬼色| 天堂成人av| 十八禁美女裸体网站| 一级毛片免费高清视频| 综合亚洲色图| 黑色丝袜高跟国产在线91| 在线观看国产精品日本不卡网| 亚洲国模精品一区| 欧美天堂在线| 国产理论一区| 狠狠做深爱婷婷久久一区| 日韩在线播放中文字幕| 中国美女**毛片录像在线| 视频一区视频二区中文精品| 在线观看精品自拍视频| 亚洲乱码在线视频| 国产精品无码在线看| 国产97视频在线观看| 亚洲免费三区| 91久久偷偷做嫩草影院精品| 久久婷婷综合色一区二区| 高清精品美女在线播放| 色首页AV在线| 久久人人97超碰人人澡爱香蕉| 在线观看免费国产| 免费不卡在线观看av| 久久这里只精品热免费99| www.日韩三级| 免费无码网站| 在线欧美日韩| 国产天天射| 成人日韩视频| 人人91人人澡人人妻人人爽 | 四虎成人在线视频| 片在线无码观看| 免费国产高清视频| 一级毛片免费不卡在线视频| 精品视频一区在线观看| 亚洲美女AV免费一区| 91精品国产自产91精品资源| 在线观看国产精品日本不卡网| 久久精品aⅴ无码中文字幕 | 中文字幕日韩丝袜一区| 国产鲁鲁视频在线观看| 欧美成人综合在线| 干中文字幕| 福利国产微拍广场一区视频在线| 日韩成人高清无码| 亚洲日本韩在线观看| 一级毛片在线播放| 国产激情无码一区二区APP| 日韩天堂视频| 亚洲第一精品福利| 欧洲免费精品视频在线| 亚洲精品欧美日本中文字幕| 国产精品手机视频| 无套av在线| 久久久受www免费人成| 一级成人欧美一区在线观看|