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

    基于EMD和形態(tài)分形維數(shù)的微震波形識別

    2017-10-14 14:02:26趙國彥鄧青林李夕兵董隴軍陳光輝張楚旋
    中南大學學報(自然科學版) 2017年1期
    關鍵詞:微震維數(shù)形態(tài)學

    趙國彥,鄧青林,李夕兵,董隴軍,陳光輝,張楚旋

    ?

    基于EMD和形態(tài)分形維數(shù)的微震波形識別

    趙國彥,鄧青林,李夕兵,董隴軍,陳光輝,張楚旋

    (中南大學資源與安全工程學院,湖南長沙,410083)

    針對現(xiàn)有礦山微震監(jiān)測系統(tǒng)信號自動識別難的問題,提出基于經(jīng)驗模態(tài)分解(EMD)和形態(tài)分形維數(shù)的識別方法。首先,采用EMD將原始信號分解為若干個本征模態(tài)分量(IMF),選擇前5個分量進行重構(gòu)得到新的信號。其次,求出處理后信號的形態(tài)學分形維數(shù),利用微震波形和爆破波形分形維數(shù)的差異進行信號識別。對50組微震波形和50組爆破波形進行試驗研究。對比未經(jīng)EMD處理的形態(tài)學分形維數(shù)以及經(jīng)EMD處理的盒維數(shù)識別結(jié)果。研究結(jié)果表明:50組微震波形和50組爆破波形在形態(tài)學分形維數(shù)為1.4時具有較高的識別率;微震波形維數(shù)主要在1.4以下,爆破波形維數(shù)則基本高于1.4;EMD結(jié)合形態(tài)學分形維數(shù)的識別效果最好,為微震監(jiān)測波形識別提供了新途徑。

    微震;EMD;形態(tài)學;分形維數(shù);波形識別

    微震監(jiān)測技術通過對采集的地震波信息進行分析處理確定破裂點的空間位置和震級,根據(jù)累積微震事件的數(shù)量和能量釋放的關系,估計因微震活動引起的非線性應變區(qū)域的發(fā)展,據(jù)此判斷巖體的穩(wěn)定性。該技術對于保障礦山人員安全及設備完整至關重要,目前已在國內(nèi)外采礦工程中推廣應用[1?5]。然而,傳感器采集的原始數(shù)據(jù)往往摻雜著大量干擾的信息,需要后期人為剔除。為減輕人工勞動量,近年來許多學者對波形的自動識別進行研究,極大地推動了微震監(jiān)測技術的發(fā)展。歸結(jié)起來,微震波形的自動識別特征主要包括頻譜特征、能量分布特征和參數(shù)化特征,分析方法主要有快速傅里葉變換、小波及小波包變換、經(jīng)驗模態(tài)分解EMD等[6?11]。然而,一些分析方法本身存在缺陷,可能會帶來不準確的分析結(jié)果,如傅里葉變換存在時域和頻域的局部化矛盾,小波及小波包變換則受小波基選擇以及信號分解層次的限制。根據(jù)信號本身的局部特征時間尺度,將1個復雜振動信號分解為若干分量之和的經(jīng)驗模態(tài)分解方法具有較大優(yōu)勢,該方法不需要對信號具備先驗知識,是一種自適應的信號處理方式。此外,地震波在巖體內(nèi)傳播時,由于巖層的散射、反射等效應以及巖體對于地震波具有吸收作用,其能量與頻率是不斷衰減的[12],因此,無論采用頻率識別方式還是能量識別方式,均存在較大的局限性和隨機性。但巖體微震與爆破振動在觸發(fā)機制上存在本質(zhì)區(qū)別,直觀表達出來的形式就是它們的波形差異。分形維數(shù)則能定性、定量地刻畫這種非線性狀態(tài),可用于識別這種差異,但要獲得準確的分形維數(shù)則需要對信號進行濾波處理。為此,本文作者在借助經(jīng)驗模態(tài)分解構(gòu)造自適應濾波器功能的前提下,提取微震及爆破波形的主要信息,并且選擇精確度較高的形態(tài)學分形算法對微震和爆破波形的分形維數(shù)進行估計。

    1 EMD方法

    對于任意的復雜信號(),HUANG等[13]指出通過其本身特征時間尺度來獲得一組本征模態(tài)函數(shù)IMF,稱為經(jīng)驗模態(tài)分解篩選的方法,其步驟如下。

    1) 確定原始信號()的所有的局部極大值點、極小值點,并分別用3次樣條曲線擬合形成上、下包 絡線。

    2) 求出原始的數(shù)據(jù)序列與上下包絡線均值()之差:

    3) 檢驗式(1)得到的()是否滿足IMF的條件。若不滿足,則將其視為新的(),重復以上步驟直到滿足IMF的條件為止,此時()為第1階IMF,記為1(),通常包含信號的最高頻成分。

    4) 求出原始信號()與1()之差得到剩余序列1():

    5) 將1()作為新的原始信號,重復以上步驟,直到第階的剩余信號R()成為單調(diào)函數(shù)或者比預期值小,不能再篩分為止,即可得到一系列IMF分量1(),2(),…,n()以及剩余項n(),這樣初始信號()就分解成若干IMF分量與1個殘余項之和,即

    (3)

    式中:余項n()為信號的趨勢項,各IMF分量則表示信號從高到低不同頻率段的成分。

    EMD將任意復雜的信號自適應地分解為一組本征模態(tài)分量IMF,需要滿足以下2個條件:

    1) 在整個時間數(shù)據(jù)序列之中,極值點的數(shù)量(包括極大值與極小值)與過零點的數(shù)量至多相差1個。

    2) 在任意1個時間點上,由局部極值點構(gòu)成的包絡線平均值為0(局部極大值點確定上包絡線,局部極小值點確定下包絡線)。

    在EMD分解過程中,因?qū)嶋H的篩分很難保證局部均值為0,而過多的重復篩選會導致IMF失去實際的物理意義,可通過限制標準偏差來實現(xiàn)篩分停止:

    式中:為信號的時間長度;I?1()和I()為2個連續(xù)處理結(jié)果的時間序列;D通常取0.2~0.3。

    2 形態(tài)學分形維數(shù)及其計算

    2.1 分形維數(shù)

    分形是自然界復雜事物存在的普遍特征。非平穩(wěn)隨機振動信號的不規(guī)則程度能夠以分形維數(shù)來描述,計算分形維數(shù)的常見方法有盒維數(shù)、網(wǎng)格維數(shù)、關聯(lián)維數(shù)等,其中基于“覆蓋”思想的算法簡單清晰。該方法假設1個實平面集合2,和分別為集合中1個非空子集和1個尺度為的閉集,則的分形維數(shù)定義為[14]

    式中:()為覆蓋所需的個數(shù)。

    2.2 形態(tài)學分形算法

    估計分形維數(shù)的關鍵是在不同尺度下對信號進行度量,利用多尺度形態(tài)學可以實現(xiàn)[15]。在集合覆蓋信號過程中,以其上確界函數(shù)即結(jié)構(gòu)元素()進行等效變換,分析尺度范圍=1, 2, …,,其算法包含膨脹及腐蝕2種基本算子。對于一維離散時間信號()(=0, 1, 2, …),每個尺度的膨脹和腐蝕結(jié)果分別為:

    (7)

    (9)

    式中:M為信號的Minkowski-Bouligand維數(shù);為常數(shù)。對上式和進行最小二乘擬合,得到的直線斜率為M,即為最終所要求得的()分形維數(shù)。

    2.3 參數(shù)選擇

    一般地,選擇單位結(jié)構(gòu)元素()為{0,0,0},因該結(jié)構(gòu)在保證維數(shù)估計不受信號幅值范圍影響的同時,還能使算法的運算量減少。最大分析尺度原則上為小于/2的正整數(shù)(為離散信號采樣點數(shù)),當數(shù)據(jù)長度比較大時,適當減小可減少計算量[16],本文取256。

    3 基于EMD和形態(tài)分形維數(shù)的微震波形識別

    3.1 原理及步驟

    從本質(zhì)上來講,礦山微震以及爆破產(chǎn)生的波形均屬于非平穩(wěn)隨機振動(如圖1和圖2所示),因而其本身已經(jīng)具備可估計分形維數(shù)的特點。然而,由于礦井下復雜作業(yè)環(huán)境的影響,系統(tǒng)采集到的波形往往受到各種噪聲的污染,因此,直接對原始波形進行分形維數(shù)估計會得到不準確的結(jié)果,必須事先對波形進行去噪預處理。

    圖1 典型的巖體微震波形

    圖2 典型的爆破振動波形

    經(jīng)驗模態(tài)分解EMD是一種無需先驗知識的非平穩(wěn)振動信號處理方法,可將信號自適應地分解為有限個IMF,IMF按信息重要程度由高到低階排列,且由于EMD算法本身還存在一定缺陷,容易產(chǎn)生虛假分量[17?18],因此,本文在結(jié)合大量試驗結(jié)果的前提下選擇前5階IMF提取信號的主要信息,作為去噪后的最終信號。綜上所述,基于EMD和形態(tài)分形維數(shù)的微震波形識別步驟為:

    1) 利用EMD將原始信號()分解為若干個IMF;

    由此計算出的典型微震與爆破波形形態(tài)學分形維數(shù)分別為1.292 9和1.543 2,可見兩者差異較大。

    3.2 實例應用

    為尋找微震波形與爆破波形的分形特征差異,確定兩者形態(tài)學維數(shù)的定量差異,分別從原始波形數(shù)據(jù)庫隨機選取50組微震事件波形、50組爆破事件波形,利用matlab編程,計算各個波形的形態(tài)學分形維數(shù),見表1。

    表1 微震波形與爆破波形形態(tài)學分形維數(shù)

    從表1可以看出:微震事件波形的形態(tài)學分形維數(shù)基本上低于爆破事件波形的形態(tài)學分形維數(shù)。因此,可以用形態(tài)學分形維數(shù)來辨識這2類波形。圖3所示為微震及爆破波形的形態(tài)學分形維數(shù)分布圖。

    圖3 微震與爆破波形的形態(tài)學分形維數(shù)分布

    從圖3可知:微震波形與爆破波形的形態(tài)學分形維數(shù)為1.4時呈現(xiàn)明顯分界,分形維數(shù)在1.4以上以爆破事件為主,1.4以下以微震事件為主。若以形態(tài)學分形維數(shù)1.4作為辨識微震與爆破事件的判別因子,那么,在抽取的樣本中,分形維數(shù)高于1.4的微震事件有4個,識別正確率為92%,而分形維數(shù)低于1.4的爆破事件則有3個,識別正確率為94%;總體事件識別的正確率達到93%。從波形復雜度的角度看,微震事件產(chǎn)生的波形大多數(shù)較單一,而井下生產(chǎn)爆破多有多段微差延時,波形會相對復雜。因此,以形態(tài)學分形維數(shù)1.4作為該礦山微震系統(tǒng)波形識別因子是可行的。

    3.3 結(jié)果驗證

    為驗證形態(tài)學分形維數(shù)辨識的效果,一方面,以未經(jīng)EMD處理的原始波形形態(tài)學分形維數(shù)比較,另一方面,選擇分形常用的盒維數(shù)計算方法進行對比。圖4和圖5所示分別為未處理的原始波形形態(tài)學分形維數(shù)分布圖以及用盒維數(shù)代替形態(tài)學分形維數(shù)的微震與爆破波形維數(shù)分布圖。

    從圖4和圖5可見:微震與爆破原始波形的形態(tài)學維數(shù)以及經(jīng)EMD的分形盒維數(shù)在總體上也呈爆破高、微震低的趨勢,但無論是以形態(tài)學分形維數(shù)1.4~1.5還是以盒維數(shù)1.30~1.32中任意1個維數(shù)作為兩者的識別因子,其總體正確識別率均低于90%,說明不經(jīng)過預處理的波形以及采用本身估計精度較低的盒維數(shù)計算方法,均可能影響最終結(jié)果,使得高維數(shù)微震波形和低維數(shù)爆破波形容易產(chǎn)生混淆。以上分析結(jié)果表明了基于EMD的形態(tài)學分形維數(shù)識別微震和爆破波形的有效性及合理性。

    圖4 未經(jīng)EMD處理的微震與爆破波形的形態(tài)學分形維數(shù)分布

    圖5 微震與爆破波形的分形盒維數(shù)分布

    4 結(jié)論

    1) EMD具有自適應濾波器的功能,能夠提取振動信號的主要特征與信息,從而提高了原始信號的信噪比。

    2) 分形維數(shù)具有刻畫不同波形特征的能力?;谛螒B(tài)學分形維數(shù)的算法估計精度高,對經(jīng)EMD處理的100組樣本進行試驗,求出微震波形與爆破波形的識別維數(shù)為1.4,且識別正確率高達93%,相比未經(jīng)EMD處理的原始波形形態(tài)學分形維數(shù)以及經(jīng)EMD處理的盒維數(shù)識別,更加合理有效,可為后續(xù)大量樣本的模式自動識別提供一種新的方法。

    [1] SALVONI M, DIGHT P M. Rock damage assessment in a large unstable slope from microseismic monitoring-MMG century mine (Queensland, Australia) case study[J]. Engineering Geology, 2016, 210: 45?56.

    [2] URBANCIC T I, TRIFU C I. Recent advances in seismic monitoring technology at Canadian mines[J]. Journal of Applied Geophysics, 2000, 45(4): 225?237.

    [3] POTVIN Y. Strategies and tactics to control seismic risks in mines[J]. The Journal of the Southern African Institute of Mining and Metallurgy, 2009, 109(3): 177?186.

    [4] 楊承祥, 羅周全, 唐禮忠. 基于微震監(jiān)測技術的深井開采地壓活動規(guī)律研究[J].巖石力學與工程學報, 2007, 26(4): 818?824.YANG Chengxiang, LUO Zhouquan, TANG Lizhong. Study on rule of geostatic activity based on microseismic monitoring technique in deep mining[J]. Chinese Journal of Rock Mechanics and Engineering, 2007, 26(4): 818?824.

    [5] 董隴軍, 李夕兵, 唐禮忠, 等. 無需預先測速的微震震源定位的數(shù)學形式及震源參數(shù)確定[J].巖石力學與工程學報, 2011, 30(10): 2057?2067. DONG Longjun, LI Xibing, TANG Lizhong, et al. Mathematical functions and parameters for microseismic source location without pre-measuring speed[J]. Chinese Journal of Rock Mechanics and Engineering, 2011, 30(10): 2057?2067.

    [6] ZHAO Guoyan, MA Ju, DONG Longjun, et al. Classification of mine blasts and microseismic events using starting-up features in seismograms[J]. Transactions of Nonferrous Metals Society of China, 2015, 25(10): 3410?3420.

    [7] 趙國彥, 鄧青林, 馬舉. 基于FSWT時頻分析的礦山微震信號分析與識別[J]. 巖土工程學報, 2015, 37(2): 306?312. ZHAO Guoyan, DENG Qinglin, MA Ju. Recognition of mine microseismic signals based on FSWT time-frequency analysis[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(2): 306?312.

    [8] 陸菜平, 竇林名, 吳興榮, 等. 巖體微震監(jiān)測的頻譜分析與信號識別[J]. 巖土工程學報, 2005, 27(7): 772?775. LU Caiping, DOU Linming, WU Xingrong, et al. Frequency spectrum analysis on microseismic monitoring and signal differentiation of rock material[J]. Chinese Journal of Geotechnical Engineering, 2005, 27(7): 772?775.

    [9] 朱權潔, 姜福興, 于正興, 等. 爆破震動與巖石破裂微震信號能量分布特征研究[J]. 巖石力學與工程學報, 2012, 31(4): 723?730. ZHU Quanjie, JIANG Fuxing, YU Zhengxing, et al. study on energy distribution characters about blasting vibration and rock fracture microseismic signal[J]. Chinese Journal of Rock Mechanics and Engineering, 2012, 31(4): 723?730.

    [10] 朱權潔, 姜福興, 尹永明, 等. 基于小波分形特征與模式識別的礦山微震波形識別研究[J]. 巖土工程學報, 2012, 34(11): 2036?2042. ZHU Quanjie, JIANG Fuxing, YIN Yongming, et al. Classification of mine microseismic events based on wavelet-fractal method and pattern recognition[J]. Chinese Journal of Geotechnical Engineering, 2012, 34(11): 2036?2042.

    [11] 吳響, 錢建生, 王海燕, 等. 微震信號多尺度非線性特征提取與辨識研究[J]. 儀器儀表學報, 2014, 35(5): 969?975. WU Xiang, QIAN Jiansheng, WANG Haiyan, et al. Study on multi-scale nonlinear feature extraction and signal identification for microseismic signal[J]. Chinese Journal of Scientific Instrument, 2014, 35(5): 969?975.

    [12] 李振春, 王清振. 地震波衰減機理及能量補償研究綜述[J]. 地球物理學進展, 2007, 22(4): 1147?1152. LI Zhenchun, WANG Qingzhen. A review of research on mechanism of seismic attenuation and energy compensation[J]. Progress in Geophysics, 2007, 22(4): 1147?1152.

    [13] HUANG N E, SHEN Z, LONG S R, et al. The empirical mode decomposition and the hilbert spectrum for nonlinear and non-stationary time series analysis[C]// Proceedings of the Royal Society of A: Mathematical, Physical and Engineering Sciences. London: The Royal Society, 1998: 903?995.

    [14] 李兵, 張培林, 任國全, 等. 基于數(shù)學形態(tài)學的分形維數(shù)計算及在軸承故障診斷中的應用[J].振動與沖擊, 2010, 29(5): 191?194. LI Bing, ZHANG Peilin, REN Guoquan, et al. Mathematic morphology based fractal dimension calculation and its application in fault diagnosis of roller bearings[J]. Journal of Vibration and Shock, 2010, 29(5): 191?194.

    [15] 張亢, 程軍圣, 楊宇. 基于局部均值分解與形態(tài)學分形維數(shù)的滾動軸承故障診斷方法[J].振動與沖擊, 2013, 32(9): 90?94. ZHANG Kang, CHENG Junsheng, YANG Yu. Roller bearing fault diagnosis based on local mean decomposition and morphological fractal dimension[J]. Journal of Vibration and Shock, 2013, 32(9): 90?94.

    [16] 王鳳利, 段樹林, 于洪亮, 等. 基于EEMD和形態(tài)學分形維數(shù)的柴油機故障診斷[J].內(nèi)燃機學報, 2012, 30(6): 557?562. WANG Fengli, DUAN Shulin, YU Hongliang, et al. Fault diagnosis of diesel engine based on EEMD and morphology fractal dimension[J]. Transactions of CSICE, 2012, 30(6): 557?562.

    [17] 韓東穎, 李庚, 時培明. 基于EMD和分形盒維數(shù)的旋轉(zhuǎn)機械耦合故障診斷方法研究[J].振動與沖擊, 2013, 32(15): 209?214. HAN Dongying, LI Geng, SHI Peiming. Diagnosis method for coupling faults of rotary machinery based on EMD and fractal box dimension[J]. Journal of Vibration and Shock, 2013, 32(15): 209?214.

    [18] 程軍圣, 于德介, 楊宇. 基于EMD和分形維數(shù)的轉(zhuǎn)子系統(tǒng)故障診斷[J].中國機械工程, 2005, 16(12): 1088?1091. CHENG Junsheng, YU Dejie, YANG Yu. Fault diagnosis for rotor system based on EMD and fractal dimension[J]. China Mechanical Engineering, 2005, 16(12): 1088?1091.

    (編輯 陳燦華)

    Recognition of microseismic waveforms based on EMD and morphological fractal dimension

    ZHAO Guoyan, DENG Qinglin, LI Xibing, DONG Longjun, CHEN Guanghui, ZHANG Chuxuan

    (School of Resources and Safety Engineering, Central South University, Changsha 410083, China)

    Considering that it is difficult to recognize signals automatically in microseismic monitoring system of mine, the new recognition method was proposed based on empirical mode decomposition (EMD) and morphological fractal dimension. Firstly, origined-signal was decomposed to some intrinsic mode functions (IMF) by EMD, and new signal was obtained by adding former five components. Then, morphological fractal dimension of the new signal was calculated to recognize microseismic waveform and blasting waveform. The 50 microseimic waveforms and 50 blasting waveforms were tested. Finally, the dimensions were calculated with two different methods, i.e. morphological fractal dimension without EMD and fractal box dimension with EMD. The results show that 50 microseimic waveforms and 50 blasting waveforms show high recognition ratio in morphological fractal dimension of 1.4. Dimensions of microseismic waveform and blasting waveform are mainly below and over 1.4, respectively. The new recognition method has the best recognition result, and provides a new way for waveform recognition.

    microseimic; EMD; morphology; fractal dimension; waveform recognition

    10.11817/j.issn.1672-7207.2017.01.023

    TD76;TD235

    A

    1672?7207(2017)01?0162?06

    2016?01?10;

    2016?03?24

    國家自然科學基金資助項目(51374244) (Project(51374244) supported by the National Natural Science Foundation of China)

    趙國彥,教授,博士生導師,從事采礦工程、礦山安全和巖石力學與工程等研究;E-mail: 810703752@qq.com

    猜你喜歡
    微震維數(shù)形態(tài)學
    β-變換中一致丟番圖逼近問題的維數(shù)理論
    淺談KJ768煤礦微震監(jiān)測系統(tǒng)的應用
    新疆鋼鐵(2021年1期)2021-10-14 08:45:48
    一類齊次Moran集的上盒維數(shù)
    長平煤業(yè)5302 綜放工作面頂板巖層移動規(guī)律研究
    魅力中國(2020年46期)2020-02-02 11:50:08
    關于齊次Moran集的packing維數(shù)結(jié)果
    基于隸屬度分析的回采面沖擊地壓微震能量閾值判定
    中國煤炭(2016年1期)2016-05-17 06:11:34
    涉及相變問題Julia集的Hausdorff維數(shù)
    醫(yī)學微觀形態(tài)學在教學改革中的應用分析
    數(shù)學形態(tài)學濾波器在轉(zhuǎn)子失衡識別中的應用
    巨核細胞數(shù)量及形態(tài)學改變在四種類型MPN中的診斷價值
    91老司机精品| 中文字幕高清在线视频| 久久精品国产亚洲av高清一级| 一区二区三区高清视频在线| 欧美另类亚洲清纯唯美| 香蕉丝袜av| 国产在线观看jvid| 国产亚洲精品av在线| 90打野战视频偷拍视频| 夜夜夜夜夜久久久久| 757午夜福利合集在线观看| 午夜免费观看网址| 天天躁夜夜躁狠狠躁躁| 亚洲无线在线观看| 国产精品一区二区免费欧美| 久久久国产欧美日韩av| 一级片免费观看大全| 国产精品野战在线观看| 国产蜜桃级精品一区二区三区| xxxwww97欧美| 亚洲精品在线美女| 午夜亚洲福利在线播放| 亚洲专区国产一区二区| 国内精品久久久久精免费| 非洲黑人性xxxx精品又粗又长| 日韩欧美三级三区| 亚洲中文字幕一区二区三区有码在线看 | 亚洲一区二区三区不卡视频| 国产精品亚洲av一区麻豆| 久久久国产精品麻豆| 女人高潮潮喷娇喘18禁视频| 精品高清国产在线一区| 在线观看免费午夜福利视频| 国产高清视频在线观看网站| 国产精品一区二区精品视频观看| 两性夫妻黄色片| 色精品久久人妻99蜜桃| 91麻豆av在线| 一边摸一边抽搐一进一小说| 欧美国产日韩亚洲一区| 国产亚洲欧美在线一区二区| 国产亚洲av高清不卡| 久久人妻av系列| 亚洲一卡2卡3卡4卡5卡精品中文| 变态另类丝袜制服| 亚洲国产高清在线一区二区三| 观看免费一级毛片| 成人手机av| 国产精品日韩av在线免费观看| 国产高清视频在线播放一区| 国内精品久久久久精免费| 高清毛片免费观看视频网站| 亚洲美女视频黄频| 国产99白浆流出| 脱女人内裤的视频| 国产成人aa在线观看| 搡老妇女老女人老熟妇| 在线观看一区二区三区| 国产亚洲精品久久久久久毛片| 一级毛片女人18水好多| 成年免费大片在线观看| 亚洲七黄色美女视频| 亚洲av美国av| 色在线成人网| 中文在线观看免费www的网站 | 免费在线观看完整版高清| 午夜福利成人在线免费观看| 亚洲成人中文字幕在线播放| 成人亚洲精品av一区二区| 999久久久精品免费观看国产| 母亲3免费完整高清在线观看| 婷婷亚洲欧美| 成人午夜高清在线视频| 日韩精品中文字幕看吧| 亚洲五月天丁香| 动漫黄色视频在线观看| 蜜桃久久精品国产亚洲av| 亚洲五月婷婷丁香| 床上黄色一级片| 亚洲欧美日韩无卡精品| 精品一区二区三区视频在线观看免费| 久久久久免费精品人妻一区二区| 在线播放国产精品三级| 国产麻豆成人av免费视频| 精品第一国产精品| 十八禁人妻一区二区| 在线免费观看的www视频| 精品少妇一区二区三区视频日本电影| 久久热在线av| 在线观看免费日韩欧美大片| 精品国产美女av久久久久小说| 国产99白浆流出| 男男h啪啪无遮挡| 午夜成年电影在线免费观看| 一级片免费观看大全| 亚洲av日韩精品久久久久久密| av在线播放免费不卡| 男人舔奶头视频| 女人高潮潮喷娇喘18禁视频| 波多野结衣高清无吗| 精品国产乱子伦一区二区三区| 1024手机看黄色片| 亚洲人成77777在线视频| 美女午夜性视频免费| 动漫黄色视频在线观看| 久久精品91蜜桃| 日韩精品中文字幕看吧| 亚洲中文日韩欧美视频| 两个人看的免费小视频| 高清在线国产一区| 亚洲精品国产一区二区精华液| 别揉我奶头~嗯~啊~动态视频| 精品久久久久久久末码| 一个人免费在线观看的高清视频| 亚洲,欧美精品.| 亚洲国产精品sss在线观看| 天天躁夜夜躁狠狠躁躁| 日本黄色视频三级网站网址| 精品福利观看| 国产精品一区二区三区四区免费观看 | 久久久久久久午夜电影| 国产又色又爽无遮挡免费看| 伊人久久大香线蕉亚洲五| 亚洲精品在线美女| 日本成人三级电影网站| 亚洲欧美日韩东京热| 国产亚洲精品综合一区在线观看 | 神马国产精品三级电影在线观看 | 少妇裸体淫交视频免费看高清 | 免费看十八禁软件| 亚洲av第一区精品v没综合| 免费电影在线观看免费观看| 国内揄拍国产精品人妻在线| 久久久久久人人人人人| 欧美日韩精品网址| 好看av亚洲va欧美ⅴa在| 黑人欧美特级aaaaaa片| 欧洲精品卡2卡3卡4卡5卡区| 国产69精品久久久久777片 | 日韩三级视频一区二区三区| 国产蜜桃级精品一区二区三区| 亚洲av成人不卡在线观看播放网| 国产精品98久久久久久宅男小说| 亚洲欧美日韩东京热| 国产亚洲精品av在线| 一级毛片高清免费大全| 久久久久精品国产欧美久久久| 欧美日韩福利视频一区二区| 国产午夜精品久久久久久| 一边摸一边做爽爽视频免费| 成人欧美大片| 中文字幕最新亚洲高清| www国产在线视频色| 国产97色在线日韩免费| 给我免费播放毛片高清在线观看| 男女之事视频高清在线观看| 波多野结衣高清作品| 国产免费av片在线观看野外av| 波多野结衣高清无吗| 日本黄大片高清| 日韩中文字幕欧美一区二区| 国产精品免费视频内射| 成年版毛片免费区| 又黄又爽又免费观看的视频| 老司机在亚洲福利影院| 亚洲一区中文字幕在线| www日本在线高清视频| 日日夜夜操网爽| 亚洲午夜精品一区,二区,三区| 国产97色在线日韩免费| 亚洲一卡2卡3卡4卡5卡精品中文| 青草久久国产| 亚洲精华国产精华精| 一个人观看的视频www高清免费观看 | 国产伦人伦偷精品视频| 亚洲18禁久久av| 男女床上黄色一级片免费看| 国产av一区二区精品久久| 美女大奶头视频| 亚洲成人中文字幕在线播放| 亚洲自拍偷在线| 在线视频色国产色| 日本五十路高清| 免费在线观看完整版高清| 亚洲aⅴ乱码一区二区在线播放 | 怎么达到女性高潮| www.精华液| 99热这里只有精品一区 | 香蕉av资源在线| 婷婷六月久久综合丁香| 国内揄拍国产精品人妻在线| 久久天堂一区二区三区四区| 国产高清视频在线观看网站| a级毛片a级免费在线| 亚洲av日韩精品久久久久久密| 波多野结衣高清作品| 免费高清视频大片| 18禁国产床啪视频网站| 在线永久观看黄色视频| 亚洲激情在线av| 国产免费av片在线观看野外av| 国产91精品成人一区二区三区| 欧美黄色片欧美黄色片| 国产精品av久久久久免费| 亚洲狠狠婷婷综合久久图片| 听说在线观看完整版免费高清| 免费看美女性在线毛片视频| 国产黄色小视频在线观看| 制服丝袜大香蕉在线| 亚洲黑人精品在线| 18禁美女被吸乳视频| 亚洲九九香蕉| 亚洲人成77777在线视频| 极品教师在线免费播放| 欧美性猛交╳xxx乱大交人| 国产伦人伦偷精品视频| 国内毛片毛片毛片毛片毛片| 欧美黑人巨大hd| 毛片女人毛片| 精品无人区乱码1区二区| 成人国语在线视频| 蜜桃久久精品国产亚洲av| 黑人欧美特级aaaaaa片| 精品一区二区三区av网在线观看| 麻豆一二三区av精品| 日韩欧美在线乱码| 精品国内亚洲2022精品成人| 亚洲中文日韩欧美视频| 精品国内亚洲2022精品成人| 男女那种视频在线观看| 国产aⅴ精品一区二区三区波| 两个人视频免费观看高清| 国产精品一区二区三区四区久久| 亚洲成人久久爱视频| 亚洲午夜精品一区,二区,三区| 亚洲第一电影网av| 在线观看午夜福利视频| 99久久国产精品久久久| 欧美丝袜亚洲另类 | 又黄又爽又免费观看的视频| 老熟妇乱子伦视频在线观看| netflix在线观看网站| 成人国产综合亚洲| 又紧又爽又黄一区二区| 亚洲熟女毛片儿| 日韩免费av在线播放| 欧美av亚洲av综合av国产av| 少妇人妻一区二区三区视频| 久久精品aⅴ一区二区三区四区| av国产免费在线观看| 国内揄拍国产精品人妻在线| 国产精品av视频在线免费观看| 久久中文看片网| 久久精品综合一区二区三区| 欧美日韩乱码在线| 国模一区二区三区四区视频 | 女生性感内裤真人,穿戴方法视频| 国产三级黄色录像| 国产精品av久久久久免费| 国产欧美日韩一区二区精品| av在线天堂中文字幕| 久久精品国产亚洲av高清一级| 国产高清有码在线观看视频 | 这个男人来自地球电影免费观看| 又大又爽又粗| 成人手机av| x7x7x7水蜜桃| 久久伊人香网站| 91大片在线观看| 一区二区三区国产精品乱码| 欧美精品亚洲一区二区| 成人av一区二区三区在线看| www.自偷自拍.com| 欧美黑人精品巨大| 亚洲精品美女久久av网站| 床上黄色一级片| 国产黄色小视频在线观看| 国产熟女午夜一区二区三区| 国产97色在线日韩免费| 国产精品永久免费网站| 欧美大码av| 99国产精品99久久久久| 岛国在线观看网站| 亚洲成人久久性| 午夜免费成人在线视频| 国产精品久久久久久精品电影| 国产亚洲欧美在线一区二区| 国产高清有码在线观看视频 | 校园春色视频在线观看| 午夜精品一区二区三区免费看| 又黄又粗又硬又大视频| 日本黄大片高清| av欧美777| 色综合站精品国产| 一边摸一边抽搐一进一小说| x7x7x7水蜜桃| 国产成人av激情在线播放| 在线观看舔阴道视频| 成人特级黄色片久久久久久久| 亚洲国产精品999在线| 亚洲av成人av| 国产精品 国内视频| 亚洲免费av在线视频| 人妻丰满熟妇av一区二区三区| 老司机在亚洲福利影院| 97碰自拍视频| 在线免费观看的www视频| 岛国在线免费视频观看| 久久香蕉激情| 久久久久精品国产欧美久久久| 午夜日韩欧美国产| 国产午夜福利久久久久久| 脱女人内裤的视频| 国产1区2区3区精品| 国产精品一区二区三区四区免费观看 | 桃红色精品国产亚洲av| 很黄的视频免费| 国内久久婷婷六月综合欲色啪| 香蕉久久夜色| 999久久久精品免费观看国产| 村上凉子中文字幕在线| 久久精品影院6| 亚洲精品中文字幕一二三四区| 18禁黄网站禁片午夜丰满| 91大片在线观看| 亚洲色图av天堂| 午夜福利在线在线| 丰满人妻一区二区三区视频av | a级毛片在线看网站| 成年免费大片在线观看| 欧美另类亚洲清纯唯美| 99精品在免费线老司机午夜| 久热爱精品视频在线9| 欧美又色又爽又黄视频| 十八禁网站免费在线| 人成视频在线观看免费观看| 国内精品久久久久久久电影| 两个人看的免费小视频| 亚洲色图av天堂| av天堂在线播放| 哪里可以看免费的av片| 日本 欧美在线| 国产单亲对白刺激| 看黄色毛片网站| 欧美又色又爽又黄视频| 琪琪午夜伦伦电影理论片6080| 日本一二三区视频观看| 在线观看免费午夜福利视频| 亚洲一区高清亚洲精品| 国产1区2区3区精品| 欧美三级亚洲精品| 亚洲精品在线美女| 欧美av亚洲av综合av国产av| 亚洲国产精品sss在线观看| www.999成人在线观看| 亚洲av第一区精品v没综合| 最新在线观看一区二区三区| 精品久久久久久久人妻蜜臀av| 精品少妇一区二区三区视频日本电影| 一级a爱片免费观看的视频| 久久人妻福利社区极品人妻图片| 午夜老司机福利片| 日日爽夜夜爽网站| 亚洲人成网站在线播放欧美日韩| 中文字幕人成人乱码亚洲影| 久久中文字幕人妻熟女| 国产成人欧美在线观看| 国产在线精品亚洲第一网站| 动漫黄色视频在线观看| 欧美黄色片欧美黄色片| 男插女下体视频免费在线播放| 欧美成人免费av一区二区三区| 9191精品国产免费久久| 神马国产精品三级电影在线观看 | 久久久国产成人精品二区| 久久久久久久久中文| 久久精品人妻少妇| 19禁男女啪啪无遮挡网站| 国产精品美女特级片免费视频播放器 | 十八禁网站免费在线| 丝袜人妻中文字幕| 国产激情久久老熟女| netflix在线观看网站| 国产三级在线视频| 国产黄色小视频在线观看| 一本精品99久久精品77| 老司机靠b影院| 国产黄a三级三级三级人| 麻豆国产av国片精品| 国产成人啪精品午夜网站| 欧美成人午夜精品| 国产一区二区激情短视频| 久久久久久亚洲精品国产蜜桃av| 国产精品免费视频内射| 黄色视频不卡| 日韩欧美在线二视频| 国产三级中文精品| 亚洲精华国产精华精| 日韩欧美在线乱码| 久久中文看片网| 久久这里只有精品19| 激情在线观看视频在线高清| 亚洲精品久久成人aⅴ小说| 成人午夜高清在线视频| 蜜桃久久精品国产亚洲av| 欧美激情久久久久久爽电影| 村上凉子中文字幕在线| 国产精品自产拍在线观看55亚洲| aaaaa片日本免费| 亚洲国产欧美网| 国产三级中文精品| 精品久久久久久久毛片微露脸| 18禁黄网站禁片午夜丰满| 91九色精品人成在线观看| 在线视频色国产色| 色播亚洲综合网| 日本在线视频免费播放| 中文亚洲av片在线观看爽| 亚洲成人国产一区在线观看| 欧美乱妇无乱码| 啦啦啦观看免费观看视频高清| 女同久久另类99精品国产91| 丰满人妻一区二区三区视频av | 91麻豆av在线| 亚洲自偷自拍图片 自拍| 男人舔女人下体高潮全视频| 国模一区二区三区四区视频 | 欧美中文综合在线视频| 小说图片视频综合网站| 久久中文看片网| 女警被强在线播放| 正在播放国产对白刺激| 日本黄大片高清| 一个人免费在线观看的高清视频| 国产视频内射| 曰老女人黄片| 久久精品aⅴ一区二区三区四区| 高清毛片免费观看视频网站| 亚洲欧美日韩高清在线视频| 好男人电影高清在线观看| 香蕉丝袜av| 美女扒开内裤让男人捅视频| 久久香蕉精品热| 精品乱码久久久久久99久播| 亚洲一卡2卡3卡4卡5卡精品中文| 黑人操中国人逼视频| 美女午夜性视频免费| 久久久久九九精品影院| 啦啦啦韩国在线观看视频| 欧美av亚洲av综合av国产av| 这个男人来自地球电影免费观看| 日日爽夜夜爽网站| 美女大奶头视频| 18禁裸乳无遮挡免费网站照片| 亚洲精品av麻豆狂野| 在线视频色国产色| 亚洲精品久久国产高清桃花| 可以在线观看毛片的网站| 久久国产乱子伦精品免费另类| 夜夜爽天天搞| 精品国内亚洲2022精品成人| 国产野战对白在线观看| 亚洲五月天丁香| 窝窝影院91人妻| 亚洲国产看品久久| 在线观看免费日韩欧美大片| 亚洲中文日韩欧美视频| 99re在线观看精品视频| 男插女下体视频免费在线播放| 看黄色毛片网站| 一a级毛片在线观看| 毛片女人毛片| 两性午夜刺激爽爽歪歪视频在线观看 | 亚洲熟妇熟女久久| 国产v大片淫在线免费观看| 亚洲熟女毛片儿| 国产精品av久久久久免费| 国产精品永久免费网站| 波多野结衣高清无吗| 99久久无色码亚洲精品果冻| 99国产精品99久久久久| 日韩 欧美 亚洲 中文字幕| 制服诱惑二区| 又黄又粗又硬又大视频| 日韩欧美一区二区三区在线观看| 国产成人一区二区三区免费视频网站| 搞女人的毛片| 免费在线观看完整版高清| 最新在线观看一区二区三区| 日本黄色视频三级网站网址| 最近视频中文字幕2019在线8| 三级男女做爰猛烈吃奶摸视频| 久久中文看片网| 国产单亲对白刺激| 动漫黄色视频在线观看| 欧美日韩亚洲综合一区二区三区_| 成人18禁在线播放| av天堂在线播放| 日韩大码丰满熟妇| 无限看片的www在线观看| 欧美黑人欧美精品刺激| 一卡2卡三卡四卡精品乱码亚洲| 每晚都被弄得嗷嗷叫到高潮| 一个人观看的视频www高清免费观看 | 午夜福利欧美成人| 日韩有码中文字幕| 757午夜福利合集在线观看| 精品国产美女av久久久久小说| 麻豆av在线久日| 久久人妻av系列| 女人爽到高潮嗷嗷叫在线视频| 黑人巨大精品欧美一区二区mp4| 亚洲 欧美 日韩 在线 免费| 国产成人系列免费观看| 久久久久九九精品影院| 黄色毛片三级朝国网站| АⅤ资源中文在线天堂| 97超级碰碰碰精品色视频在线观看| 熟妇人妻久久中文字幕3abv| 91麻豆精品激情在线观看国产| 国产69精品久久久久777片 | 香蕉久久夜色| 国产精品一区二区三区四区久久| 久久久久久亚洲精品国产蜜桃av| 久久婷婷成人综合色麻豆| 18禁美女被吸乳视频| 久久热在线av| 国产精品香港三级国产av潘金莲| 狂野欧美白嫩少妇大欣赏| 亚洲人成网站在线播放欧美日韩| 国产亚洲精品久久久久5区| 看片在线看免费视频| 国产精品一及| 校园春色视频在线观看| av天堂在线播放| 欧美乱妇无乱码| 变态另类成人亚洲欧美熟女| 国产亚洲精品久久久久5区| 欧美一级a爱片免费观看看 | 无人区码免费观看不卡| 亚洲精品一卡2卡三卡4卡5卡| 国产黄a三级三级三级人| 伦理电影免费视频| 国产精品久久久av美女十八| 成人av在线播放网站| 一边摸一边做爽爽视频免费| 一个人免费在线观看电影 | 国产精品 欧美亚洲| av天堂在线播放| 亚洲成av人片免费观看| 国产精品一区二区精品视频观看| 国产精品av久久久久免费| 精品久久久久久久人妻蜜臀av| 美女黄网站色视频| 熟女少妇亚洲综合色aaa.| 男女午夜视频在线观看| 欧美人与性动交α欧美精品济南到| 熟女电影av网| 久久精品91无色码中文字幕| 岛国在线观看网站| 欧美在线一区亚洲| 又黄又爽又免费观看的视频| 亚洲七黄色美女视频| 亚洲国产精品合色在线| 亚洲国产精品成人综合色| 在线观看免费视频日本深夜| 色精品久久人妻99蜜桃| 激情在线观看视频在线高清| 最近最新中文字幕大全免费视频| 中文资源天堂在线| 天天添夜夜摸| 天堂av国产一区二区熟女人妻 | 啪啪无遮挡十八禁网站| www日本在线高清视频| 免费观看人在逋| 无限看片的www在线观看| 日韩欧美免费精品| 91九色精品人成在线观看| 久久久精品欧美日韩精品| 国产亚洲精品综合一区在线观看 | 精品第一国产精品| 宅男免费午夜| 成人18禁高潮啪啪吃奶动态图| 两个人免费观看高清视频| √禁漫天堂资源中文www| 男女做爰动态图高潮gif福利片| 欧美中文综合在线视频| 一区福利在线观看| 日本三级黄在线观看| 久久精品国产亚洲av香蕉五月| www.精华液| 国产蜜桃级精品一区二区三区| 亚洲第一电影网av| 91麻豆精品激情在线观看国产| 欧美大码av| 88av欧美| 脱女人内裤的视频| 亚洲精品国产精品久久久不卡| 成人欧美大片| 久久精品综合一区二区三区| 一级毛片高清免费大全| 嫩草影院精品99| 国产午夜精品论理片| 午夜亚洲福利在线播放| www日本在线高清视频| 国产成人精品无人区| 黄色成人免费大全| 久久天躁狠狠躁夜夜2o2o| 在线a可以看的网站| 青草久久国产| 一进一出抽搐gif免费好疼| 热99re8久久精品国产| 午夜福利18| 色av中文字幕| 在线永久观看黄色视频|