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

    基于主成分SVM 的防腐層缺陷分類識別算法研究*

    2020-10-20 09:13:16呂瑞宏趙柏山楊佳怡
    微處理機 2020年5期
    關(guān)鍵詞:特征參數(shù)時域防腐

    呂瑞宏,趙 晗,趙柏山,楊佳怡

    (沈陽工業(yè)大學(xué)信息科學(xué)與工程學(xué)院, 沈陽110870)

    1 引 言

    管道防腐層缺陷檢測與識別是保障管道防腐層安全運行的重要研究內(nèi)容,利用超聲導(dǎo)波檢測防腐層缺陷具有明顯優(yōu)勢[1-3]。然而防腐層缺陷信號具有多模態(tài)、頻散特性,回波信號的識別、分析和處理難度較大[4-5]。為了更好地研究防腐層缺陷與回波信號的對應(yīng)關(guān)系,針對超聲導(dǎo)波回波信號的成分分析和處理尤為重要[6]。

    目前,針對檢測信號的識別分析和診斷主要集中在防腐層缺陷信號的時頻特征分析[7]上。通過對檢測信號進行時頻分析和模式識別,能夠有效實現(xiàn)信號的識別分效率以及準確率。胡明宇等[8]將小波包去噪和模式識別應(yīng)用于局部放電信號的研究;孫潔娣等[9]將LMD 局域均值分解應(yīng)用于天然氣管道泄漏特征的識別分析,取得較好的成果;姜久亮等[10]通過內(nèi)積延拓LMD 結(jié)合SVM 方法,對軸承故障診斷方法進行研究;蹇清平等[11]基于SVM 方法,對油管內(nèi)外表面缺陷進行了識別,有效提升了識別準確率;吳建鑫[12]等基于線性回歸將SVM 學(xué)習(xí)應(yīng)用于大規(guī)模分類;儲茂祥[13]基于雙支持向量機方法對鋼表面缺陷進行多信息分類,使其理論更加完善。

    基于上述研究,在此針對SVM 對處理大樣本缺陷回波信號大樣本特征參數(shù)、防腐層缺陷多分類問題和SVM 核函數(shù)選擇的局限性,以及SVM 處理多維回波信號特征參數(shù)的優(yōu)勢,提出基于主成分分析的SVM 防腐層缺陷分類方法。通過對防腐層回波信號進行局域均值分解,獲得噪聲主導(dǎo)與信號主導(dǎo)分量的轉(zhuǎn)折點,針對防腐層缺陷信號峰度、偏度、離散、形狀和小波包能量系數(shù)的特征矩陣,進行主成分分析降維SVM 模型特征指標的相關(guān)性,并引入動態(tài)擾動粒子群算法,解決傳統(tǒng)SVM 懲罰系數(shù)和核函數(shù)無法動態(tài)尋優(yōu)的問題,構(gòu)建基于主成分分析的SVM 動態(tài)粒子群優(yōu)化模型,實現(xiàn)管道防腐層缺陷類型的識別與分類。

    2 防腐層缺陷信號主成分分解模型

    針對超聲導(dǎo)波的頻散和多模態(tài)特性, 以及實驗檢測過程中引入的噪聲影響, 使得接收到的超聲回波信號非常復(fù)雜,需要對非線性超聲回波信號進行時頻分析。

    2.1 防腐層缺陷信號自適應(yīng)分解

    針對復(fù)雜、非線性、不平穩(wěn)、多分量的防腐層缺陷信號, 采用局域均值分解對缺陷回波信號進行時頻分析,將其分解為一系列回波信號分量之和?;夭ㄐ盘柗至慷加煞栏瘜尤毕莅j(luò)信號及防腐層缺陷純調(diào)頻信號相乘組成。通過回波信號分量的瞬時幅值和瞬時頻率, 可對防腐層缺陷信號進行自適應(yīng)分解,從而實現(xiàn)防腐層缺陷信號的時頻分布。

    對于防腐層缺陷信號進行均值處理,并求得其包絡(luò)值,如下二式所示:

    式中ni、ni+1為相鄰局域極值點。

    用直線連接所有缺陷均值點mi,并對連接后的均值波形進行平滑處理,從而得到缺陷均值函數(shù);用直線連接所有缺陷包絡(luò)值ai,對連接后的包絡(luò)值平滑處理,得到缺陷包絡(luò)函數(shù)。并從信號Y(t)中分離缺陷均值函數(shù),得到:

    式中,m11(t)為缺陷均值函數(shù)。對h11(t)解調(diào),可得:

    重復(fù)進行以上迭代,直到S1n(t)為理想條件下的純調(diào)頻缺陷信號。由于受到分解效果、迭代次數(shù)、速度等影響,選擇迭代停止條件為:a1n(t)≈1。迭代終止時,將此過程中產(chǎn)生的所有缺陷包絡(luò)函數(shù)相乘,得到初始回波信號分量的缺陷包絡(luò)信號,初始缺陷信號分量PF1(t)為:

    從信號Y(t)中分離初始缺陷信號分量,得到信號u1(t),將u1(t)作為初始缺陷信號重復(fù)以上迭代過程,直到uk(t)為單調(diào)函數(shù)。最后,防腐層缺陷回波信號X(t)被分解為k 個單分量與一個余量uk(t)之和:

    經(jīng)LMD 分解過的信號,會得到一系列的PF 分量,而噪聲主要存在于前面的高頻分量中,往后的大部分都是信號的主導(dǎo)部分,用連續(xù)均方誤差噪聲確定噪聲與信號主導(dǎo)分量之間的區(qū)分點:

    因在PF 分量分解過程中噪聲能量是逐步減少的,所以將PF 分量能量第一次出現(xiàn)轉(zhuǎn)折的地方定為區(qū)分點q,作為CMSE 的局部最小值點。

    首先對防腐層缺陷信號Y(t)進行LMD 分解,依次得出由高頻到低頻排列的PF 分量,并通過式(7)計算PF 分量的CMSE。通過獲取各CMSE 中的區(qū)分點q,將該點之前的PF 量提出并進行軟閾值去噪,并將去噪后的PF 分量與剩余的PF 分量和余量相加重構(gòu),則可獲得去噪后防腐層缺陷信號如下式:

    2.2 缺陷回波信號時頻域特征選取

    超聲導(dǎo)波與不同管道防腐層缺陷相互作用時,缺陷回波信號的時域特征參數(shù)會有明顯的差異。對于相同幅值的激勵信號,防腐層缺陷截面積越大,防腐層缺陷回波信號的幅值越高,因此有量綱的時域特征不能用來區(qū)分不同類型的管道防腐層缺陷。故此,選用峰度、偏度、離散和形狀系數(shù)四個時域無量綱特征參數(shù)在時域上區(qū)分管道防腐層缺陷類型。

    峰度系數(shù)(Ku)是用來描繪缺陷回波信號的尖銳程度的無量綱特征參數(shù),可以用來區(qū)分不同類型的防腐層缺陷。峰度系數(shù)的定義公式為:

    其中,N 為缺陷波形數(shù)據(jù)點的個數(shù),F(xiàn)i為缺陷波形在各數(shù)據(jù)點處的幅值,F(xiàn) 為缺陷波形平均幅值。

    偏度系數(shù)(Sk)被用來描述缺陷波形的分布正態(tài)性。由于超聲導(dǎo)波在傳播過程中與不同類型管道防腐層缺陷相互作用,會使得防腐層缺陷信號產(chǎn)生一定程度的偏斜,偏度系數(shù)可以用來區(qū)分不同類型的防腐層缺陷。偏度系數(shù)定義公式如下:

    離散系數(shù)(v)被用來描述防腐層缺陷信號的離散程度。離散系數(shù)越大表示防腐層缺陷信號的分布離散程度越大,可以此來區(qū)分不同類型的防腐層缺陷。離散系數(shù)定義公式如下式:

    形狀系數(shù)(Sc)是時域特征參數(shù)中能夠較好區(qū)分不同類型防腐層缺陷的特征參數(shù),能夠表征防腐層缺陷的跨度。形狀系數(shù)定義公式如下:

    其中,L 為防腐層缺陷信號所占據(jù)的長度,S 為防腐層缺陷信號與橫坐標之間的包絡(luò)面積,e 為防腐層缺陷信號的幅值。

    然而,在檢測過程中,由于作用于管道防腐層的超聲導(dǎo)波激勵頻率不同,僅通過時域特征參數(shù)難以表征不同類型的防腐層缺陷。而對不同防腐層缺陷通過使用相同頻率的激勵信號進行檢測時,缺陷回波信號的頻譜會存在差異,此時便不能僅根據(jù)缺陷回波信號頻譜進行分類,而是需要量化不同防腐層缺陷信號在頻譜上的差異。

    小波包分析可以對每一層缺陷信號所有頻率成分進行分解,從而獲得更高的防腐層缺陷信號時頻分辨率。從小波包分解后的防腐層缺陷信號提取頻段能量作為頻域特征可以反映防腐層缺陷信號在頻域上的能量分布差異。對防腐層缺陷信號進行小波包分解后,防腐層缺陷信號的能量被分解到不同的頻段,如果管道防腐層中存在缺陷,頻段能量會聚焦在防腐層缺陷處,則其他頻段的能量也會相應(yīng)的發(fā)生變化。將歸一化處理后的各頻段的能量構(gòu)成小波包缺陷能量譜特征集,歸一化后的頻域特征參數(shù)如下式所示:

    式中,za,b為第a 層上第b 個節(jié)點頻段的能量;Za為第a 層上各個頻段能量的總和。

    2.3 防腐層缺陷信號特征降維

    為了在時頻域充分反映防腐層缺陷信息,需提取多維度特征參數(shù),而多維度特征參數(shù)之間的相關(guān)性會對后續(xù)的缺陷識別分類產(chǎn)生影響。同時多維度的特征參數(shù)會對缺陷識別的計算過程造成負擔(dān),因此選用合適的特征降維方法能夠有效提高缺陷的分類準確程度。

    在此,選用峰度、偏度、離散和形狀系數(shù)四個時域無量綱特征參數(shù)在時域上區(qū)分管道防腐層缺陷類型,并對防腐層缺陷信號進行7 層小波包分解,選擇前20 個小波包能量系數(shù)在頻域上區(qū)分管道防腐層缺陷類型。

    防腐層缺陷信號特征向量矩陣為x(m·n),為了消除防腐層缺陷信號時頻域特征參數(shù)不同量綱以及大差異數(shù)值產(chǎn)生的影響,對特征參數(shù)進行標準化處理如下:

    式中,xj為缺陷特征均值,sj為缺陷特征標準差,i=1,2,...,m,j=1,2,...,n,為標準化后的缺陷特征。標準化特征參數(shù)的相關(guān)系數(shù)矩陣為:

    其中:

    求解矩陣T 的特征值和對應(yīng)的特征向量,并求解累計貢獻率:

    一般選取φ(j)在85%~95%的整數(shù),就確定了前j 個主成分來表征85%~95%原始數(shù)據(jù)信息,在保證信息完整的情況下,達到數(shù)據(jù)降維的目的。

    3 基于擾動粒子群優(yōu)化SVM 分類模型

    3.1 擾動粒子群算法參數(shù)尋優(yōu)

    粒子群算法通過粒子的速度迭代和位移迭代進行尋優(yōu),算法原理簡易,僅需調(diào)節(jié)較少的參數(shù),具有良好的記憶功能和參數(shù)優(yōu)化穩(wěn)定性,同時也具有優(yōu)良的全局收斂性和尋優(yōu)速率,但在尋優(yōu)過程中易出現(xiàn)陷入局部最優(yōu)參數(shù)值和未完成尋優(yōu)即輸出參數(shù)的情況,而且種群多樣性差,搜索范圍受限較大。這一情況在多維度參數(shù)尋優(yōu)過程中更為明顯。

    設(shè)在d 維空間進行參數(shù)尋優(yōu),m 個元素組成整個種群,第i 個元素位置為=(xi1,xi2,...,xid),速度為,其中i=1,2,...,m;最優(yōu)參數(shù)位置為,整個種群的最優(yōu)位置為;第i 個元素的速度按下式更新:

    第i 個粒子的位置更新公式為:

    其中,t 表示更新的代數(shù),c1和c2表示加速度系數(shù),最小值取0,最大值取2.5。此處選取c1=1.5,c2=1.7;r1和r2為在0 和1 之間任取的隨機數(shù),用來表示尋優(yōu)過程的記憶功能。

    為加快傳統(tǒng)參數(shù)尋優(yōu)過程的收斂速度、提高參數(shù)尋優(yōu)精度,并解決粒子群算法易陷入局部最優(yōu)參數(shù)值的問題,在此對算法的速度更新改進為:

    對于t 時刻粒子的速度更新公式,引入動態(tài)變化的速度更新權(quán)重,此權(quán)重的形式為指數(shù)函數(shù)與一個0 和1 之間任取隨機數(shù)的乘積,其定義公式為:

    式(22)在迭代初期,為了增強參數(shù)種群的全局尋優(yōu)能力,使ω 盡可能取較大值;到了迭代后期,為了提高參數(shù)種群的局部尋優(yōu)能力,使ω 盡可能取較小值。雖然在整個尋優(yōu)過程中,ω 整體保持下降趨勢,但并不要求前一次尋優(yōu)時的迭代慣性權(quán)重要比后一次的取值大。這樣提高了參數(shù)種群多樣性,還增強了粒子跳出局部最優(yōu)參數(shù)值的能力,有助于參數(shù)取值有效趨近全局最優(yōu)參數(shù)值。

    式(22)中δ 為尋優(yōu)擾動,考慮到全局最優(yōu)參數(shù)和局部最優(yōu)參數(shù)的值域相鄰處存在較優(yōu)參數(shù),采取在全局最優(yōu)參數(shù)值域搜索的過程中,在全局最優(yōu)參數(shù)和局部最優(yōu)參數(shù)中加入尋優(yōu)擾動,這樣可以避免尋優(yōu)過程陷入局部最優(yōu)參數(shù)值,提高參數(shù)在全局最優(yōu)參數(shù)的尋優(yōu)能力。經(jīng)過大量仿真研究發(fā)現(xiàn)δ=0.01時可有效提高局部尋優(yōu)和全局尋優(yōu)的效果。

    3.2 支持向量機缺陷分類

    支持向量機理論是以樣本統(tǒng)計進行模式識別為基礎(chǔ),強調(diào)以結(jié)構(gòu)風(fēng)險的最小化為前提,能夠解決例如神經(jīng)網(wǎng)絡(luò)等統(tǒng)計學(xué)習(xí)理論在模式識別方面的缺點。樣本統(tǒng)計學(xué)習(xí)方法也是針對高維數(shù)據(jù)統(tǒng)計學(xué)習(xí)方法的分類學(xué)習(xí)理論。此理論以高維數(shù)據(jù)統(tǒng)計方向?qū)π碌膶W(xué)習(xí)理論進行研究,在此種統(tǒng)計學(xué)習(xí)理論下的模式識別強調(diào)對于對漸近性能的提高,同時在高維度特征參數(shù)的分類條件下還能獲得最優(yōu)的預(yù)測分類精度。

    目前的理論研究主要集中在優(yōu)化基本SVM 算法中參數(shù)選擇的難點問題。傳統(tǒng)SVM 理論算法中的懲罰系數(shù)C 和RBF 核函數(shù)g 值, 對其基本要靠人工經(jīng)驗進行取值,這導(dǎo)致了學(xué)習(xí)過程過于繁雜,且缺乏理論根據(jù),通常測試集分類效果不佳。本研究針對基本SVM 算法中參數(shù)選擇的難點問題通過采用擾動粒子群優(yōu)化算法,對SVM 理論算法中所參數(shù)值進行尋優(yōu),使預(yù)測效果達到理想效果。

    擾動粒子群算法的支持向量機缺陷分類步驟為:首先進行參數(shù)初始化,設(shè)置種群大小,設(shè)置學(xué)習(xí)因子c1和c2,預(yù)置尋優(yōu)次數(shù)T、慣性權(quán)重的取值范圍及擾動范圍,設(shè)置參數(shù)C 和g 的取值范圍以及初始參數(shù)尋優(yōu)速度和參數(shù)位置;然后計算各次尋優(yōu)的適應(yīng)度值并計算每次尋優(yōu)的平均適應(yīng)度值,將適應(yīng)度值排序后求出個體和全局極值點,并對初始參數(shù)尋優(yōu)速度和參數(shù)位置進行更新;按式(22)進行每個參數(shù)的尋優(yōu)速度更新,按式(21)進行參數(shù)的位置更新,對參數(shù)進行自適應(yīng)的變異計算(變異率為10%),檢驗是否滿足尋優(yōu)終止條件;最后將C 和g 的尋優(yōu)結(jié)果輸入SVM 進行分類,通過對參數(shù)訓(xùn)練集進行基于統(tǒng)計學(xué)習(xí)參數(shù)訓(xùn)練,對測試數(shù)據(jù)集進行測試,得到防腐層缺陷分類結(jié)果,并對分類結(jié)果的準確率做分析。

    4 實驗研究

    4.1 實驗數(shù)據(jù)來源

    實驗對象為COMSOL 軟件預(yù)制的含若干防腐層缺陷的管道,管道為內(nèi)徑35mm、壁厚5mm 的結(jié)構(gòu)鋼, 防腐層為厚度4mm 的石英材料, 管道長度為300mm。對實驗對象進行超聲導(dǎo)波缺陷檢測,通過COMSOL 軟件處理并獲取檢測到的缺陷信號,共獲取178 組信號,其中裂紋60 組,凹坑59 組,孔洞59組,并使用MATLAB 進行仿真分析。含缺陷的管道防腐層建模如圖1 所示。

    圖1 含缺陷的管道防腐層COMSOL 建模

    仿真中加載頻率為2.5MHz、幅值為0.1mm 的正弦交變力作為激勵源。超聲波激勵信號為:

    加入式(24)所示的激勵信號,對三種波形的缺陷回波信號進行采集,得到缺陷回波信號的時域特征參數(shù)如表1 所示。

    表1 防腐層缺陷時域特征參數(shù)

    由表1 數(shù)據(jù)可知,在相同的激勵條件下,不同缺陷的回波信號在時域上的特征參數(shù)存在一定差異,時域特征參數(shù)可以作為表征不同缺陷的特征參數(shù)。

    4.2 LMD 缺陷信號時頻分析

    原超聲導(dǎo)波檢測信號及其經(jīng)LMD 局域均值分解的8 階PF 分量和殘余分量如圖2 所示。由圖可知,原信號X(t)由于在檢測過程中混入噪聲,在分解出的PF 分量中,隨著分解的階數(shù)升高,噪聲的混疊明顯減弱,其攜帶的原始信號信息隨之變得明顯。由于噪聲主要存在于前面的高頻分量中,往后的大部分都是信號的主導(dǎo)部分,用連續(xù)均方誤差噪聲確定噪聲與信號主導(dǎo)分量之間的區(qū)分點n 為4。實驗得出前4 階PF 分量為噪聲主導(dǎo)部分,5 階到8 階PF 分量為信號主導(dǎo)部分。LMD 分解前4 階噪聲主導(dǎo)PF 分量如圖2(a)所示;5~8 階信號主導(dǎo)PF 分量及余量如圖2(b)所示。

    圖2 LMD 局域均值分解PF 分量及余量

    由圖2(a)可知,對于防腐層缺陷回波信號進行LMD 分解,提取出的前4 階PF 分量雜亂,其噪聲主要存在于高頻成分中;由圖2(b)可知,對于防腐層缺陷回波信號進行LMD 分解,隨著分解階數(shù)的增大,信號對于缺陷信息的表征變得更為清晰,即信號主導(dǎo)部分分量。對前4 階PF 分量進行小波閾值去噪,結(jié)果如圖3 所示。

    圖3 對前4 階PF 分量小波閾值去噪

    將去噪后的PF 分量與剩余的PF 分量及余量相加重構(gòu),得到去噪后信號。原始與信號重構(gòu)信號如圖4 所示。由圖中對缺陷回波信號預(yù)處理前后的對比情況可知,對信號進行LMD 局域均值分解,通過連續(xù)均方誤差噪聲確定噪聲與信號主導(dǎo)分量之間的區(qū)分點,并對噪聲主導(dǎo)成分進行小波閾值去噪,相加重構(gòu)后,可以有效去除回波信號中的噪聲。

    圖4 防腐層缺陷原始信號與重構(gòu)信號

    4.3 SVM 缺陷分類結(jié)果分析

    實驗對缺陷信號的特征系數(shù)進行PCA 主成分分析,將提取的主成分通過支持向量機進行缺陷分類。缺陷類型分別為孔洞、凹坑和裂紋,并通過改進粒子群算法對SVM 方法中所涉及的懲罰系數(shù)C 和RBF 核函數(shù)g 值尋優(yōu),得到BestC=0.1,g=654.8023。詳細對比結(jié)果如表2 示。在178 組測試集中,共11組錯分,分類準確率達到了93.8%。

    表2 防腐層缺陷分類結(jié)果

    5 結(jié) 束 語

    利用非線性超聲導(dǎo)波對管道防腐層進行內(nèi)檢測,對回波信號進行時頻分析,采用主成分分析提取回波信號的特征參數(shù),并基于改進的粒子群算法優(yōu)化的支持向量機對缺陷類型進行分類。實驗結(jié)果表明,利用改進后的粒子群算法對SVM 方法中所涉及的懲罰系數(shù)C 和RBF 核函數(shù)g 值尋優(yōu),使其分類準確率達到了93.8%,相比傳統(tǒng)SVM 分類在準確率上有12.9%的提高,改善效果明顯,為管道防腐層缺陷檢測提供了一種有效的分類方法。

    猜你喜歡
    特征參數(shù)時域防腐
    故障診斷中信號特征參數(shù)擇取方法
    基于特征參數(shù)化的木工CAD/CAM系統(tǒng)
    基于防腐層電阻率的埋地管道防腐層退化規(guī)律
    煤氣與熱力(2021年3期)2021-06-09 06:16:26
    基于時域信號的三電平逆變器復(fù)合故障診斷
    基于PSO-VMD的齒輪特征參數(shù)提取方法研究
    基于極大似然準則與滾動時域估計的自適應(yīng)UKF算法
    壓力容器補強圈設(shè)計與防腐對策
    基于時域逆濾波的寬帶脈沖聲生成技術(shù)
    事故車維修中的車身防腐(三)
    統(tǒng)計特征參數(shù)及多分類SVM的局部放電類型識別
    電測與儀表(2015年7期)2015-04-09 11:40:04
    国产精品.久久久| 国产av一区二区精品久久| 高清不卡的av网站| 青春草亚洲视频在线观看| 99热国产这里只有精品6| 久久天躁狠狠躁夜夜2o2o | 亚洲国产毛片av蜜桃av| 国产免费一区二区三区四区乱码| 精品国产国语对白av| 午夜激情av网站| 免费在线观看影片大全网站 | 久久人妻福利社区极品人妻图片 | 好男人视频免费观看在线| 欧美成人午夜精品| 激情五月婷婷亚洲| 亚洲 欧美一区二区三区| 老司机亚洲免费影院| 成年人黄色毛片网站| 少妇人妻 视频| 国产又爽黄色视频| 国产免费视频播放在线视频| 美女扒开内裤让男人捅视频| 亚洲av男天堂| 久久99精品国语久久久| 日本黄色日本黄色录像| 亚洲第一av免费看| 又紧又爽又黄一区二区| av有码第一页| 国产欧美日韩一区二区三 | 麻豆国产av国片精品| 又大又爽又粗| 伊人亚洲综合成人网| 最近最新中文字幕大全免费视频 | 日日摸夜夜添夜夜爱| 亚洲国产欧美在线一区| 国产一区亚洲一区在线观看| 久久精品久久久久久噜噜老黄| 亚洲伊人久久精品综合| 日韩中文字幕视频在线看片| 亚洲色图 男人天堂 中文字幕| 亚洲一区二区三区欧美精品| 视频在线观看一区二区三区| 纵有疾风起免费观看全集完整版| 成人亚洲欧美一区二区av| 看十八女毛片水多多多| 日韩欧美一区视频在线观看| 欧美日韩精品网址| √禁漫天堂资源中文www| 国产一区二区三区av在线| 久久精品人人爽人人爽视色| 黄片小视频在线播放| 一边摸一边抽搐一进一出视频| 国产亚洲欧美精品永久| 五月开心婷婷网| 亚洲伊人色综图| 欧美日韩国产mv在线观看视频| 国产精品欧美亚洲77777| 中文欧美无线码| 国产男女超爽视频在线观看| 久久精品aⅴ一区二区三区四区| 亚洲五月色婷婷综合| 人人妻人人添人人爽欧美一区卜| 精品国产一区二区久久| 国产精品av久久久久免费| 嫩草影视91久久| 一级a爱视频在线免费观看| 免费一级毛片在线播放高清视频 | 一区在线观看完整版| 香蕉丝袜av| 国产午夜精品一二区理论片| 丝瓜视频免费看黄片| 美女视频免费永久观看网站| 亚洲国产中文字幕在线视频| 制服诱惑二区| 性色av一级| 亚洲 欧美一区二区三区| 国产成人一区二区三区免费视频网站 | 19禁男女啪啪无遮挡网站| 日韩 欧美 亚洲 中文字幕| 亚洲午夜精品一区,二区,三区| 国产男人的电影天堂91| 天堂中文最新版在线下载| 尾随美女入室| 男人舔女人的私密视频| 国产精品国产三级国产专区5o| 日本av手机在线免费观看| 成人国语在线视频| 精品熟女少妇八av免费久了| 香蕉国产在线看| 永久免费av网站大全| 亚洲,欧美,日韩| 国产深夜福利视频在线观看| 一级黄色大片毛片| 亚洲欧美精品综合一区二区三区| 黄色 视频免费看| 国产真人三级小视频在线观看| 人妻一区二区av| 亚洲一码二码三码区别大吗| 午夜av观看不卡| 亚洲av成人精品一二三区| 欧美久久黑人一区二区| 亚洲国产欧美网| 久久久精品区二区三区| videos熟女内射| 久久鲁丝午夜福利片| 国产免费福利视频在线观看| 别揉我奶头~嗯~啊~动态视频 | 亚洲伊人久久精品综合| 国产精品国产av在线观看| 一区二区日韩欧美中文字幕| 亚洲人成77777在线视频| 如日韩欧美国产精品一区二区三区| 国产高清videossex| 50天的宝宝边吃奶边哭怎么回事| 少妇 在线观看| 新久久久久国产一级毛片| 性高湖久久久久久久久免费观看| 香蕉国产在线看| 欧美少妇被猛烈插入视频| 老司机亚洲免费影院| 久久久国产欧美日韩av| 成年人免费黄色播放视频| 中文字幕另类日韩欧美亚洲嫩草| 亚洲人成电影免费在线| 欧美中文综合在线视频| 精品卡一卡二卡四卡免费| 亚洲熟女毛片儿| 亚洲国产看品久久| 国产免费福利视频在线观看| 又大又黄又爽视频免费| 一区二区三区四区激情视频| 在线观看www视频免费| 国产日韩欧美在线精品| 国产av精品麻豆| 国产麻豆69| 91老司机精品| 一边摸一边抽搐一进一出视频| 天天躁夜夜躁狠狠躁躁| av欧美777| 99re6热这里在线精品视频| 久久天堂一区二区三区四区| 男女之事视频高清在线观看 | 99热国产这里只有精品6| 国产日韩一区二区三区精品不卡| svipshipincom国产片| 美女扒开内裤让男人捅视频| 美国免费a级毛片| 国产高清视频在线播放一区 | 国产视频一区二区在线看| 久热这里只有精品99| 少妇裸体淫交视频免费看高清 | 黄网站色视频无遮挡免费观看| 久久久久久久国产电影| 久久精品aⅴ一区二区三区四区| 欧美日韩亚洲综合一区二区三区_| 男的添女的下面高潮视频| 青草久久国产| 桃花免费在线播放| 精品一区二区三区av网在线观看 | 亚洲成色77777| 最新的欧美精品一区二区| 欧美人与性动交α欧美精品济南到| 黄色片一级片一级黄色片| 一级毛片女人18水好多 | 国产伦人伦偷精品视频| 国产成人一区二区三区免费视频网站 | 国产97色在线日韩免费| 18禁黄网站禁片午夜丰满| 搡老乐熟女国产| 亚洲精品成人av观看孕妇| 精品视频人人做人人爽| 在线天堂中文资源库| 日本av免费视频播放| 亚洲av电影在线观看一区二区三区| 久久 成人 亚洲| 国产片内射在线| 又黄又粗又硬又大视频| 久久国产精品男人的天堂亚洲| 国产福利在线免费观看视频| 丰满人妻熟妇乱又伦精品不卡| 手机成人av网站| 欧美变态另类bdsm刘玥| 一级,二级,三级黄色视频| 精品一区在线观看国产| 免费看av在线观看网站| 国产亚洲精品第一综合不卡| 日日夜夜操网爽| 久久人人97超碰香蕉20202| 人人妻人人澡人人爽人人夜夜| 亚洲欧洲国产日韩| 好男人视频免费观看在线| 欧美+亚洲+日韩+国产| 热99久久久久精品小说推荐| 亚洲黑人精品在线| 亚洲成国产人片在线观看| 久久精品久久精品一区二区三区| 老汉色∧v一级毛片| 捣出白浆h1v1| 性少妇av在线| 一边摸一边抽搐一进一出视频| av欧美777| 国产精品av久久久久免费| 久久九九热精品免费| 丰满迷人的少妇在线观看| 天堂俺去俺来也www色官网| 国产爽快片一区二区三区| 老司机亚洲免费影院| 韩国精品一区二区三区| 在线天堂中文资源库| 亚洲国产欧美一区二区综合| 免费不卡黄色视频| 国产免费一区二区三区四区乱码| 超碰97精品在线观看| 久久天堂一区二区三区四区| 久久久久久久久久久久大奶| 亚洲人成电影免费在线| 亚洲av电影在线观看一区二区三区| 丰满饥渴人妻一区二区三| av视频免费观看在线观看| 天天操日日干夜夜撸| 久久久精品国产亚洲av高清涩受| 中国国产av一级| 老司机深夜福利视频在线观看 | 香蕉国产在线看| 两个人看的免费小视频| 亚洲欧美激情在线| 考比视频在线观看| 捣出白浆h1v1| 1024视频免费在线观看| 在线av久久热| 99久久综合免费| 韩国精品一区二区三区| 亚洲,欧美,日韩| 国产精品.久久久| 交换朋友夫妻互换小说| 国产色视频综合| www.精华液| 美国免费a级毛片| 欧美+亚洲+日韩+国产| 久热这里只有精品99| videosex国产| 最黄视频免费看| 激情五月婷婷亚洲| 国产老妇伦熟女老妇高清| 十八禁高潮呻吟视频| 亚洲第一青青草原| 亚洲第一青青草原| www.av在线官网国产| 90打野战视频偷拍视频| 成人黄色视频免费在线看| 精品国产一区二区久久| 国产精品三级大全| 亚洲国产最新在线播放| 美女国产高潮福利片在线看| 亚洲欧美一区二区三区国产| 97在线人人人人妻| 亚洲av日韩精品久久久久久密 | 国产日韩一区二区三区精品不卡| 国产老妇伦熟女老妇高清| 丁香六月天网| 老汉色av国产亚洲站长工具| 日韩大码丰满熟妇| 一边亲一边摸免费视频| 亚洲精品日本国产第一区| 精品国产一区二区三区四区第35| 在线看a的网站| 久久久久视频综合| 免费在线观看视频国产中文字幕亚洲 | 成人国产av品久久久| 97精品久久久久久久久久精品| 午夜免费观看性视频| 亚洲av日韩在线播放| 首页视频小说图片口味搜索 | 香蕉丝袜av| 一本色道久久久久久精品综合| 久久精品国产综合久久久| 男人舔女人的私密视频| 91精品三级在线观看| 亚洲av国产av综合av卡| 熟女少妇亚洲综合色aaa.| 人妻人人澡人人爽人人| 乱人伦中国视频| 亚洲av美国av| 这个男人来自地球电影免费观看| 精品一区在线观看国产| 丝袜脚勾引网站| 国产高清视频在线播放一区 | 少妇被粗大的猛进出69影院| 免费高清在线观看日韩| 免费不卡黄色视频| 国产爽快片一区二区三区| 国产在线免费精品| 久久久精品国产亚洲av高清涩受| 久久精品久久精品一区二区三区| 亚洲专区中文字幕在线| 视频区欧美日本亚洲| 免费不卡黄色视频| 成年人午夜在线观看视频| 国产精品香港三级国产av潘金莲 | 9191精品国产免费久久| 一个人免费看片子| 少妇裸体淫交视频免费看高清 | 欧美日韩亚洲国产一区二区在线观看 | 两性夫妻黄色片| 久久久久久亚洲精品国产蜜桃av| 精品卡一卡二卡四卡免费| 一个人免费看片子| 免费看十八禁软件| 我的亚洲天堂| 少妇粗大呻吟视频| 午夜福利在线免费观看网站| 大码成人一级视频| 午夜激情av网站| 亚洲国产精品999| 国产熟女午夜一区二区三区| netflix在线观看网站| 亚洲国产看品久久| 日韩大码丰满熟妇| 国产成人免费无遮挡视频| 少妇精品久久久久久久| 亚洲七黄色美女视频| 久久综合国产亚洲精品| 亚洲精品国产av蜜桃| 国产精品九九99| 日本猛色少妇xxxxx猛交久久| 亚洲欧美精品综合一区二区三区| 18禁国产床啪视频网站| 中文字幕人妻熟女乱码| 一级毛片我不卡| 国产午夜精品一二区理论片| 操出白浆在线播放| 一级毛片我不卡| 亚洲av美国av| 又大又黄又爽视频免费| 欧美人与性动交α欧美软件| 国产日韩欧美亚洲二区| 欧美激情高清一区二区三区| 美女扒开内裤让男人捅视频| 亚洲国产精品成人久久小说| 精品福利永久在线观看| 日韩av免费高清视频| 啦啦啦视频在线资源免费观看| 好男人电影高清在线观看| 最近中文字幕2019免费版| 免费不卡黄色视频| 久久精品国产综合久久久| 久久毛片免费看一区二区三区| 视频在线观看一区二区三区| 搡老乐熟女国产| 老汉色∧v一级毛片| 看免费成人av毛片| 亚洲欧洲国产日韩| 91精品国产国语对白视频| 国产在线一区二区三区精| 叶爱在线成人免费视频播放| 观看av在线不卡| 亚洲成人手机| 国产欧美日韩一区二区三区在线| 欧美黑人精品巨大| 在线观看免费视频网站a站| 一级片免费观看大全| 久久综合国产亚洲精品| 午夜精品国产一区二区电影| 90打野战视频偷拍视频| 中文精品一卡2卡3卡4更新| 欧美人与性动交α欧美软件| 国产精品久久久av美女十八| 国产高清不卡午夜福利| 高清av免费在线| 91精品国产国语对白视频| 极品少妇高潮喷水抽搐| 国产成人精品无人区| 黄色a级毛片大全视频| 晚上一个人看的免费电影| 国产一区二区激情短视频 | 最黄视频免费看| 99九九在线精品视频| kizo精华| 久久99精品国语久久久| 久久午夜综合久久蜜桃| 久久精品国产综合久久久| 精品亚洲乱码少妇综合久久| 大话2 男鬼变身卡| 啦啦啦视频在线资源免费观看| 两个人免费观看高清视频| 日本午夜av视频| 国产精品一区二区精品视频观看| 捣出白浆h1v1| 多毛熟女@视频| 另类精品久久| 亚洲精品美女久久av网站| 精品欧美一区二区三区在线| 99热全是精品| 美女高潮到喷水免费观看| 激情五月婷婷亚洲| 亚洲精品第二区| 美女高潮到喷水免费观看| 男女国产视频网站| 久久久久久免费高清国产稀缺| 久久国产精品大桥未久av| 久久国产精品人妻蜜桃| 中文字幕亚洲精品专区| 电影成人av| 人人妻人人添人人爽欧美一区卜| 国产亚洲精品第一综合不卡| 成人黄色视频免费在线看| 免费一级毛片在线播放高清视频 | 男人添女人高潮全过程视频| 国产色视频综合| 亚洲国产欧美网| av视频免费观看在线观看| 亚洲av国产av综合av卡| 在线观看免费高清a一片| 精品国产一区二区三区四区第35| 十八禁人妻一区二区| 丝袜美足系列| 日本91视频免费播放| 久久九九热精品免费| 中文乱码字字幕精品一区二区三区| 亚洲人成网站在线观看播放| 麻豆国产av国片精品| 亚洲av片天天在线观看| av视频免费观看在线观看| 99久久精品国产亚洲精品| 狂野欧美激情性bbbbbb| 免费高清在线观看日韩| 亚洲国产中文字幕在线视频| 18禁观看日本| 国产一区二区三区av在线| 啦啦啦在线免费观看视频4| 国产午夜精品一二区理论片| 亚洲av国产av综合av卡| 啦啦啦视频在线资源免费观看| 在线观看人妻少妇| 1024视频免费在线观看| 日本欧美国产在线视频| 99国产精品免费福利视频| av一本久久久久| 日本五十路高清| 国产精品国产三级专区第一集| 午夜福利影视在线免费观看| 一个人免费看片子| 91麻豆av在线| 水蜜桃什么品种好| 亚洲成色77777| 久久热在线av| 亚洲精品国产av成人精品| 国精品久久久久久国模美| 日本欧美视频一区| 成人黄色视频免费在线看| 丁香六月欧美| 国产一卡二卡三卡精品| 亚洲精品乱久久久久久| 国产免费又黄又爽又色| 黄色片一级片一级黄色片| 男人舔女人的私密视频| 别揉我奶头~嗯~啊~动态视频 | 99热网站在线观看| 水蜜桃什么品种好| 美女中出高潮动态图| 国产免费现黄频在线看| 国产精品亚洲av一区麻豆| 亚洲午夜精品一区,二区,三区| 国产日韩欧美在线精品| 成人国语在线视频| 国产爽快片一区二区三区| 男人操女人黄网站| 一边摸一边做爽爽视频免费| 亚洲精品自拍成人| 日韩 亚洲 欧美在线| 亚洲av日韩精品久久久久久密 | 丝袜在线中文字幕| 在线观看免费日韩欧美大片| 建设人人有责人人尽责人人享有的| 久久天躁狠狠躁夜夜2o2o | 亚洲欧美一区二区三区黑人| 国精品久久久久久国模美| 亚洲成人免费av在线播放| 亚洲激情五月婷婷啪啪| 男女之事视频高清在线观看 | av国产久精品久网站免费入址| 宅男免费午夜| 首页视频小说图片口味搜索 | 免费高清在线观看视频在线观看| 免费一级毛片在线播放高清视频 | 久9热在线精品视频| 一级黄片播放器| 亚洲第一青青草原| 亚洲国产av新网站| 国产欧美日韩一区二区三区在线| 免费看十八禁软件| 精品福利观看| 国产伦人伦偷精品视频| 精品国产国语对白av| 国产欧美亚洲国产| 九色亚洲精品在线播放| 久久国产精品影院| 热99久久久久精品小说推荐| 91国产中文字幕| 一边摸一边做爽爽视频免费| 国产精品国产三级专区第一集| 国产老妇伦熟女老妇高清| 久久久国产精品麻豆| 亚洲精品久久午夜乱码| 乱人伦中国视频| 新久久久久国产一级毛片| 香蕉国产在线看| 亚洲欧美日韩另类电影网站| 亚洲七黄色美女视频| 天天躁日日躁夜夜躁夜夜| 人人妻人人澡人人爽人人夜夜| 丁香六月欧美| 欧美变态另类bdsm刘玥| xxx大片免费视频| av在线app专区| 亚洲五月婷婷丁香| 热re99久久精品国产66热6| 电影成人av| 色婷婷久久久亚洲欧美| 欧美激情高清一区二区三区| 久久综合国产亚洲精品| 人人妻人人爽人人添夜夜欢视频| 91老司机精品| 在线亚洲精品国产二区图片欧美| 欧美亚洲 丝袜 人妻 在线| 免费在线观看视频国产中文字幕亚洲 | 久久精品久久精品一区二区三区| 亚洲成人国产一区在线观看 | 欧美+亚洲+日韩+国产| 国产成人精品久久二区二区91| 老司机亚洲免费影院| 男的添女的下面高潮视频| 99国产精品99久久久久| 女人久久www免费人成看片| 天堂8中文在线网| 亚洲视频免费观看视频| 高清av免费在线| 国产午夜精品一二区理论片| 高清av免费在线| 久久狼人影院| √禁漫天堂资源中文www| 午夜视频精品福利| 久久国产精品大桥未久av| 日日爽夜夜爽网站| 国产免费现黄频在线看| 咕卡用的链子| 视频区欧美日本亚洲| 国产精品国产三级专区第一集| 国产一区亚洲一区在线观看| 国产视频一区二区在线看| 丝袜人妻中文字幕| 下体分泌物呈黄色| 久久女婷五月综合色啪小说| 久久影院123| 国产主播在线观看一区二区 | 老汉色∧v一级毛片| 久久狼人影院| 最近最新中文字幕大全免费视频 | 亚洲成av片中文字幕在线观看| 夫妻午夜视频| 亚洲精品久久成人aⅴ小说| 久久午夜综合久久蜜桃| 美女中出高潮动态图| 国产亚洲午夜精品一区二区久久| 日韩制服丝袜自拍偷拍| 国产精品秋霞免费鲁丝片| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲精品国产av蜜桃| 精品熟女少妇八av免费久了| 老司机靠b影院| 女人精品久久久久毛片| 亚洲欧美色中文字幕在线| 亚洲精品乱久久久久久| 搡老岳熟女国产| 日本a在线网址| 精品少妇内射三级| 99国产综合亚洲精品| 亚洲午夜精品一区,二区,三区| 久久天躁狠狠躁夜夜2o2o | 精品久久久久久久毛片微露脸 | 自拍欧美九色日韩亚洲蝌蚪91| 97在线人人人人妻| 日韩制服骚丝袜av| 精品国产乱码久久久久久男人| 亚洲,一卡二卡三卡| 国产精品九九99| 亚洲 欧美一区二区三区| 久久精品国产亚洲av涩爱| 啦啦啦中文免费视频观看日本| 国产免费现黄频在线看| 王馨瑶露胸无遮挡在线观看| 19禁男女啪啪无遮挡网站| 国产免费视频播放在线视频| 亚洲,欧美精品.| 国产99久久九九免费精品| 午夜福利免费观看在线| 欧美日韩国产mv在线观看视频| 欧美激情高清一区二区三区| 国产精品一区二区精品视频观看| 久久性视频一级片| 老汉色∧v一级毛片| 国产精品秋霞免费鲁丝片| 男女午夜视频在线观看| 亚洲免费av在线视频| 精品亚洲成国产av| 国产精品亚洲av一区麻豆| 另类亚洲欧美激情| 男女国产视频网站| 亚洲,欧美,日韩| av国产精品久久久久影院| 九草在线视频观看| 丰满少妇做爰视频| 丝袜美腿诱惑在线| 2021少妇久久久久久久久久久| 51午夜福利影视在线观看|