• <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.免费观看| 在线播放国产精品三级| 国产高清不卡午夜福利| 草草在线视频免费看| 国产精品自产拍在线观看55亚洲| 中文字幕精品亚洲无线码一区| 蜜桃久久精品国产亚洲av| 精品99又大又爽又粗少妇毛片 | 男女之事视频高清在线观看| 最新在线观看一区二区三区| 在线播放无遮挡| 一进一出抽搐gif免费好疼| 一进一出好大好爽视频| 一边摸一边抽搐一进一小说| 日本一本二区三区精品| 精品一区二区三区人妻视频| 亚洲一级一片aⅴ在线观看| 嫩草影院入口| 午夜久久久久精精品| 国语自产精品视频在线第100页| 亚洲国产精品sss在线观看| 综合色av麻豆| 最新在线观看一区二区三区| 99久久九九国产精品国产免费| 亚洲乱码一区二区免费版| 亚洲av日韩精品久久久久久密| 精品一区二区三区视频在线| 成人综合一区亚洲| 嫩草影视91久久| 国产成人福利小说| 久久久久性生活片| 在线国产一区二区在线| 成人av一区二区三区在线看| a级毛片a级免费在线| 精品一区二区三区视频在线| 窝窝影院91人妻| 12—13女人毛片做爰片一| 深夜精品福利| 99久久精品国产国产毛片| 久久人妻av系列| 免费无遮挡裸体视频| 欧美bdsm另类| 久久亚洲真实| 九色国产91popny在线| 日本成人三级电影网站| 美女黄网站色视频| 久久午夜亚洲精品久久| 日韩大尺度精品在线看网址| 91午夜精品亚洲一区二区三区 | 两性午夜刺激爽爽歪歪视频在线观看| 在线观看一区二区三区| 毛片一级片免费看久久久久 | 自拍偷自拍亚洲精品老妇| 噜噜噜噜噜久久久久久91| 可以在线观看的亚洲视频| 国产av麻豆久久久久久久| 国产精品综合久久久久久久免费| 亚洲最大成人av| 一级黄片播放器| 亚洲中文日韩欧美视频| 国产 一区精品| 日本免费a在线| 久久久国产成人精品二区| 欧美日韩国产亚洲二区| 日本 av在线| 毛片一级片免费看久久久久 | 亚洲精品成人久久久久久| 一个人免费在线观看电影| 又爽又黄无遮挡网站| 久久久久久久久久成人| 久久这里只有精品中国| 毛片女人毛片| 看黄色毛片网站| 日韩欧美国产一区二区入口| 少妇猛男粗大的猛烈进出视频 | 久久草成人影院| 久久精品国产亚洲av涩爱 | 97碰自拍视频| 两个人视频免费观看高清| 亚洲美女搞黄在线观看 | 高清日韩中文字幕在线| 1024手机看黄色片| 国产亚洲精品综合一区在线观看| 欧美日韩精品成人综合77777| 欧美成人性av电影在线观看| 男女那种视频在线观看| 97碰自拍视频| 少妇猛男粗大的猛烈进出视频 | 女的被弄到高潮叫床怎么办 | av女优亚洲男人天堂| 啦啦啦韩国在线观看视频| 国产精品综合久久久久久久免费| 美女xxoo啪啪120秒动态图| 给我免费播放毛片高清在线观看| 日本一本二区三区精品| 亚洲av中文av极速乱 | 亚洲专区国产一区二区| 国产国拍精品亚洲av在线观看| 午夜精品久久久久久毛片777| av专区在线播放| 嫩草影视91久久| 亚洲av一区综合| 婷婷丁香在线五月| 亚洲中文字幕一区二区三区有码在线看| 精品久久久久久久久av| 一级av片app| 欧美精品啪啪一区二区三区| 欧美日本亚洲视频在线播放| 午夜老司机福利剧场| 九色国产91popny在线| 91久久精品电影网| 我要搜黄色片| 久久久久精品国产欧美久久久| 国产熟女欧美一区二区| 亚洲av美国av| 国产日本99.免费观看| 亚洲欧美日韩东京热| 久久精品夜夜夜夜夜久久蜜豆| 一边摸一边抽搐一进一小说| 一卡2卡三卡四卡精品乱码亚洲| a级一级毛片免费在线观看| 亚洲成人免费电影在线观看| 美女免费视频网站| 国产精品一区二区三区四区久久| 欧美成人性av电影在线观看| 日韩高清综合在线| 欧美色视频一区免费| 亚洲人成伊人成综合网2020| 欧美成人免费av一区二区三区| 真实男女啪啪啪动态图| 人人妻人人看人人澡| 在线观看66精品国产| 亚洲精品在线观看二区| 一区福利在线观看| 欧美一区二区国产精品久久精品| 欧美xxxx黑人xx丫x性爽| 国产成人av教育| 中文字幕高清在线视频| 久久精品夜夜夜夜夜久久蜜豆| 男插女下体视频免费在线播放| 99精品在免费线老司机午夜| av国产免费在线观看| 最近最新免费中文字幕在线| 欧美最黄视频在线播放免费| 欧美日韩精品成人综合77777| www.www免费av| 精品久久久噜噜| 日日啪夜夜撸| 亚洲精品久久国产高清桃花| 色5月婷婷丁香| avwww免费| 日本一二三区视频观看| 欧美人与善性xxx| 草草在线视频免费看| 欧美色视频一区免费| 在线播放无遮挡| 国内揄拍国产精品人妻在线| 亚洲人与动物交配视频| av福利片在线观看| 超碰av人人做人人爽久久| 精品久久久久久成人av| 欧美一区二区精品小视频在线| 乱人视频在线观看| 亚洲不卡免费看| 日韩欧美 国产精品| 国产伦人伦偷精品视频| 成人综合一区亚洲| 搡老岳熟女国产| 日韩高清综合在线| 欧美绝顶高潮抽搐喷水| 欧美黑人欧美精品刺激| 精品午夜福利视频在线观看一区| 成人av一区二区三区在线看| 亚洲av第一区精品v没综合| av黄色大香蕉| 亚洲成人精品中文字幕电影| 一边摸一边抽搐一进一小说| 国产伦在线观看视频一区| 久久精品国产亚洲av香蕉五月| 成人av在线播放网站| 内射极品少妇av片p| 免费电影在线观看免费观看| 成人二区视频| 永久网站在线| 亚洲黑人精品在线| 尾随美女入室| 国产男人的电影天堂91| 99久国产av精品| 日本黄色视频三级网站网址| 一夜夜www| 精品久久久久久久久久免费视频| av天堂在线播放| 男女视频在线观看网站免费| 成人欧美大片| 精品久久久久久久久亚洲 | 99国产极品粉嫩在线观看| 丝袜美腿在线中文| 国产伦精品一区二区三区四那| 精品久久久久久,| 国内毛片毛片毛片毛片毛片| 国产精品久久久久久久久免| 国产毛片a区久久久久| 成人国产一区最新在线观看| 99热网站在线观看| 日日干狠狠操夜夜爽| 99热6这里只有精品| 成年免费大片在线观看| 联通29元200g的流量卡| 婷婷亚洲欧美| 亚洲精品日韩av片在线观看| 欧美性猛交黑人性爽| а√天堂www在线а√下载| 热99re8久久精品国产| 啦啦啦观看免费观看视频高清| 婷婷精品国产亚洲av| 婷婷六月久久综合丁香| 国产日本99.免费观看| 99视频精品全部免费 在线| 综合色av麻豆| 久久精品国产亚洲av天美| 日韩精品青青久久久久久| 一个人观看的视频www高清免费观看| 国产伦精品一区二区三区四那| 热99在线观看视频| 久久99热6这里只有精品| 日韩中文字幕欧美一区二区| 免费看美女性在线毛片视频| 美女被艹到高潮喷水动态| 一区福利在线观看| 免费看日本二区| 国产精品久久久久久久电影| 一级毛片久久久久久久久女| 小蜜桃在线观看免费完整版高清| 国产亚洲91精品色在线| 99热这里只有精品一区| 人人妻人人澡欧美一区二区| 国产人妻一区二区三区在| 国产乱人伦免费视频| 国内久久婷婷六月综合欲色啪| 精品一区二区三区人妻视频| 18禁裸乳无遮挡免费网站照片| 欧美+日韩+精品| 波多野结衣高清作品| 校园春色视频在线观看| av在线观看视频网站免费| 免费观看的影片在线观看| 老司机深夜福利视频在线观看| 天堂动漫精品| 国产又黄又爽又无遮挡在线| 一区福利在线观看| 精品午夜福利视频在线观看一区| 三级男女做爰猛烈吃奶摸视频| 久久中文看片网| 国产欧美日韩一区二区精品| 国产伦精品一区二区三区四那| 久久人人精品亚洲av| 国产伦人伦偷精品视频| 极品教师在线视频| 人人妻人人看人人澡| 人人妻人人澡欧美一区二区| 亚洲精品粉嫩美女一区| 成年人黄色毛片网站| 国产aⅴ精品一区二区三区波| 熟妇人妻久久中文字幕3abv| 人人妻,人人澡人人爽秒播| 久久久精品欧美日韩精品| 麻豆av噜噜一区二区三区| 黄色欧美视频在线观看| 成年女人看的毛片在线观看| 国产精品av视频在线免费观看| 乱人视频在线观看| 国产视频内射| 国产精品一区二区三区四区免费观看 | 日本一本二区三区精品| 成年版毛片免费区| 91av网一区二区| 天天躁日日操中文字幕| 美女高潮的动态| 色5月婷婷丁香| 欧美日本亚洲视频在线播放| 午夜亚洲福利在线播放| 美女 人体艺术 gogo| 尾随美女入室| av在线观看视频网站免费| 偷拍熟女少妇极品色| 欧美极品一区二区三区四区| .国产精品久久| 91久久精品国产一区二区三区| 精品久久久久久久久久免费视频| 婷婷六月久久综合丁香| 亚洲国产精品久久男人天堂| 午夜福利欧美成人| 免费看日本二区| 婷婷精品国产亚洲av在线| 久久久国产成人免费| 又紧又爽又黄一区二区| 97热精品久久久久久| 日韩精品中文字幕看吧| 亚洲av第一区精品v没综合| 在线播放无遮挡| 亚洲精品粉嫩美女一区| 日韩一区二区视频免费看| 深夜精品福利| 成人美女网站在线观看视频| 国产精品人妻久久久影院| 国内久久婷婷六月综合欲色啪| 精华霜和精华液先用哪个| 非洲黑人性xxxx精品又粗又长| 亚洲中文字幕日韩| 免费大片18禁| 99久久成人亚洲精品观看| 人妻夜夜爽99麻豆av| 国产精品爽爽va在线观看网站| 色视频www国产| 伦精品一区二区三区| 色综合亚洲欧美另类图片| 毛片女人毛片| 亚洲中文日韩欧美视频| 国产三级在线视频| www.www免费av| 精品国内亚洲2022精品成人| 亚洲成人久久性| 俄罗斯特黄特色一大片| 蜜桃亚洲精品一区二区三区| 免费观看在线日韩| 日韩强制内射视频| 成熟少妇高潮喷水视频| 12—13女人毛片做爰片一| 一夜夜www| 亚洲无线观看免费| 免费大片18禁| 亚洲va日本ⅴa欧美va伊人久久| 国产精品女同一区二区软件 | 成人午夜高清在线视频| 伦理电影大哥的女人| 少妇裸体淫交视频免费看高清| 亚洲欧美日韩高清在线视频| 国产 一区精品| 亚洲最大成人手机在线| 免费av观看视频| 亚洲精品成人久久久久久| 高清毛片免费观看视频网站| 九色国产91popny在线| 我的老师免费观看完整版| 日本免费a在线| 国内久久婷婷六月综合欲色啪| 久久精品91蜜桃| 中文字幕熟女人妻在线| 中文字幕免费在线视频6| 看十八女毛片水多多多| 美女黄网站色视频| 久久久久九九精品影院| 久久久精品大字幕| 此物有八面人人有两片| 嫩草影院入口| 九九在线视频观看精品| a在线观看视频网站| 亚洲人成网站在线播放欧美日韩| 国产成人影院久久av| 国产亚洲精品久久久久久毛片| 性插视频无遮挡在线免费观看| 精品福利观看| 精品久久久久久久末码| 网址你懂的国产日韩在线| 亚洲欧美日韩卡通动漫| 亚洲专区国产一区二区| 免费观看的影片在线观看| 成年女人看的毛片在线观看| 在线免费十八禁| 久久热精品热| 两人在一起打扑克的视频| 此物有八面人人有两片| 欧美中文日本在线观看视频| 成人三级黄色视频| 内地一区二区视频在线| 国产私拍福利视频在线观看| 啦啦啦观看免费观看视频高清| 久久精品影院6| 国产av麻豆久久久久久久| 琪琪午夜伦伦电影理论片6080| 免费一级毛片在线播放高清视频| 狂野欧美激情性xxxx在线观看| 欧美高清成人免费视频www| 听说在线观看完整版免费高清| 国产精品人妻久久久久久| 国产精品精品国产色婷婷| 波多野结衣巨乳人妻| 狂野欧美激情性xxxx在线观看| 在线播放无遮挡| 亚洲av一区综合| 亚洲欧美日韩东京热| 欧美在线一区亚洲| 中文亚洲av片在线观看爽| 亚洲最大成人手机在线| 日韩一本色道免费dvd| 一级黄片播放器| 亚洲乱码一区二区免费版| 91精品国产九色| 亚洲三级黄色毛片| АⅤ资源中文在线天堂| 国产欧美日韩精品一区二区| 99九九线精品视频在线观看视频| 又黄又爽又免费观看的视频| 亚洲精华国产精华精| 国产高清不卡午夜福利| 99在线视频只有这里精品首页| 久久热精品热| 性欧美人与动物交配| 色噜噜av男人的天堂激情| 国内久久婷婷六月综合欲色啪| 国产精品,欧美在线| 十八禁国产超污无遮挡网站| 国产一区二区激情短视频| 在线播放无遮挡| 最近中文字幕高清免费大全6 | 欧美性感艳星| 老熟妇乱子伦视频在线观看| 俺也久久电影网| 免费在线观看影片大全网站| 在线免费观看不下载黄p国产 | 久久精品国产亚洲av天美| 国产私拍福利视频在线观看| 欧美潮喷喷水| 欧美一区二区亚洲| h日本视频在线播放| 春色校园在线视频观看| 日韩人妻高清精品专区| 蜜桃亚洲精品一区二区三区| 亚洲av.av天堂| 大又大粗又爽又黄少妇毛片口| 亚洲精品粉嫩美女一区| 国产主播在线观看一区二区| 男女视频在线观看网站免费| 日日夜夜操网爽| 色噜噜av男人的天堂激情| 韩国av在线不卡| 国内精品美女久久久久久| 亚洲av熟女| 中文资源天堂在线| 国产男人的电影天堂91| 高清在线国产一区| 亚洲乱码一区二区免费版| 简卡轻食公司| 国产精品亚洲一级av第二区| 国产伦人伦偷精品视频| 午夜日韩欧美国产| 深爱激情五月婷婷| 波野结衣二区三区在线| 国产真实乱freesex| 欧美日本亚洲视频在线播放| 欧美国产日韩亚洲一区| 久久人人爽人人爽人人片va| 欧美日韩国产亚洲二区| 中国美女看黄片| 我的老师免费观看完整版| 在线看三级毛片| 18禁黄网站禁片免费观看直播| 精品久久久久久久久久久久久| 99在线人妻在线中文字幕| 亚洲无线观看免费| videossex国产| 亚洲一区二区三区色噜噜| 精品人妻熟女av久视频| 婷婷色综合大香蕉| 精品久久久久久久久久免费视频| 综合色av麻豆| 看十八女毛片水多多多| 在线免费十八禁| 无遮挡黄片免费观看| 亚洲第一电影网av| 国产精品日韩av在线免费观看| 中国美白少妇内射xxxbb| 免费人成视频x8x8入口观看| 国产午夜福利久久久久久| 韩国av在线不卡| 热99在线观看视频| 国产欧美日韩精品一区二区| 欧美成人性av电影在线观看| 欧美日韩黄片免| 一级黄片播放器| 久久精品夜夜夜夜夜久久蜜豆| 久久久精品大字幕| 免费观看人在逋| 九九在线视频观看精品| 亚洲人成伊人成综合网2020| 九九久久精品国产亚洲av麻豆| 啦啦啦韩国在线观看视频| 国产黄色小视频在线观看| 亚洲成人免费电影在线观看| 99热这里只有是精品在线观看| 婷婷六月久久综合丁香| 亚洲av中文字字幕乱码综合| 欧美成人性av电影在线观看| 人妻夜夜爽99麻豆av| 成人二区视频| 亚洲18禁久久av| 国产伦一二天堂av在线观看| 嫁个100分男人电影在线观看| 嫩草影视91久久| 国产精品久久电影中文字幕| 精品久久国产蜜桃| 欧美不卡视频在线免费观看| 国产精品嫩草影院av在线观看 | 国产精品美女特级片免费视频播放器| 国产一区二区在线观看日韩| 精品国内亚洲2022精品成人| 成人午夜高清在线视频| 99久久久亚洲精品蜜臀av| 国产v大片淫在线免费观看| 深夜精品福利| 不卡视频在线观看欧美| 男插女下体视频免费在线播放| 国产av在哪里看| av黄色大香蕉| 欧美性猛交╳xxx乱大交人| 女同久久另类99精品国产91| 97超级碰碰碰精品色视频在线观看| 免费av毛片视频| 精品久久久久久久人妻蜜臀av| 99久久无色码亚洲精品果冻| 久久精品人妻少妇| 九九爱精品视频在线观看| 成年女人永久免费观看视频| 国产毛片a区久久久久| 在线免费十八禁| 一a级毛片在线观看| 国产精品电影一区二区三区| 色播亚洲综合网| 日本免费一区二区三区高清不卡| 欧美成人免费av一区二区三区| АⅤ资源中文在线天堂| 免费看光身美女| 最近在线观看免费完整版| 亚洲欧美日韩无卡精品| 日韩大尺度精品在线看网址| 欧美zozozo另类| 精品人妻偷拍中文字幕| 欧美国产日韩亚洲一区| 国产一区二区三区在线臀色熟女| 九九在线视频观看精品| 成人鲁丝片一二三区免费| 99热只有精品国产| 伦精品一区二区三区| 在线观看午夜福利视频| 日本黄色视频三级网站网址| 能在线免费观看的黄片| 国产av在哪里看| 成年女人毛片免费观看观看9| 美女大奶头视频| 久久99热6这里只有精品| 国产乱人视频| 午夜激情欧美在线| 亚洲精品亚洲一区二区| 亚洲欧美激情综合另类| 免费看a级黄色片| 亚洲欧美激情综合另类| 亚洲成人精品中文字幕电影| 三级毛片av免费| 国产三级在线视频| 性欧美人与动物交配| 久久久午夜欧美精品| 亚洲人成伊人成综合网2020| 日日啪夜夜撸| 桃色一区二区三区在线观看| 久久人人精品亚洲av| 极品教师在线视频| 国产精品av视频在线免费观看| 亚洲美女搞黄在线观看 | 长腿黑丝高跟| 色综合婷婷激情| 日韩欧美精品免费久久| 精品人妻视频免费看| 欧美一级a爱片免费观看看| 欧美在线一区亚洲| 久久久午夜欧美精品| 国产高清三级在线| 欧美激情在线99| 此物有八面人人有两片| 国产色爽女视频免费观看| 一区二区三区四区激情视频 | 国产白丝娇喘喷水9色精品| 久久国内精品自在自线图片| 午夜福利高清视频| 老司机午夜福利在线观看视频| 黄色女人牲交| 淫妇啪啪啪对白视频| 91久久精品国产一区二区三区| 免费看光身美女| 一边摸一边抽搐一进一小说| 亚洲一区二区三区色噜噜| 亚洲人成伊人成综合网2020| 五月伊人婷婷丁香| 老女人水多毛片| 亚洲无线在线观看| 亚洲国产精品合色在线| 男插女下体视频免费在线播放| 欧美+亚洲+日韩+国产| 女的被弄到高潮叫床怎么办 | 观看免费一级毛片| 国产国拍精品亚洲av在线观看| 国产午夜精品论理片| 午夜福利在线观看吧| 成人毛片a级毛片在线播放| 又爽又黄无遮挡网站| 香蕉av资源在线| 午夜视频国产福利| 天美传媒精品一区二区|