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

    頂煤放落過程煤矸聲信號特征提取與分類方法

    2021-11-06 05:54:16袁源汪嘉文朱德昇王家臣王統(tǒng)海楊克虎
    礦業(yè)科學學報 2021年6期
    關鍵詞:煤矸放頂頻域

    袁源汪嘉文朱德昇王家臣王統(tǒng)海楊克虎

    1.中國礦業(yè)大學(北京)機電與信息工程學院,北京 100083;2.中國礦業(yè)大學(北京)能源與礦業(yè)學院,北京 100083;3.山東能源臨沂礦業(yè)集團有限責任公司,山東臨沂 276022

    推進煤礦智能開采是我國能源技術的發(fā)展趨勢,是信息革命的時代要求。 相比其他采煤法,放頂煤采煤法具有高產高效、巷道掘進率低等優(yōu)勢而被廣泛應用。 優(yōu)化改進放煤工藝是實現煤礦智能開采的一個重要環(huán)節(jié)[1-2]。 傳統(tǒng)放頂煤開采過程環(huán)境感知困難,煤礦工人工作環(huán)境惡劣、粉塵多、視野窄,加之放煤工藝煩瑣,容易造成過放、欠放的現象,無法保證資源的回采率與煤質的優(yōu)良率。 煤矸識別方法的引入可以讓煤礦工人更好地掌控放煤口的開關時機,為煤礦智能開采提供技術支持。

    學術界對煤矸識別方法進行了大量的研究。針對放落煤和矸石的紋理形態(tài)、灰度特征,文獻[3-5]研究了基于圖像的煤矸識別方法,該方法因樣本集易受光照、粉塵干擾,識別精度受到了一定程度的影響;針對粒子波在煤層與矸石層傳播的差異性,文獻[6]研究了基于電磁波的煤矸識別方法,該方法適用于頂煤探測,由于存在波長衰減,探測范圍和煤矸識別精度被頂煤厚度所制約;針對煤與矸石所包含元素放射性強度的不同,文獻[7]研究了基于γ射線的煤矸識別方法,適用于頂煤探測,由于受到矸石占比變化的影響,回采率較低,且存在放射性危害;針對采煤機割煤時滾筒扭矩的變化,文獻[8]研究了基于扭矩信號的煤矸識別方法,該方法適用于采煤機割煤過程的煤矸識別,但識別裝置對傳感器要求較高;針對煤矸普氏系數的差異,文獻[9]研究了基于振動信號的煤矸識別方法,該方法識別效果受到了工作面復雜程度和機械振動的影響,抗干擾能力較差;針對煤和矸石放落到液壓支架尾梁和后部刮板輸送機上的聲信號物理特性差別,文獻[10-11]研究了基于聲信號的煤矸識別方法,該方法抗干擾能力較強,但是目前采集與識別系統(tǒng)尚未完善。

    基于聲信號的煤矸識別方法,因其成本低、標注便捷而廣受重視。 文獻[12]采用了獨立分量分析法對頂煤放落過程的聲信號進行分析處理,根據頻譜差異初步判斷了頂煤中的含矸率;文獻[13]考慮了環(huán)境噪聲,采用含噪的超完備獨立分量分析方法解決放頂煤過程中聲信號分離的問題;文獻[14]采用小波包變換的分析方法,發(fā)現頻譜會隨煤矸混合比例不同而產生差異;文獻[15]經實驗分析得出煤矸聲波時域信號中的峰峰值、方差和峭度系數區(qū)分度明顯;文獻[16]通過獨立分量分析將放煤聲信號分離,并提取其梅爾倒譜系數特征作為神經網絡的輸入,進行煤矸識別;文獻[17]通過希爾波特-黃變換的方法提取聲振信號特征,基于雙峰深度學習網絡實現煤矸分類。

    目前,放頂煤聲信號特征提取和分類算法的研究對比較少,且大部分測試沒有現場數據的支持。本文完成了工作面放頂煤聲信號的采集,通過人工標注建立放頂煤聲信號樣本庫來研究聲信號的機理,綜合分析樣本幀長、特征差異、算法性能,對機器學習分類模型進行評估。

    1 放頂煤聲信號特征提取方法

    放頂煤聲信號的特征可以用音調、響度、音色來描述。 聲信號的高低稱為音調,聲信號振動的頻率越高,即聲源物體振動得越快,音調越高;聲信號的強弱稱為響度,聲信號的振幅越大,響度越大;聲信號的波形特征表現稱為音色,聲信號波形的各個諧波頻譜與包絡決定了音色的特性。 本文通過濾波去噪,分析振幅、相位與時間、頻率的對應關系來提取聲信號的特征。

    1.1 時域特征提取

    放頂煤聲信號的本質為非平穩(wěn)信號[18],可使用非平穩(wěn)信號的時域分析方法進行特征提取,放頂煤聲信號的時域特征及計算公式見表1。

    表1 放頂煤聲信號的時域特征及其計算公式Tab.1 Time-domain characteristics of caving coal acoustic signal and its calculation formula

    表1 選取了有效值、標準差、偏度、峭度、峰峰值、波形因數、脈沖因數和峰值因數8 種特征。 其中,x(n)為放頂煤聲信號;N為序列長度;Xrms為信號的有效值;Xstd為信號的標準差;Xske為信號的偏度;Xkur為信號的峭度;Xp-p為信號的峰峰值;CF為信號的波形因數;SF為信號的脈沖因數;IF為信號的峰值因數。

    1.2 頻域特征提取

    放頂煤聲信號是數字采樣信號,可使用離散傅里葉變換(Discrete Fourier Transform, DFT)方法進行頻域特征提取。 離散傅里葉變換公式為

    式中,X(k)為頻率系數;WN為蝶形因子。

    基于蝶形因子的周期性和對稱性,在DFT 基礎上改進,得到快速傅里葉變換(Fast Fourier Transform, FFT)方法,將復數乘法運算轉換為復數加法運算,大大降低了復雜程度與運算量,非常適用于計算機運算。 進一步計算得到聲信號的頻域特征及計算公式見表2。 其中,XFC為信號的重心頻率;XMSF為信號的均方頻率;XVF為信號的頻率方差。

    表2 放頂煤聲信號的頻域特征及其計算公式Tab.2 Frequency-domain characteristics of caving coal acoustic signal and its calculation formula

    1.3 時頻域特征提取

    放頂煤聲信號局部特征區(qū)別較大,可使用小波變換(Wavelet Transform, WT)[19]的方法進行時頻域特征提取。 小波變換可以完成多尺度化的時頻分析,但是其分辨率隨著頻率增加而降低,高頻處難以區(qū)分。 為對其適用性進行改進,衍生出小波包分解(Wavelet Packet Decomposition, WPD)[20]的方法,小波包分解引入了不同中心頻率的帶通濾波器,以獲得各個節(jié)點的信號分量。

    各頻帶分量的能量:

    式中,E(i,j)為第j層上的第i個小波包頻帶能量。

    對各頻帶能量歸一化后,得到各頻帶能量占比:

    2 機器學習分類算法的理論基礎

    執(zhí)行機器學習算法來完成對放頂煤聲信號樣本集數據的學習過程,獲得分類模型。 該模型本質上是對聲信號數據集的劃分。

    2.1 支持向量機算法

    支持向量機 ( Support Vector Machine,SVM)[21]算法的目標是找到一個最優(yōu)的分割面,使得樣本分類間隔最大。 在樣本空間中,劃分超平面為ωTx+b=0,線性支持向量機的目標函數與約束條件為

    式中,ω為超平面法向量;b為偏移項;(xi,yi)為樣本坐標點;γ為樣本點到超平面的距離,通常取值為1。

    支持向量機算法具有良好的學習能力,泛化能力強,可以將多分類問題以二分類的方法解決,但它對參數選擇要求較高。

    2.2 K 近鄰算法

    K 近鄰算法(K-Nearest-Neighbor, KNN)[22]無須使用訓練集進行訓練,只需計算測試樣本到訓練樣本的距離,選取距離測試集最近的k個訓練樣本,通過加權投票法得到測試樣本的類別。

    常用的距離衡量標準有LP 距離、歐氏距離、切比雪夫距離和曼哈頓距離。 其中歐氏距離為

    式中,(xj,yj)為空間坐標點;d為空間維數。

    K 近鄰算法理解簡單、易于實現,可用于分類、回歸問題,對異常值不敏感,但計算太復雜,屬于懶惰學習,需要大量的存儲空間。

    2.3 樸素貝葉斯算法

    樸素貝葉斯算法(Naive Bayes Classifier,NBC)[23]是在貝葉斯定理的基礎上增加了條件特征獨立性假設的算法。

    樸素貝葉斯的推斷標準為

    式中,c為類別;M為屬性數目;xm為x在第m個屬性特征上的取值。

    在實際分類問題中很難保證各項特征條件相互獨立,樸素貝葉斯模型精度受到了特征相關性的影響,但算法本身比較簡單,常用于多分類問題。

    2.4 決策樹算法

    決策樹算法(Decision Tress, DT)[24]應用了經典的樹形結構。 ID3(Iterative Dichotomiser 3)算法以信息增益作為特征。 C4.5 算法是ID3 的改進算法,支持對連續(xù)值輸入樣本節(jié)點的劃分,滿足對缺失值的處理要求。 CART(Classification And Regression Tree)算法為二叉樹結構。

    信息增益可表示為

    式中,Ent(D)為信息熵;D為樣本集合;a為離散屬性。

    有V個可能的取值,其中第v個分支節(jié)點包含了D中所有在屬性a上取值為av的樣本,記作Dv。

    決策樹算法易于實現,可在短時間內對大量數據做出可行性較高的判斷,但解決多分類問題的效果不理想,且難以預測連續(xù)值。

    2.5 隨機森林算法

    隨機森林算法(Radom Forests, RF)[25]是決策樹的集合,它的核心思想是增加樹的多樣性。

    RF 模型包含z輪訓練,即z個決策樹模型,其分類序列可表示為[DT1(ζ),DT2(ζ),…,DTz(ζ)],再進行投票表決。 分類決策為

    式中,M(ζ)為RF 分類模型;DTz(ζ)為單個分類模型;I為隨機變量的信息;T為輸出變量。

    隨機森林算法常用于多分類任務,具備較高的抗過擬合能力,對訓練集樣本數量與特征數量有著較高要求。

    2.6 神經網絡算法

    神經網絡算法[26],即人工神經網絡(Artificial Neural Network, ANN)算法。 常用的模型為多層感知器(Multi-Layer Perceptron, MLP)。 神經網絡算法模型的復雜程度隨著參數的增加而增加,復雜的網絡結構可以進一步完成更加復雜的學習任務,但該算法存在訓練效率低、易于過擬合等問題。

    3 放頂煤聲信號樣本集構建

    放頂煤聲信號樣本數據主要分為實驗模擬采集數據與工作面現場采集數據。 實驗室未考慮綜放工作面惡劣環(huán)境與復雜工況,模擬論證階段的數據較為理想;工作面現場采集的數據反映了真實的放煤過程,但數據量較少,且大都是人工參與采集,采集裝置成本較高、適用性較差。 缺乏大量現場數據作為支撐,使得現有的煤矸識別方法難以投入實際生產中。

    針對這一現象,本文設計了一種低成本、無須人工干預、可長時間采集、使用方便、適用性高的放頂煤液壓支架尾梁聲振數據采集裝置。 該裝置安裝于液壓支架尾梁內側,采用內部電池供電,覆蓋防水透氣膜并置于防水外殼中,通過對比加速度傳感器信號變化,結合運動識別算法,判斷液壓支架是否為放煤狀態(tài)。 若是,則進一步采集放頂煤聲振數據,并存儲于裝置中。

    數據采集裝置電路結構如圖1 所示,控制模塊選取了低功耗微控制器STM32L051 和高性能微控制器STM32F407,麥克風為MEMS 麥克風ICS-43432,振動傳模塊選取了MEMS 加速度計ADXL1002,存儲器容量為16 GB,采樣頻率為48 kHz,采樣深度為24 bit。

    圖1 數據采集裝置電路結構Fig.1 Circuit structure diagram of data acquisition device

    本實驗數據均采集自山東能源古城煤礦3106綜放工作面現場,數據采集裝置實物如圖2 所示。將采集到的音頻文件進行人工標注,劃分為不同類型的聲信號。 圖3 為放煤聲信號,圖4 為放矸聲信號。

    圖2 放頂煤液壓支架尾梁聲振數據采集裝置Fig.2 Acoustic and vibration data acquisition device for tail beam of caving coal hydraulic support

    圖3 不同噪聲影響下的放煤聲信號Fig.3 Coal sound signal under the influence of different noises

    圖4 不同噪聲影響下的放矸聲信號Fig.4 Gangue sound signal under the influence of different noises

    放頂煤聲信號樣本庫共有1 012.97 MB,總時長為5 825 s。 其中,放煤聲信號樣本共912.64 MB,總時長為5 296 s;放矸聲信號樣本共100.33 MB,總時長為529 s。

    放頂煤聲信號時域特征提取,是采用了非平穩(wěn)信號時域分析的方法,獲得有效值、方差、偏度、峭度、峰峰值、波形因數、脈沖因數、峰值因數作為輸入的8 維特征向量。 頻域特征提取是通過FFT 算法對原始信號進行運算,得到重心頻率、均方頻率、頻率方差作為輸入的3 維特征向量。 時頻域特征提取是采用小波包分解的方法,設置小波基為sym8、最大分解層數為6,根據能量守恒定律計算出每個頻帶分量的能量占比,將其作為輸入的64維特征向量。 放煤聲信號和放矸聲信號的小波包分解頻帶能量占比如圖5、圖6 所示。

    圖5 放煤聲信號小波包分解頻帶能量占比Fig.5 Frequency band energy ratio of wavelet packet decomposition of coal caving acoustic signal

    圖6 放矸聲信號小波包分解頻帶能量占比Fig.6 Frequency band energy ratio of wavelet packet decomposition of gangue acoustic signal

    4 機器學習分類模型訓練

    本實驗環(huán)境為Windows10 系統(tǒng),CPU 為Intel Core i7-8750H 2.20 GHz,運行內存為8 G,數據處理、模型訓練與方法驗證均在PyCharm 上運行。

    支持向量機、隨機森林、K 近鄰算法、樸素貝葉斯、決策樹和神經網絡均使用預處理后的聲信號樣本進行訓練,輸入特征分別為時域特征、頻域特征和時頻域特征。 考慮RF 模型性能與系統(tǒng)運行效率,設定隨機森林子樹數量n_estimators=100,通過oob_score 屬性打印模型的oob 分數,設定oob_score = True。 神經網絡模型采用擬牛頓法優(yōu)化器,設定lbfgs= quasi-Newton,正則化參數alpha 設置為0.000 01,隱藏層參數設置為(5,5)。 其余模型參數均為默認。

    使用交叉驗證的方法對煤矸識別準確率進行評估,采用k-fold 交叉驗證法,設置參數cv=10。分類模型的參數見表3。 每一種模型的識別時間均在3 ms 左右,滿足工作面的實時性要求。

    表3 多種分類模型的參數Tab.3 Parameters of multiple classification models

    5 分類模型應用與性能評估

    基于放頂煤聲信號樣本集,運用不同的特征提取方法與分類算法完成模型分類應用實驗與性能評估。

    5.1 樣本幀長對分類準確率的影響

    不同幀長的樣本會影響煤矸識別的準確率,本文選取幀長為50 ms、100 ms、150 ms、200 ms、250 ms、300 ms 的6 種情況,幀移默認為1/2 的幀長,針對不同特征與不同分類器進行對比實驗。

    對于時域特征,SVM、MLP 的模型分類準確率隨著幀長的增加而降低,如圖7 所示;RF、KNN、DT的模型分類準確率隨著幀長增加而增加,在200 ms 之后趨于平緩,準確率基本保持不變;NBC 的模型分類準確率波動較大,在100 ms 時取得最大值。 RF、DT 模型表現較好,RF 模型在幀長為200 ms 時,最大值為91.60% 。

    圖7 時域特征提取下不同幀長的煤矸識別準確率Fig.7 Accuracy of coal gangue recognition in different frame lengths under time-domain feature extraction

    對于頻域特征,RF 模型的分類準確率隨著幀長的增加而緩慢增加(圖8),其余5 種算法的準確率均有較大波動。 NBC 模型表現較差;KNN、RF、MLP 模型表現較好;KNN 模型在幀長為50 ms 時,煤矸分類準確率為87.69% ,隨著幀長的增加準確率有所下降, 在幀長為150 ms 時最大值為87.61% ,較為接近幀長50 ms 的表現。

    圖8 頻域特征提取下不同幀長的煤矸識別準確率Fig.8 Accuracy of coal gangue recognition in different frame lengths under frequency-domain feature extraction

    對于時頻域特征,隨著幀長的增加,SVM 模型分類準確率持續(xù)下降(圖9),其余5 種模型均有上升的趨勢。 RF、KNN、DT、MLP 模型表現較好;RF模型在幀長150 ms 時煤矸識別準確率為93.05% ,之后增長趨于平緩,受限于樣本數量,其準確率甚至有所下降,200 ms 幀長的煤矸識別準確率,最大值為93.35% 。

    圖9 時頻域特征提取下不同幀長的煤矸識別準確率Fig.9 Accuracy of coal gangue recognition in different frame lengths under timedomain feature extraction

    綜合考慮,本次實驗最佳幀長為200 ms,幀移為100 ms。 6 種機器學習算法模型中,SVM、NBC模型分類表現較差,之后的實驗主要是對其余4 種分類模型進行分析。

    5.2 特征選取對于分類準確率的影響

    不同的特征選取會影響到煤矸識別的準確率,本文研究不同特征向量對煤矸識別準確率的影響?,F采用矩形窗函數對聲信號進行分幀,設置幀長為200 ms,比較RF、KNN、DT、MLP 在不同特征選取下的分類準確率(圖10)。

    圖10 不同特征的煤矸識別準確率Fig.10 Recognition accuracy of coal and gangue with different characteristics

    RF、KNN、MLP 選取64 維時頻域特征作為輸入可獲得較高的準確率, 分別為 93.06% 、90.02% 、90.92% 。 DT 選取時域特征作為輸入可獲得較高的準確率,為89.08% ,但未達到90% ,不作考慮。 因此,選取時頻域特征,對聲信號樣本進行小波包分解,可提取更具備區(qū)分性的特征,有助于提升煤矸識別準確率。

    5.3 特征維度對于分類準確率的影響

    時頻域特征向量的維度會影響分類的準確率,本文利用主成分分析( Principal Component Analysis, PCA)[27]對時頻域特征向量進行降維,將高維的特征向量映射到低維空間中,增大樣本采樣密度,降低運算耗時,減弱信息重疊現象對分類模型性能的影響。

    本實驗設置窗函數為矩形窗,最佳幀長為200 ms,以時頻域特征作為輸入特征向量,共有64 維,對其進行主成分分析降維的結果共有64種。 輸入不同維度的特征向量,基于不同的機器學習算法進行模型訓練,采用交叉驗證法測試煤矸分類準確率。 圖11 為不同維度下多種分類模型的煤矸識別準確率。 由圖11 可知,RF、KNN、DT、MLP 模型分類準確率在降維的過程中均有提升,在64 維特征向量中存在較多無區(qū)分性特征,準確率分別為93.06% 、90.56% 、85.87% 、91.05% 。在降維至1 維時,由于損失過多的有效特征,準確率分別為54.91% 、58.36% 、54.89% 、59.00% 。 4種模型在降維過程中的煤矸分類準確率最大值分別為94.51% 、90.85% 、88.27% 、91.96% ,RF 模型分類準確率略高,可有效提升煤矸分類準確率。 當特征維度為20 維時,采用RF 分類模型,煤矸分類準確率最高,為94.51% 。

    圖11 不同維度下多種分類模型的煤矸識別準確率Fig.11 Accuracy of coal and gangue recognition by various classification models in different dimensions

    6 結 論

    本文研究了放頂煤聲信號的特征提取方法,對比了主流機器學習算法的分類表現,得出以下結論:

    (1) 不同幀長的樣本選取會對煤矸分類的準確率產生影響。 選擇幀長為200 ms、幀移為100 ms 時,基于時域特征提取與時頻域特征提取的分類準確率最高,分別為91.60%與93.35%。

    (2) 不同的特征選擇適用于不同的算法模型。本文的時頻域特征提取方法適用于RF、KNN 與MLP 分類器,煤矸分類準確率分別為93.06% 、90.05% 、91.02% 。 相比于時域特征和頻域特征,時頻域特征能提高分類準確率。

    (3) 對于時頻域特征向量,降維有助于提升分類準確率。 RF、KNN、DT、MLP 4 種模型的分類準確率分別提升了1.45% 、0.29% 、2.40% 、0.91% ,其中RF 分類模型準確率最高,為94.51% ,對應的時頻域特征向量維度為20。

    (4) RF、KNN、DT、MLP 在獲得較高分類準確率的同時,模型大小分別為6.60 MB、6.98 MB、0.05 MB 和0.01 MB。 綜合考慮分類準確率與模型大小,MLP 模型更具備低功耗平臺的適用性。

    猜你喜歡
    煤矸放頂頻域
    基于多通道數據耦合技術的智能機器人系統(tǒng)在煤矸分選中的研究與應用
    選煤技術(2023年4期)2024-01-02 03:11:16
    X 射線透射煤矸智能識別方法
    選煤廠封閉式自動煤矸分選機設計應用
    急傾斜特厚煤層水平分層綜采放頂煤工作面旋轉開采理論分析與實踐
    智富時代(2018年8期)2018-09-28 00:12:34
    頻域稀疏毫米波人體安檢成像處理和快速成像稀疏陣列設計
    雷達學報(2018年3期)2018-07-18 02:41:34
    基于改進Radon-Wigner變換的目標和拖曳式誘餌頻域分離
    一種基于頻域的QPSK窄帶干擾抑制算法
    放頂煤液壓支架尾梁裝置設計改進
    基于頻域伸縮的改進DFT算法
    電測與儀表(2015年3期)2015-04-09 11:37:24
    綜采放頂煤工作面快速拆除技術
    河南科技(2014年23期)2014-02-27 14:18:47
    成人18禁高潮啪啪吃奶动态图| 亚洲自拍偷在线| 两人在一起打扑克的视频| 波多野结衣高清作品| 变态另类成人亚洲欧美熟女| 成人18禁在线播放| 日韩av在线大香蕉| 男女做爰动态图高潮gif福利片| 免费无遮挡裸体视频| 日本在线视频免费播放| 一区二区三区精品91| 日本精品一区二区三区蜜桃| 日韩欧美国产一区二区入口| 国产精品乱码一区二三区的特点| 亚洲一区高清亚洲精品| 亚洲国产精品成人综合色| 色综合亚洲欧美另类图片| 亚洲人成网站在线播放欧美日韩| 中文在线观看免费www的网站 | 一级a爱视频在线免费观看| 亚洲 欧美一区二区三区| 高清毛片免费观看视频网站| 琪琪午夜伦伦电影理论片6080| 天天躁狠狠躁夜夜躁狠狠躁| 欧美日韩精品网址| 午夜福利成人在线免费观看| 国产精品久久视频播放| 岛国视频午夜一区免费看| 日本免费一区二区三区高清不卡| 搡老熟女国产l中国老女人| 午夜福利成人在线免费观看| 精品欧美国产一区二区三| 欧美成人性av电影在线观看| 黄片小视频在线播放| 男女午夜视频在线观看| 国产精品98久久久久久宅男小说| 亚洲av片天天在线观看| 国产精品电影一区二区三区| 中文字幕精品亚洲无线码一区 | 三级毛片av免费| 国产成人一区二区三区免费视频网站| 欧美人与性动交α欧美精品济南到| 中文字幕人妻熟女乱码| 99热6这里只有精品| 亚洲人成电影免费在线| 日韩有码中文字幕| 免费在线观看成人毛片| 久久精品国产99精品国产亚洲性色| 校园春色视频在线观看| 在线播放国产精品三级| 真人做人爱边吃奶动态| 搞女人的毛片| 亚洲一区中文字幕在线| 人人妻人人澡欧美一区二区| 久久中文字幕人妻熟女| 免费在线观看成人毛片| 一级毛片高清免费大全| 好男人电影高清在线观看| 亚洲成av片中文字幕在线观看| 中亚洲国语对白在线视频| 色av中文字幕| 侵犯人妻中文字幕一二三四区| 精华霜和精华液先用哪个| 国产精品爽爽va在线观看网站 | 国产精品98久久久久久宅男小说| 成人三级黄色视频| 人成视频在线观看免费观看| 免费看日本二区| 欧美日本视频| 女人高潮潮喷娇喘18禁视频| 91av网站免费观看| 人妻久久中文字幕网| 无限看片的www在线观看| 欧美日本亚洲视频在线播放| 精品免费久久久久久久清纯| 老司机福利观看| 成人一区二区视频在线观看| 午夜福利高清视频| a级毛片a级免费在线| 一级作爱视频免费观看| 亚洲av中文字字幕乱码综合 | 老司机午夜十八禁免费视频| 9191精品国产免费久久| 在线观看66精品国产| 婷婷六月久久综合丁香| 中国美女看黄片| 一卡2卡三卡四卡精品乱码亚洲| 日本熟妇午夜| 久久99热这里只有精品18| 日韩av在线大香蕉| 一卡2卡三卡四卡精品乱码亚洲| 欧美日韩瑟瑟在线播放| 久久99热这里只有精品18| 欧美成人性av电影在线观看| www.熟女人妻精品国产| 黑人巨大精品欧美一区二区mp4| 一区福利在线观看| 丝袜人妻中文字幕| 夜夜躁狠狠躁天天躁| 在线观看免费视频日本深夜| 在线永久观看黄色视频| www.www免费av| 91大片在线观看| 成人国产综合亚洲| 欧美午夜高清在线| 成年人黄色毛片网站| 两个人看的免费小视频| 老熟妇仑乱视频hdxx| 国产av不卡久久| 成年人黄色毛片网站| 最近最新中文字幕大全电影3 | 亚洲狠狠婷婷综合久久图片| 欧美色欧美亚洲另类二区| 免费一级毛片在线播放高清视频| 亚洲av成人不卡在线观看播放网| 国产精品av久久久久免费| 性欧美人与动物交配| 搞女人的毛片| 午夜a级毛片| 国产成人精品无人区| 精品欧美一区二区三区在线| 一二三四在线观看免费中文在| 亚洲片人在线观看| 国产精品 国内视频| 午夜久久久在线观看| 又黄又爽又免费观看的视频| 国产高清激情床上av| 久久天躁狠狠躁夜夜2o2o| 成年女人毛片免费观看观看9| 91在线观看av| 丁香欧美五月| 一级片免费观看大全| 91麻豆av在线| 又黄又爽又免费观看的视频| 在线免费观看的www视频| 俺也久久电影网| 久久久久精品国产欧美久久久| 亚洲va日本ⅴa欧美va伊人久久| 人人妻人人澡欧美一区二区| 制服丝袜大香蕉在线| 成年免费大片在线观看| 久久久久国产精品人妻aⅴ院| 国产真人三级小视频在线观看| 国产高清有码在线观看视频 | 淫妇啪啪啪对白视频| 久久久久久免费高清国产稀缺| 亚洲,欧美精品.| 亚洲色图 男人天堂 中文字幕| 怎么达到女性高潮| 国产三级在线视频| 一级片免费观看大全| 中文在线观看免费www的网站 | 国产成年人精品一区二区| 在线观看免费日韩欧美大片| 中文字幕高清在线视频| av在线播放免费不卡| a级毛片在线看网站| 久久久久精品国产欧美久久久| 日本黄色视频三级网站网址| 露出奶头的视频| 热99re8久久精品国产| 91在线观看av| 亚洲片人在线观看| 又大又爽又粗| 啦啦啦观看免费观看视频高清| 天天躁夜夜躁狠狠躁躁| 两个人视频免费观看高清| 可以免费在线观看a视频的电影网站| 久久久久免费精品人妻一区二区 | 久久久精品国产亚洲av高清涩受| 国产成+人综合+亚洲专区| 亚洲自拍偷在线| 黄频高清免费视频| 黄色丝袜av网址大全| 亚洲 欧美 日韩 在线 免费| 欧美日韩中文字幕国产精品一区二区三区| 少妇的丰满在线观看| 激情在线观看视频在线高清| 一二三四社区在线视频社区8| 十八禁人妻一区二区| 观看免费一级毛片| 国产高清视频在线播放一区| 淫秽高清视频在线观看| 1024手机看黄色片| 观看免费一级毛片| 韩国av一区二区三区四区| 亚洲精品国产区一区二| 波多野结衣高清作品| 美女 人体艺术 gogo| 看免费av毛片| 国产亚洲精品第一综合不卡| 亚洲av第一区精品v没综合| 精品日产1卡2卡| 欧美 亚洲 国产 日韩一| 国产三级在线视频| 一边摸一边抽搐一进一小说| 中文字幕av电影在线播放| 亚洲专区中文字幕在线| 日韩欧美免费精品| 色尼玛亚洲综合影院| 国产精品亚洲一级av第二区| 欧美激情 高清一区二区三区| 欧美人与性动交α欧美精品济南到| 亚洲av中文字字幕乱码综合 | 法律面前人人平等表现在哪些方面| 一级a爱视频在线免费观看| 99精品欧美一区二区三区四区| 制服诱惑二区| av中文乱码字幕在线| 久久精品国产亚洲av高清一级| 在线播放国产精品三级| 午夜福利在线在线| 麻豆av在线久日| 国产精品野战在线观看| 久久国产精品男人的天堂亚洲| 一边摸一边抽搐一进一小说| 国产伦在线观看视频一区| 国产主播在线观看一区二区| 欧美乱妇无乱码| 热re99久久国产66热| 午夜福利欧美成人| 久久九九热精品免费| 国产1区2区3区精品| 精品久久久久久,| 视频在线观看一区二区三区| 久久久久久久久中文| 国产一区在线观看成人免费| 一级作爱视频免费观看| 老司机深夜福利视频在线观看| 国产精品久久久人人做人人爽| 欧美人与性动交α欧美精品济南到| 91九色精品人成在线观看| 国产精品久久久av美女十八| 成年人黄色毛片网站| netflix在线观看网站| 波多野结衣av一区二区av| 18美女黄网站色大片免费观看| 人妻丰满熟妇av一区二区三区| 日本五十路高清| 久久精品成人免费网站| 亚洲一区二区三区不卡视频| 亚洲av美国av| 国产精品美女特级片免费视频播放器 | 深夜精品福利| 亚洲欧美激情综合另类| 亚洲精品在线美女| 精品国产乱码久久久久久男人| 亚洲精品国产一区二区精华液| 动漫黄色视频在线观看| 成在线人永久免费视频| 日本三级黄在线观看| 国产高清视频在线播放一区| 免费看十八禁软件| 一区二区三区国产精品乱码| 亚洲片人在线观看| 99精品久久久久人妻精品| 一本大道久久a久久精品| 国产欧美日韩一区二区三| 麻豆国产av国片精品| 中文资源天堂在线| 99在线人妻在线中文字幕| 亚洲成国产人片在线观看| 亚洲 欧美 日韩 在线 免费| 日韩av在线大香蕉| 韩国av一区二区三区四区| 91麻豆精品激情在线观看国产| 99热6这里只有精品| 变态另类成人亚洲欧美熟女| 欧美激情高清一区二区三区| 久久青草综合色| 午夜福利欧美成人| 黄色片一级片一级黄色片| 免费无遮挡裸体视频| 欧美性猛交黑人性爽| 国产久久久一区二区三区| 亚洲专区中文字幕在线| 大香蕉久久成人网| 欧美激情久久久久久爽电影| a级毛片a级免费在线| 欧美大码av| 国产亚洲精品第一综合不卡| 真人一进一出gif抽搐免费| 天天一区二区日本电影三级| 琪琪午夜伦伦电影理论片6080| 亚洲avbb在线观看| 成人av一区二区三区在线看| 亚洲精品中文字幕一二三四区| 免费在线观看日本一区| 一区福利在线观看| 欧美丝袜亚洲另类 | 久久精品国产综合久久久| 午夜免费观看网址| 国产单亲对白刺激| 欧美成人午夜精品| 精品免费久久久久久久清纯| 久久精品国产综合久久久| 美女高潮喷水抽搐中文字幕| 中出人妻视频一区二区| 少妇粗大呻吟视频| 1024视频免费在线观看| 女人爽到高潮嗷嗷叫在线视频| 国产精品免费一区二区三区在线| 精品欧美一区二区三区在线| 成人国产一区最新在线观看| 满18在线观看网站| 成人精品一区二区免费| 亚洲男人天堂网一区| 国产精品亚洲av一区麻豆| 精品国产超薄肉色丝袜足j| videosex国产| 可以免费在线观看a视频的电影网站| 国产乱人伦免费视频| 亚洲人成77777在线视频| 91大片在线观看| www日本黄色视频网| 一级作爱视频免费观看| 欧美日韩精品网址| 哪里可以看免费的av片| 亚洲色图 男人天堂 中文字幕| 欧美人与性动交α欧美精品济南到| 亚洲成a人片在线一区二区| 亚洲专区字幕在线| 亚洲欧美日韩无卡精品| 19禁男女啪啪无遮挡网站| 在线观看免费日韩欧美大片| 国产激情欧美一区二区| 久久久久亚洲av毛片大全| 成人午夜高清在线视频 | 欧美在线黄色| 久久久久国内视频| 两人在一起打扑克的视频| 高潮久久久久久久久久久不卡| 亚洲 国产 在线| 免费看美女性在线毛片视频| 精品午夜福利视频在线观看一区| 天天一区二区日本电影三级| 亚洲国产精品sss在线观看| www日本黄色视频网| 精品少妇一区二区三区视频日本电影| 午夜福利在线在线| 深夜精品福利| 久久久国产成人免费| 国产成人啪精品午夜网站| 老司机深夜福利视频在线观看| 亚洲专区中文字幕在线| 听说在线观看完整版免费高清| 久久精品国产99精品国产亚洲性色| 老司机在亚洲福利影院| 亚洲av美国av| 午夜激情福利司机影院| 亚洲国产日韩欧美精品在线观看 | 免费看美女性在线毛片视频| 亚洲一区高清亚洲精品| 成人国产综合亚洲| 国产精品二区激情视频| 亚洲成人久久爱视频| 国产熟女xx| 啦啦啦免费观看视频1| 这个男人来自地球电影免费观看| 一级毛片精品| 啦啦啦免费观看视频1| 亚洲国产欧美一区二区综合| 高清毛片免费观看视频网站| 欧美日韩中文字幕国产精品一区二区三区| 国内精品久久久久久久电影| 岛国在线观看网站| 欧美黄色片欧美黄色片| 可以在线观看的亚洲视频| 中文资源天堂在线| 一区二区日韩欧美中文字幕| 亚洲七黄色美女视频| 成人三级做爰电影| 欧美绝顶高潮抽搐喷水| 亚洲国产看品久久| 午夜精品在线福利| 搡老熟女国产l中国老女人| av在线播放免费不卡| 这个男人来自地球电影免费观看| 日日摸夜夜添夜夜添小说| 波多野结衣av一区二区av| 日本免费a在线| 级片在线观看| 18禁观看日本| 亚洲七黄色美女视频| 伦理电影免费视频| 免费在线观看影片大全网站| 黑丝袜美女国产一区| 国产在线精品亚洲第一网站| 国产在线观看jvid| 好男人在线观看高清免费视频 | 久久婷婷人人爽人人干人人爱| 97超级碰碰碰精品色视频在线观看| 悠悠久久av| 黄片大片在线免费观看| 又紧又爽又黄一区二区| 国产精品久久久久久人妻精品电影| 日韩国内少妇激情av| 制服丝袜大香蕉在线| 黑人操中国人逼视频| 一二三四社区在线视频社区8| 九色国产91popny在线| 国产欧美日韩一区二区三| 精品欧美一区二区三区在线| 美女午夜性视频免费| 老司机福利观看| 亚洲欧美日韩无卡精品| 中出人妻视频一区二区| 人人妻,人人澡人人爽秒播| 精品熟女少妇八av免费久了| 久久精品91蜜桃| 啪啪无遮挡十八禁网站| 欧美日韩瑟瑟在线播放| 黄频高清免费视频| 一本大道久久a久久精品| 欧美性猛交黑人性爽| 窝窝影院91人妻| 宅男免费午夜| 天堂影院成人在线观看| 欧美黑人精品巨大| 国产又黄又爽又无遮挡在线| 欧美黑人欧美精品刺激| 国产在线精品亚洲第一网站| 在线免费观看的www视频| 级片在线观看| 少妇被粗大的猛进出69影院| 色尼玛亚洲综合影院| 亚洲片人在线观看| 正在播放国产对白刺激| 国产一卡二卡三卡精品| 最近最新中文字幕大全电影3 | 久久这里只有精品19| 一区二区三区高清视频在线| 成人手机av| 亚洲中文字幕日韩| 欧美色视频一区免费| 亚洲中文字幕一区二区三区有码在线看 | 精品久久久久久久久久久久久 | av福利片在线| 丝袜在线中文字幕| 制服诱惑二区| 老熟妇仑乱视频hdxx| 亚洲人成电影免费在线| 日韩免费av在线播放| 欧美人与性动交α欧美精品济南到| 999精品在线视频| 真人一进一出gif抽搐免费| 久久婷婷成人综合色麻豆| 麻豆成人午夜福利视频| 久久婷婷成人综合色麻豆| 男人操女人黄网站| 自线自在国产av| 欧美日韩一级在线毛片| 久久久久亚洲av毛片大全| 欧美黄色淫秽网站| 美女大奶头视频| 两人在一起打扑克的视频| 亚洲自拍偷在线| www国产在线视频色| 久久国产精品男人的天堂亚洲| 亚洲国产中文字幕在线视频| 午夜福利18| 国产v大片淫在线免费观看| 美女大奶头视频| 国产v大片淫在线免费观看| 色在线成人网| 国产单亲对白刺激| 日本 av在线| 在线播放国产精品三级| 国产欧美日韩一区二区三| 久久午夜综合久久蜜桃| 天堂动漫精品| 人人妻人人看人人澡| 国产精品精品国产色婷婷| 熟妇人妻久久中文字幕3abv| 男女视频在线观看网站免费 | 91麻豆精品激情在线观看国产| 精品午夜福利视频在线观看一区| 欧美黑人精品巨大| 国产精品一区二区免费欧美| 女生性感内裤真人,穿戴方法视频| 高清在线国产一区| 亚洲男人天堂网一区| 日日夜夜操网爽| 欧美日韩中文字幕国产精品一区二区三区| 91在线观看av| 亚洲av第一区精品v没综合| 欧美性长视频在线观看| 成人午夜高清在线视频 | 精品欧美国产一区二区三| 99久久久亚洲精品蜜臀av| 亚洲色图 男人天堂 中文字幕| 每晚都被弄得嗷嗷叫到高潮| 在线看三级毛片| www日本黄色视频网| 一卡2卡三卡四卡精品乱码亚洲| 久久99热这里只有精品18| 国产在线精品亚洲第一网站| 日韩三级视频一区二区三区| av电影中文网址| 日韩av在线大香蕉| 999久久久精品免费观看国产| 99精品久久久久人妻精品| 成人国产一区最新在线观看| 国产精品亚洲av一区麻豆| 国产精品香港三级国产av潘金莲| 最近最新中文字幕大全电影3 | 亚洲欧洲精品一区二区精品久久久| 免费电影在线观看免费观看| xxx96com| a级毛片a级免费在线| 99国产精品一区二区三区| 麻豆av在线久日| 成人av一区二区三区在线看| 色综合婷婷激情| 嫩草影院精品99| 亚洲国产日韩欧美精品在线观看 | 妹子高潮喷水视频| 在线观看免费午夜福利视频| 精品人妻1区二区| 99久久综合精品五月天人人| 丰满的人妻完整版| 亚洲精品av麻豆狂野| 久久久国产精品麻豆| 搡老熟女国产l中国老女人| 亚洲av熟女| 国产视频一区二区在线看| 91国产中文字幕| 亚洲国产欧美日韩在线播放| 黄色片一级片一级黄色片| 一级a爱片免费观看的视频| 高清在线国产一区| 成人免费观看视频高清| 欧美绝顶高潮抽搐喷水| 视频在线观看一区二区三区| 国内少妇人妻偷人精品xxx网站 | 久久精品成人免费网站| 老汉色∧v一级毛片| 国产精品亚洲av一区麻豆| 国产主播在线观看一区二区| 级片在线观看| 亚洲av日韩精品久久久久久密| www日本在线高清视频| a在线观看视频网站| 91字幕亚洲| 色av中文字幕| 十八禁网站免费在线| 午夜日韩欧美国产| 老鸭窝网址在线观看| 悠悠久久av| av片东京热男人的天堂| 久久婷婷人人爽人人干人人爱| 久久 成人 亚洲| 日本免费a在线| 在线观看66精品国产| 久久香蕉激情| 日韩欧美国产一区二区入口| 狂野欧美激情性xxxx| 婷婷六月久久综合丁香| 国产激情欧美一区二区| 久久精品国产亚洲av高清一级| 国产视频一区二区在线看| 精品无人区乱码1区二区| 婷婷六月久久综合丁香| 波多野结衣巨乳人妻| 亚洲一码二码三码区别大吗| 国产一区二区三区视频了| 亚洲一区高清亚洲精品| 久久久久久久久久黄片| 亚洲国产日韩欧美精品在线观看 | 免费在线观看影片大全网站| 18禁黄网站禁片午夜丰满| 午夜福利在线在线| 欧美丝袜亚洲另类 | 两个人看的免费小视频| 亚洲av成人av| 波多野结衣巨乳人妻| 国产精品精品国产色婷婷| 欧美绝顶高潮抽搐喷水| 麻豆av在线久日| 午夜两性在线视频| 欧美日韩福利视频一区二区| 免费看a级黄色片| e午夜精品久久久久久久| 欧美日本亚洲视频在线播放| 欧美一级毛片孕妇| 国产一区在线观看成人免费| 正在播放国产对白刺激| 19禁男女啪啪无遮挡网站| 亚洲国产日韩欧美精品在线观看 | 18禁黄网站禁片午夜丰满| 91字幕亚洲| 狂野欧美激情性xxxx| 国产三级在线视频| 女警被强在线播放| 999久久久国产精品视频| 成人三级做爰电影| 一进一出好大好爽视频| 成人一区二区视频在线观看| 久久人妻福利社区极品人妻图片| 国产乱人伦免费视频| 在线观看午夜福利视频| 黄色 视频免费看| 午夜激情福利司机影院| 啦啦啦免费观看视频1| 日韩欧美国产在线观看| 国产av不卡久久| 不卡av一区二区三区| 中亚洲国语对白在线视频| 亚洲国产欧美日韩在线播放| 欧美中文日本在线观看视频| 草草在线视频免费看| 成人18禁在线播放| 久久久水蜜桃国产精品网| 欧美国产精品va在线观看不卡| 18禁裸乳无遮挡免费网站照片 |