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

    基于希爾伯特-黃變換的雙馬赫-曾德分布式光纖傳感振動(dòng)定位方法

    2021-11-25 08:03:26楊文晨秦增光劉兆軍徐演平叢振華王澤群
    中國(guó)光學(xué) 2021年6期
    關(guān)鍵詞:定位精度傳感探測(cè)器

    楊文晨,秦增光,劉兆軍,徐演平,李 釗,渠 帥,叢振華,王澤群

    (1.山東大學(xué) 信息科學(xué)與工程學(xué)院,山東 青島 266237;2.山東省激光技術(shù)與應(yīng)用重點(diǎn)實(shí)驗(yàn)室,山東 青島 266237;3.山東大學(xué) 激光與紅外系統(tǒng)集成技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室,山東 青島 266237)

    1 引言

    近年來(lái),光纖通信技術(shù)發(fā)展迅速并且得到廣泛應(yīng)用,光纖傳感技術(shù)也受到了廣泛的關(guān)注,在越來(lái)越多的領(lǐng)域得以應(yīng)用,大型建筑、橋梁、管道的安全監(jiān)測(cè),飛機(jī)、高鐵的運(yùn)行安全監(jiān)測(cè),電纜的溫度監(jiān)測(cè),周界安全的防護(hù)監(jiān)測(cè),海底地質(zhì)災(zāi)害預(yù)測(cè),甚至石油產(chǎn)業(yè)、電力行業(yè)、機(jī)械加工及醫(yī)療行業(yè)等領(lǐng)域都與其息息相關(guān)[1-6]。作為光纖傳感領(lǐng)域最重要的一類,分布式光纖傳感主要分為基于后向散射光的分布式光纖傳感系統(tǒng)和基于干涉儀結(jié)構(gòu)的分布式光纖傳感系統(tǒng)?;诤笙蛏⑸涞姆植际焦饫w傳感系統(tǒng)又根據(jù)后向散射光的不同,可以分為拉曼光時(shí)域反射計(jì)(Raman Optical Time Domain Reflectometer,R-OTDR)、布里淵光時(shí)域反射計(jì)(Brillouin Optical Time Domain Reflectometer,B-OTDR)和相位敏感型光時(shí)域反射計(jì)(Phasesensitive Optical Time Domain Reflectometer, φ-OTDR)?;诟缮鎯x結(jié)構(gòu)的分布式光纖傳感系統(tǒng)根據(jù)其裝置的不同可分為薩格納克(Sagnac)干涉儀、馬赫-曾德(Mach-Zehnder,M-Z)干涉儀和邁克爾遜(Michelson)干涉儀。

    基于后向散射型的分布式光纖傳感技術(shù)原理是當(dāng)光纖中的光受到外界擾動(dòng)時(shí)會(huì)產(chǎn)生相移,使經(jīng)過(guò)耦合器的光的干涉狀態(tài)發(fā)生改變。光電探測(cè)器可以把光信號(hào)轉(zhuǎn)化成電信號(hào)進(jìn)行信號(hào)采集,對(duì)采集信號(hào)進(jìn)行處理即可得到擾動(dòng)信息。然而,基于后向散射光的光纖傳感系統(tǒng)對(duì)光源、探測(cè)器等硬件要求較高,光路較為復(fù)雜,成本較高。同時(shí),受脈沖重復(fù)頻率的限制,其在長(zhǎng)距離傳感中可測(cè)量的頻率響應(yīng)范圍較低,這也極大地限制了其在實(shí)際中的應(yīng)用。

    干涉型分布式光纖傳感系統(tǒng)可以彌補(bǔ)后向散射型分布式光纖傳感系統(tǒng)的缺點(diǎn)。同時(shí),干涉型分布式光纖傳感系統(tǒng)具有頻率響應(yīng)寬、對(duì)光源要求較低、光路較為簡(jiǎn)單、成本低等優(yōu)勢(shì)。其中,基于雙M-Z結(jié)構(gòu)的分布式光纖傳感系統(tǒng)因其相位敏感性高、響應(yīng)速度快、定位精確,逐漸成為分布式光纖傳感系統(tǒng)的研究熱點(diǎn),并已在周界安全防范、石油管道安全預(yù)警、輸電線網(wǎng)安全監(jiān)控及海底線纜監(jiān)控等領(lǐng)域有廣泛的應(yīng)用[7]。國(guó)內(nèi)外有關(guān)公司已對(duì)定位系統(tǒng)產(chǎn)品進(jìn)行設(shè)計(jì)及應(yīng)用,但由于定位精度低,很難進(jìn)行大面積推廣。隨著雙MZ應(yīng)用需求的逐漸增加,對(duì)其定位精度也有了越來(lái)越高的要求。

    傳統(tǒng)的雙M-Z振動(dòng)定位方法是直接對(duì)信號(hào)做互相關(guān)運(yùn)算,通過(guò)計(jì)算時(shí)延得到振動(dòng)位置,陳偉民等使用傳統(tǒng)M-Z振動(dòng)定位方法在長(zhǎng)20 km光纜、10 MHz采樣率下進(jìn)行計(jì)算,定位誤差為149 m[8]。直接互相關(guān)計(jì)算會(huì)把大量的噪聲信息也計(jì)算其中,使得振動(dòng)定位精度下降。此外,在信號(hào)中難以確定振動(dòng)起始點(diǎn)也會(huì)對(duì)精度和處理時(shí)間帶來(lái)不利的影響[9-11]。為了提高雙M-Z振動(dòng)系統(tǒng)的定位精度,陳沁楠等采用過(guò)零率得到擾動(dòng)的起始點(diǎn),并通過(guò)互相關(guān)算法處理,在2.25 km光纖實(shí)現(xiàn)了±20 m定位精度[7]。該方法的準(zhǔn)確度過(guò)于依賴于對(duì)過(guò)零窗口大小的選擇,窗口選擇較大會(huì)引入噪聲誤差,反之則會(huì)丟失重要信息。馬春宇等采用雙波長(zhǎng)M-Z系統(tǒng)結(jié)合連續(xù)小波變換技術(shù)通過(guò)提取采集信號(hào)包絡(luò)尋找瞬時(shí)頻率較高部分以減小定位誤差,在85 km的傳感光纖上平均誤差為24.3 m[12]。謝尚然等采用離散小波尋找振動(dòng)起始點(diǎn)進(jìn)行定位,在長(zhǎng)度為30.498 km、2 MHz采樣率下,定位誤差達(dá)50 m[13]。然而由于小波分析的結(jié)果主要取決于小波變換的基本功能,無(wú)法以相同的精度對(duì)時(shí)間或頻率相關(guān)信息進(jìn)行準(zhǔn)確分析。因此尋找有效的檢測(cè)方法不僅可以提高定位精度,而且可以減少處理時(shí)間,對(duì)于干涉行儀分布式光纖傳感具有重要意義。上述研究的定位精度還需要提高,限制了其應(yīng)用。

    為了實(shí)現(xiàn)高精度的振動(dòng)定位,本文提出了一種基于希爾伯特-黃變換(Hilbert-Huang Transform,HHT)的雙M-Z型分布式光纖傳感振動(dòng)定位方法。通過(guò)HHT方法,可以準(zhǔn)確判斷瞬時(shí)高頻率信息,從而提高互相關(guān)計(jì)算的準(zhǔn)確度,減小雙M-Z傳感系統(tǒng)誤差。

    2 實(shí)驗(yàn)系統(tǒng)及原理

    雙M-Z型分布式光纖傳感實(shí)驗(yàn)系統(tǒng)裝置圖如圖1所示,由兩個(gè)M-Z干涉儀按照順時(shí)針和逆時(shí)針?lè)较蚪M合而成。光源發(fā)出的光經(jīng)隔離器到達(dá)耦合器C1分為兩束,其中一束光經(jīng)耦合器C2分為兩束后進(jìn)入傳感光纖的兩個(gè)干涉臂,在耦合器C3處干涉,并經(jīng)傳導(dǎo)光纖進(jìn)入耦合器C4被光電探測(cè)器PD2接收;另一束光經(jīng)耦合器C4進(jìn)入傳導(dǎo)光纖,在經(jīng)過(guò)耦合器C3進(jìn)入干涉儀的兩個(gè)干涉臂,兩束光在耦合器C2處發(fā)生干涉,被探測(cè)器PD1接收。兩個(gè)探測(cè)器把探測(cè)到的光信號(hào)轉(zhuǎn)為電信號(hào)后進(jìn)入采集處理系統(tǒng)進(jìn)行定位計(jì)算。

    圖1 雙M-Z分布式光纖傳感系統(tǒng)原理圖Fig.1 Schematic diagram of dual M-Z distributed optical fiber sensor

    當(dāng)有振動(dòng)作用在M點(diǎn)時(shí),攜帶振動(dòng)的信號(hào)沿光纖順時(shí)針和逆時(shí)針?lè)较虻竭_(dá)兩個(gè)探測(cè)器存在時(shí)間延遲τ:

    式中,tL?x為信號(hào)從振動(dòng)位置傳播到耦合器C3的時(shí)間,tL為振動(dòng)信號(hào)經(jīng)過(guò)傳導(dǎo)光纖到達(dá)探測(cè)器PD1的時(shí)間,tx為振動(dòng)信號(hào)從振動(dòng)位置傳到探測(cè)器PD2的時(shí)間。從而可以計(jì)算出振動(dòng)的位置:

    式中,L為傳感光纖和傳導(dǎo)光纖的長(zhǎng)度,τ為順逆兩束光到達(dá)平衡探測(cè)器間的時(shí)間延遲,c為真空中的光速,n為光纖折射率,x為振動(dòng)位置M到耦合器C1的距離。

    由于兩路探測(cè)器探測(cè)到的振動(dòng)信號(hào)是由同一振動(dòng)源產(chǎn)生,因此兩路信號(hào)具有高度的相關(guān)性,通過(guò)對(duì)兩路信號(hào)做互相關(guān)運(yùn)算可以得出時(shí)間延遲τ,進(jìn)而可以得到擾動(dòng)位置。

    分布式光纖傳感技術(shù)中空間分辨率是衡量一個(gè)系統(tǒng)性能優(yōu)良的關(guān)鍵指標(biāo)。光干涉后被兩個(gè)光電探測(cè)器接收,光信號(hào)轉(zhuǎn)化為電信號(hào),后將電信號(hào)經(jīng)過(guò)采集卡采樣之后送入計(jì)算機(jī)進(jìn)行處理?;诜菍?duì)稱雙 M-Z結(jié)構(gòu)的光纖擾動(dòng)傳感利用兩路干涉信號(hào)的時(shí)延計(jì)算擾動(dòng)位置。因此,采集卡的采樣頻率直接決定著傳感器的理論定位精度。假設(shè)采樣頻率為fs,此時(shí)相鄰采樣數(shù)據(jù)的時(shí)間間隔為?=1/fs, 任意時(shí)延均可以寫(xiě)成d=N?t的形式(N為正整數(shù))。因此可以得到理想條件下傳感器的空間分辨率為:

    由于可以認(rèn)為光纖折射率是不變的,因此定位空間分辨率只與采樣率有關(guān),采樣率越高,理論空間分辨率就越高。本文使用的采集卡的采樣率為10 MS/s,因此,本文的理論空間分辨率為10 m。

    3 希爾伯特-黃變換算法

    希爾伯特-黃變換包含經(jīng)驗(yàn)?zāi)B(tài)分解(Empirical Mode Decomposition,EMD)和希爾伯特譜分析(Hilbert Spectrum Analysis,HSA)[14-18]。

    3.1 經(jīng)驗(yàn)?zāi)B(tài)分解

    EMD的基本思想是把一個(gè)頻率變化無(wú)規(guī)律的波轉(zhuǎn)化成多個(gè)單一頻率波+殘余波的形式。具體是將一個(gè)原始信號(hào)x(t)的極大值點(diǎn)用3次Hermite差值方式擬合出上包絡(luò)線,再用同樣的方式得出下包絡(luò)線;上下包絡(luò)線取均值,得到平均包絡(luò)線,記為pl;用原始信號(hào)x(t)減去平均包絡(luò)線pl,得到一個(gè)新的信號(hào)m(t),即:

    若新信號(hào)m(t)還存在負(fù)的局部極大值和正的局部極小值,則繼續(xù)按上述操作分解,由此原始信號(hào)x(t)分解為n個(gè)本征模函數(shù)(Intrinsic Mode Function,IMF)和余波r(t),可表示為:

    IMF滿足以下兩個(gè)條件:(1) 信號(hào)極值點(diǎn)的數(shù)量與零點(diǎn)數(shù)相等或者相差1;(2) 由極大值定義的上包絡(luò)和極小值定義的下包絡(luò)的局部均值是零。

    作為燒結(jié)配料:鋼渣主要化學(xué)成分為CaO、FeO等,當(dāng)CaO含量較高時(shí),鋼渣可作燒結(jié)礦助熔劑替代部分石灰熔劑,適量配入鋼渣可改善燒結(jié)礦質(zhì)量,使轉(zhuǎn)鼓指數(shù)和燒結(jié)率提高,并且由于鋼渣中FeO的氧化放熱,降低燒結(jié)礦燃耗。

    3.2 希爾伯特譜分析

    信號(hào)x(t)的希爾伯特變換y(t)表示為:

    式中,p代表柯西主值,由此確定的解析信號(hào)為:

    式中,a(t)為瞬時(shí)幅值, θ(t)為瞬時(shí)相位,由式(8)決定。

    相位函數(shù)求導(dǎo)得到信號(hào)的瞬時(shí)頻率函數(shù):

    對(duì)所有IMF進(jìn)行Hilbert變換,得:

    這里省略了余波r(t),式中,H(f,t)是隨時(shí)間和頻率變化的信號(hào)幅度。定義Hilbert譜:

    4 實(shí)驗(yàn)過(guò)程

    搭建如圖1所示的實(shí)驗(yàn)系統(tǒng),采用1550 nm窄線寬激光器(NLL)作為傳感光源,經(jīng)過(guò)隔離器(ISO)后進(jìn)入50∶50耦合器分束,分別沿順時(shí)針和逆時(shí)針?lè)较騻鞑ィ髯园l(fā)生干涉后被光電探測(cè)器(PD)采集,進(jìn)入采集處理系統(tǒng)進(jìn)行數(shù)據(jù)分析。信號(hào)采集使用NI PCI-5114采集卡,通過(guò)LabVIEW采集信號(hào),信號(hào)處理使用Matlab進(jìn)行定位計(jì)算。采用的光纖為2 km雙芯鎧裝光纜,抗環(huán)境噪聲的能力較強(qiáng),同時(shí)更適合實(shí)際應(yīng)用。

    實(shí)驗(yàn)中,在距離C2耦合器1800 m處用腳踩光纖的方式施加振動(dòng),采樣率選用10 MHz/s,采樣時(shí)間為0.1 s,振動(dòng)時(shí)間選取0.01 s。采集到的信號(hào)如圖2(彩圖見(jiàn)期刊電子版)所示,可見(jiàn)兩路信號(hào)具有很強(qiáng)的相關(guān)性。由于光纖的雙折射因素、探測(cè)器增益帶寬和光源噪聲造成了兩路信號(hào)大小的不一致。由于兩路信號(hào)反映的都是同一振動(dòng)源的振動(dòng),因此兩路信號(hào)形狀上具有高度相關(guān)性,通過(guò)互相關(guān)計(jì)算可以得到兩路相關(guān)信號(hào)的時(shí)間延遲,這個(gè)結(jié)果主要考慮兩路信號(hào)的形狀相似度,信號(hào)幅值的大小對(duì)其影響較小。雙M-Z的定位結(jié)果與得到的時(shí)間延遲有關(guān),因此還是可以通過(guò)式(2)對(duì)振動(dòng)位置最終進(jìn)行定位。

    圖2 干涉信號(hào)時(shí)域圖Fig.2 Time domain diagram of interference signal

    該信號(hào)處理方法流程圖如圖3所示,具體運(yùn)算如下:

    圖3 定位計(jì)算方法流程圖Fig.3 Flow chart of positioning calculation method

    (1)對(duì)振動(dòng)光纖信號(hào)按分段3次Hermite差值多項(xiàng)式做EMD,得到IMF,分解層數(shù)選取10層,如圖4所示,這里給出了5層分解。

    圖4 干涉信號(hào)EMD分解圖Fig.4 EMD decomposition diagram of interference signal

    (2)對(duì)分解得到的每層IMF做Hilbert變換,結(jié)果反映了原始信號(hào)的瞬時(shí)頻率。

    (3)將所有IMF做Hilbert變換后的結(jié)果進(jìn)行疊加,得到Hilbert譜,清晰直觀地反映原始信號(hào)的頻率變化情況,如圖5所示。

    (4)根據(jù)Hilbert譜峰值位置提取原始信號(hào)的信號(hào)段,取最高峰值處附近0.01 s的數(shù)據(jù)做互相關(guān)運(yùn)算及定位運(yùn)算得出振動(dòng)位置信息。從圖5可以看出,880000采樣點(diǎn)處頻率最高,選取800000~900000采樣點(diǎn)處的區(qū)域信號(hào)作為振動(dòng)的有效信號(hào),對(duì)此做互相關(guān)運(yùn)算得出兩路信號(hào)的時(shí)間延遲,通過(guò)計(jì)算得出定位為1810 m,與實(shí)際振動(dòng)位置1800 m相比,誤差為10 m。

    圖5 干涉信號(hào)的希爾伯特譜圖Fig.5 Hilbert spectrum of interference signal

    為了說(shuō)明該方法的有效性,分別取不同頻率的區(qū)域信號(hào)做定位運(yùn)算,比較定位誤差的大小,如表1所示。通過(guò)對(duì)不同區(qū)域信號(hào)定位誤差的對(duì)比可知,本文所提出的定位方法比傳統(tǒng)的直接互相關(guān)方法定位誤差小,同時(shí)利用該方法可精確定位出最優(yōu)振動(dòng)區(qū)域信號(hào),取高頻段的信號(hào)定位誤差比低頻段信號(hào)定位誤差小。

    表1 不同區(qū)域信號(hào)定位誤差對(duì)比Tab.1 Comparison of signal positioning errors in different regions

    5 結(jié)論

    本文利用HHT來(lái)確定振動(dòng)的區(qū)域信號(hào)位置,從而提取出包含振動(dòng)的有效區(qū)域信號(hào),振動(dòng)區(qū)域?yàn)轭l率最高的信號(hào)區(qū)域,對(duì)該區(qū)域信號(hào)進(jìn)行互相關(guān)運(yùn)算,得出兩路信號(hào)的時(shí)間延遲,再進(jìn)行定位計(jì)算,減少了運(yùn)算量,并且極大地提高了定位精度,在2 km的傳感長(zhǎng)度、10 MHz采樣率下,定位精度達(dá)10 m。

    猜你喜歡
    定位精度傳感探測(cè)器
    北斗定位精度可達(dá)兩三米
    軍事文摘(2023年4期)2023-04-05 13:57:35
    《傳感技術(shù)學(xué)報(bào)》期刊征訂
    新型無(wú)酶便攜式傳感平臺(tái) 兩秒內(nèi)測(cè)出果蔬農(nóng)藥殘留
    第二章 探測(cè)器有反應(yīng)
    EN菌的引力波探測(cè)器
    GPS定位精度研究
    第二章 探測(cè)器有反應(yīng)
    IPv6與ZigBee無(wú)線傳感網(wǎng)互聯(lián)網(wǎng)關(guān)的研究
    電子制作(2018年23期)2018-12-26 01:01:26
    組合導(dǎo)航的AGV定位精度的改善
    有7顆彗星已經(jīng)被探測(cè)器造訪過(guò)
    太空探索(2014年9期)2014-07-10 13:06:30
    免费一级毛片在线播放高清视频 | 亚洲精品国产色婷婷电影| 日韩大码丰满熟妇| 你懂的网址亚洲精品在线观看| 国产欧美日韩一区二区三 | 久久午夜综合久久蜜桃| 波野结衣二区三区在线| 国产精品久久久久久精品电影小说| 亚洲成国产人片在线观看| 深夜精品福利| 亚洲五月色婷婷综合| 欧美精品一区二区免费开放| 老司机影院毛片| 天天影视国产精品| 一级黄色大片毛片| 高清不卡的av网站| 国产亚洲午夜精品一区二区久久| 国产精品欧美亚洲77777| 欧美精品一区二区免费开放| 亚洲,欧美,日韩| 在线亚洲精品国产二区图片欧美| 我要看黄色一级片免费的| 好男人视频免费观看在线| 亚洲精品国产色婷婷电影| 女人精品久久久久毛片| 亚洲精品日韩在线中文字幕| 精品久久久久久久毛片微露脸 | 色婷婷av一区二区三区视频| 美女午夜性视频免费| 欧美日韩亚洲国产一区二区在线观看 | 狠狠婷婷综合久久久久久88av| 亚洲国产欧美一区二区综合| 国产成人免费无遮挡视频| 国产亚洲精品久久久久5区| 亚洲精品第二区| 免费人妻精品一区二区三区视频| 好男人电影高清在线观看| 激情五月婷婷亚洲| 搡老岳熟女国产| 啦啦啦中文免费视频观看日本| 一区二区av电影网| 操美女的视频在线观看| 欧美人与性动交α欧美精品济南到| 国产精品国产三级专区第一集| 亚洲人成77777在线视频| 欧美在线黄色| 婷婷成人精品国产| 五月天丁香电影| 制服人妻中文乱码| 成人18禁高潮啪啪吃奶动态图| 操美女的视频在线观看| 亚洲av成人不卡在线观看播放网 | 少妇精品久久久久久久| 免费久久久久久久精品成人欧美视频| 亚洲国产欧美一区二区综合| 大陆偷拍与自拍| 十八禁网站网址无遮挡| 免费黄频网站在线观看国产| 欧美成人午夜精品| 高清欧美精品videossex| 男女国产视频网站| 一级毛片黄色毛片免费观看视频| 欧美黑人精品巨大| 自线自在国产av| 精品少妇一区二区三区视频日本电影| 国精品久久久久久国模美| 啦啦啦在线观看免费高清www| 一级,二级,三级黄色视频| 老汉色av国产亚洲站长工具| 精品一区二区三卡| 在线观看www视频免费| 老司机亚洲免费影院| 免费在线观看视频国产中文字幕亚洲 | 中文字幕亚洲精品专区| 晚上一个人看的免费电影| 美女视频免费永久观看网站| 亚洲国产精品成人久久小说| 女性被躁到高潮视频| 久久精品国产综合久久久| 亚洲精品久久午夜乱码| av又黄又爽大尺度在线免费看| 一区在线观看完整版| 日韩熟女老妇一区二区性免费视频| 久久狼人影院| 久久久久久免费高清国产稀缺| 2021少妇久久久久久久久久久| 日本午夜av视频| a级片在线免费高清观看视频| 纯流量卡能插随身wifi吗| 久久热在线av| 亚洲男人天堂网一区| 91九色精品人成在线观看| 两人在一起打扑克的视频| 亚洲五月婷婷丁香| 精品福利永久在线观看| 女人久久www免费人成看片| 国产欧美日韩一区二区三 | 成年美女黄网站色视频大全免费| 国产日韩欧美在线精品| 美女视频免费永久观看网站| 欧美精品亚洲一区二区| 亚洲欧美清纯卡通| 高清不卡的av网站| 久久青草综合色| 国产精品一区二区免费欧美 | 午夜激情久久久久久久| 在线av久久热| 国产97色在线日韩免费| 精品福利永久在线观看| 国产精品秋霞免费鲁丝片| 亚洲一区二区三区欧美精品| 波野结衣二区三区在线| 人妻一区二区av| 无遮挡黄片免费观看| 男女无遮挡免费网站观看| 久久国产精品影院| av天堂久久9| 日韩中文字幕视频在线看片| 国产国语露脸激情在线看| 肉色欧美久久久久久久蜜桃| 国产成人精品在线电影| 欧美日韩亚洲综合一区二区三区_| 1024香蕉在线观看| av有码第一页| 1024香蕉在线观看| 五月开心婷婷网| www.av在线官网国产| 亚洲国产日韩一区二区| 久久鲁丝午夜福利片| 午夜免费男女啪啪视频观看| 国产xxxxx性猛交| av在线app专区| 自线自在国产av| 成人国产一区最新在线观看 | 亚洲精品一二三| 熟女av电影| 亚洲精品一二三| 精品人妻熟女毛片av久久网站| 亚洲专区国产一区二区| 各种免费的搞黄视频| 中文字幕人妻丝袜制服| 免费高清在线观看日韩| 女人久久www免费人成看片| h视频一区二区三区| 新久久久久国产一级毛片| 91成人精品电影| 久久国产精品男人的天堂亚洲| 波多野结衣av一区二区av| 亚洲成人免费av在线播放| 国产一区有黄有色的免费视频| 色婷婷久久久亚洲欧美| 天天躁日日躁夜夜躁夜夜| 91精品伊人久久大香线蕉| 在线观看免费高清a一片| 色视频在线一区二区三区| 久久热在线av| 久久久久精品国产欧美久久久 | 成人午夜精彩视频在线观看| 亚洲视频免费观看视频| 男女床上黄色一级片免费看| 女性生殖器流出的白浆| 大话2 男鬼变身卡| 久久免费观看电影| 欧美精品一区二区免费开放| 婷婷色综合大香蕉| 一级毛片 在线播放| 又大又爽又粗| 叶爱在线成人免费视频播放| 亚洲人成网站在线观看播放| 另类亚洲欧美激情| 国产亚洲av片在线观看秒播厂| 好男人电影高清在线观看| 在线观看人妻少妇| 国产精品香港三级国产av潘金莲 | 亚洲欧美一区二区三区黑人| 久久中文字幕一级| 色精品久久人妻99蜜桃| 免费女性裸体啪啪无遮挡网站| 国产亚洲精品久久久久5区| 51午夜福利影视在线观看| a级毛片在线看网站| 一区二区av电影网| 亚洲精品国产一区二区精华液| av国产精品久久久久影院| 国产精品一区二区免费欧美 | 只有这里有精品99| 99久久精品国产亚洲精品| 中文字幕色久视频| 黄片播放在线免费| 丰满少妇做爰视频| 欧美精品亚洲一区二区| 一区二区三区激情视频| 精品福利观看| 又紧又爽又黄一区二区| 在线观看一区二区三区激情| 亚洲,欧美,日韩| 国产主播在线观看一区二区 | 五月开心婷婷网| 看十八女毛片水多多多| 满18在线观看网站| 一级,二级,三级黄色视频| 在线观看人妻少妇| 久久国产精品影院| 欧美成人精品欧美一级黄| 日本五十路高清| 热99久久久久精品小说推荐| 人妻一区二区av| 叶爱在线成人免费视频播放| 国产女主播在线喷水免费视频网站| 日韩人妻精品一区2区三区| 首页视频小说图片口味搜索 | av视频免费观看在线观看| 久久这里只有精品19| 一边亲一边摸免费视频| 中文乱码字字幕精品一区二区三区| 天天影视国产精品| 满18在线观看网站| 伦理电影免费视频| 亚洲中文字幕日韩| 宅男免费午夜| 日本av免费视频播放| 国产免费视频播放在线视频| 涩涩av久久男人的天堂| 90打野战视频偷拍视频| 国产成人免费观看mmmm| 亚洲av片天天在线观看| 在线亚洲精品国产二区图片欧美| 精品人妻熟女毛片av久久网站| 免费观看a级毛片全部| 脱女人内裤的视频| 啦啦啦啦在线视频资源| 国产不卡av网站在线观看| 国产成人影院久久av| 免费av中文字幕在线| 性色av一级| 女性生殖器流出的白浆| 欧美黄色淫秽网站| a 毛片基地| 91麻豆精品激情在线观看国产 | 色婷婷久久久亚洲欧美| 免费日韩欧美在线观看| 亚洲av电影在线观看一区二区三区| 欧美精品啪啪一区二区三区 | 在线观看免费视频网站a站| 亚洲精品成人av观看孕妇| 人妻一区二区av| 一区在线观看完整版| 成人18禁高潮啪啪吃奶动态图| av不卡在线播放| 精品视频人人做人人爽| 久久鲁丝午夜福利片| 日本午夜av视频| 亚洲av电影在线观看一区二区三区| 国产日韩欧美视频二区| 亚洲国产毛片av蜜桃av| 欧美人与善性xxx| 国产成人精品无人区| 99热国产这里只有精品6| 亚洲情色 制服丝袜| 中文字幕精品免费在线观看视频| 亚洲熟女毛片儿| 丁香六月天网| 99国产精品免费福利视频| 免费女性裸体啪啪无遮挡网站| 两性夫妻黄色片| 色综合欧美亚洲国产小说| 高清不卡的av网站| 美女视频免费永久观看网站| 国产淫语在线视频| 90打野战视频偷拍视频| 久9热在线精品视频| 国产人伦9x9x在线观看| 欧美 亚洲 国产 日韩一| 国产高清视频在线播放一区 | 熟女av电影| 女警被强在线播放| 9191精品国产免费久久| 91字幕亚洲| 叶爱在线成人免费视频播放| av又黄又爽大尺度在线免费看| 9热在线视频观看99| 午夜免费成人在线视频| 成年美女黄网站色视频大全免费| 一本—道久久a久久精品蜜桃钙片| 久久av网站| 99精国产麻豆久久婷婷| 2018国产大陆天天弄谢| 亚洲欧美精品综合一区二区三区| 两个人免费观看高清视频| 亚洲精品自拍成人| 久久国产亚洲av麻豆专区| av在线app专区| 免费观看人在逋| 国产一区亚洲一区在线观看| 亚洲国产最新在线播放| 悠悠久久av| 日韩视频在线欧美| 成人亚洲欧美一区二区av| 国产精品一国产av| 国产男女内射视频| 一区二区三区乱码不卡18| 婷婷成人精品国产| 80岁老熟妇乱子伦牲交| 纵有疾风起免费观看全集完整版| 久久久久久亚洲精品国产蜜桃av| 91精品国产国语对白视频| 国产日韩欧美视频二区| 亚洲欧美色中文字幕在线| 夫妻性生交免费视频一级片| 69精品国产乱码久久久| 男女国产视频网站| 最近手机中文字幕大全| tube8黄色片| 国产熟女欧美一区二区| 一区二区av电影网| 日日摸夜夜添夜夜爱| 精品少妇一区二区三区视频日本电影| 国产免费又黄又爽又色| 成人免费观看视频高清| 欧美日韩视频精品一区| svipshipincom国产片| 美国免费a级毛片| 乱人伦中国视频| 美女高潮到喷水免费观看| 国产成人影院久久av| 日本猛色少妇xxxxx猛交久久| 大片免费播放器 马上看| 女人被躁到高潮嗷嗷叫费观| 国产亚洲av片在线观看秒播厂| 男女床上黄色一级片免费看| 色精品久久人妻99蜜桃| 久久久久国产一级毛片高清牌| 大片免费播放器 马上看| 超碰97精品在线观看| 亚洲成人国产一区在线观看 | 亚洲精品久久午夜乱码| 欧美成人午夜精品| 日韩 欧美 亚洲 中文字幕| 国产有黄有色有爽视频| 国产极品粉嫩免费观看在线| 久久热在线av| 日韩一区二区三区影片| 欧美黑人精品巨大| 在线看a的网站| 中文字幕色久视频| 99九九在线精品视频| www.自偷自拍.com| 成人免费观看视频高清| 一区在线观看完整版| 久久午夜综合久久蜜桃| 日韩 亚洲 欧美在线| 国产又色又爽无遮挡免| 视频区图区小说| 亚洲欧美一区二区三区国产| 热99久久久久精品小说推荐| 午夜老司机福利片| 欧美精品亚洲一区二区| 亚洲av日韩精品久久久久久密 | 免费日韩欧美在线观看| 亚洲精品一二三| 又大又黄又爽视频免费| 亚洲午夜精品一区,二区,三区| 一二三四社区在线视频社区8| 日本vs欧美在线观看视频| 一区二区三区激情视频| 欧美日韩综合久久久久久| 一级毛片电影观看| 国产爽快片一区二区三区| 久久免费观看电影| 高清不卡的av网站| 国产免费又黄又爽又色| 欧美亚洲 丝袜 人妻 在线| 丁香六月欧美| 女人被躁到高潮嗷嗷叫费观| 亚洲欧美一区二区三区黑人| 久久久久精品国产欧美久久久 | 极品人妻少妇av视频| 欧美乱码精品一区二区三区| 男女床上黄色一级片免费看| 操出白浆在线播放| 王馨瑶露胸无遮挡在线观看| 在线观看免费视频网站a站| 国产女主播在线喷水免费视频网站| 久久人人爽av亚洲精品天堂| 黄色 视频免费看| 国产一区二区激情短视频 | 亚洲天堂av无毛| 老司机在亚洲福利影院| 日本黄色日本黄色录像| 99国产精品一区二区蜜桃av | 一级毛片 在线播放| 亚洲精品中文字幕在线视频| 黑丝袜美女国产一区| 一级毛片电影观看| 日韩制服丝袜自拍偷拍| 亚洲一区二区三区欧美精品| 亚洲情色 制服丝袜| 好男人电影高清在线观看| 男女下面插进去视频免费观看| 两个人看的免费小视频| 精品人妻1区二区| 男女边摸边吃奶| 热re99久久国产66热| h视频一区二区三区| 99精品久久久久人妻精品| 老司机影院成人| 亚洲久久久国产精品| 亚洲精品乱久久久久久| 蜜桃国产av成人99| 男的添女的下面高潮视频| 激情视频va一区二区三区| 国产三级黄色录像| 日韩av免费高清视频| 亚洲av在线观看美女高潮| 两性夫妻黄色片| 丝袜脚勾引网站| 最近中文字幕2019免费版| 青青草视频在线视频观看| 国产一区有黄有色的免费视频| 日本色播在线视频| 午夜激情av网站| 51午夜福利影视在线观看| 亚洲成人手机| 日本黄色日本黄色录像| 国产高清videossex| 国产xxxxx性猛交| 日本vs欧美在线观看视频| 纯流量卡能插随身wifi吗| 婷婷色麻豆天堂久久| 丰满少妇做爰视频| 日韩,欧美,国产一区二区三区| 久久久精品94久久精品| videos熟女内射| 欧美日韩福利视频一区二区| 国产亚洲av高清不卡| 国产免费福利视频在线观看| 丝袜美足系列| 亚洲av日韩精品久久久久久密 | 亚洲欧美日韩高清在线视频 | 中文字幕高清在线视频| 欧美人与善性xxx| 欧美少妇被猛烈插入视频| 黑人巨大精品欧美一区二区蜜桃| 精品福利永久在线观看| 少妇猛男粗大的猛烈进出视频| 久久久欧美国产精品| 国产精品亚洲av一区麻豆| 亚洲第一青青草原| 一区二区三区精品91| 男男h啪啪无遮挡| 欧美人与善性xxx| 国产精品一区二区在线不卡| 国产在线观看jvid| 精品久久久久久电影网| 免费看av在线观看网站| 人人妻人人澡人人爽人人夜夜| 亚洲国产精品一区三区| 国产日韩欧美在线精品| 国产一区二区在线观看av| 老汉色∧v一级毛片| 1024视频免费在线观看| 国产精品久久久久成人av| 国产一卡二卡三卡精品| 国产成人啪精品午夜网站| 亚洲伊人色综图| 午夜av观看不卡| 日韩 亚洲 欧美在线| 黑丝袜美女国产一区| 亚洲欧洲精品一区二区精品久久久| 午夜福利视频在线观看免费| 国产精品秋霞免费鲁丝片| 欧美亚洲日本最大视频资源| 99国产精品一区二区蜜桃av | 亚洲欧美一区二区三区久久| 久久精品国产亚洲av涩爱| 欧美大码av| 赤兔流量卡办理| 狂野欧美激情性xxxx| 一本综合久久免费| 男人爽女人下面视频在线观看| 精品少妇黑人巨大在线播放| 国产亚洲欧美精品永久| 午夜视频精品福利| 中文欧美无线码| 五月天丁香电影| 亚洲国产精品999| 别揉我奶头~嗯~啊~动态视频 | 中文精品一卡2卡3卡4更新| 欧美日韩黄片免| 天天躁夜夜躁狠狠久久av| 老司机深夜福利视频在线观看 | 久久久国产精品麻豆| 亚洲欧美日韩高清在线视频 | 色精品久久人妻99蜜桃| 国产熟女午夜一区二区三区| 免费日韩欧美在线观看| 日韩,欧美,国产一区二区三区| 国产亚洲午夜精品一区二区久久| 极品少妇高潮喷水抽搐| 欧美人与性动交α欧美精品济南到| 水蜜桃什么品种好| 伊人亚洲综合成人网| 一级毛片我不卡| 国产亚洲欧美在线一区二区| 亚洲国产精品一区三区| 亚洲专区中文字幕在线| 久久99精品国语久久久| 2021少妇久久久久久久久久久| 色网站视频免费| 人人妻人人添人人爽欧美一区卜| 国产精品香港三级国产av潘金莲 | 一区二区三区乱码不卡18| 午夜福利,免费看| 国产一区亚洲一区在线观看| 中文字幕制服av| 1024视频免费在线观看| 一区二区三区乱码不卡18| 男女床上黄色一级片免费看| av天堂久久9| 午夜免费男女啪啪视频观看| 一边亲一边摸免费视频| 日韩欧美一区视频在线观看| 91麻豆av在线| 国产精品麻豆人妻色哟哟久久| 在线av久久热| 精品人妻1区二区| 亚洲伊人久久精品综合| 韩国精品一区二区三区| 男女免费视频国产| 制服诱惑二区| 国产亚洲一区二区精品| 日韩av在线免费看完整版不卡| 无限看片的www在线观看| 国产av精品麻豆| 王馨瑶露胸无遮挡在线观看| 一级毛片黄色毛片免费观看视频| 一本久久精品| 大型av网站在线播放| 成年人午夜在线观看视频| 肉色欧美久久久久久久蜜桃| 午夜91福利影院| 看十八女毛片水多多多| 中文字幕av电影在线播放| 亚洲av美国av| 久久综合国产亚洲精品| 91麻豆精品激情在线观看国产 | 19禁男女啪啪无遮挡网站| 亚洲欧洲精品一区二区精品久久久| 日韩欧美一区视频在线观看| 亚洲综合色网址| 两个人免费观看高清视频| 伦理电影免费视频| 国产不卡av网站在线观看| 伦理电影免费视频| 精品人妻在线不人妻| 在线观看国产h片| 日韩免费高清中文字幕av| 精品人妻熟女毛片av久久网站| 一级毛片女人18水好多 | 午夜老司机福利片| 涩涩av久久男人的天堂| 亚洲精品乱久久久久久| 精品视频人人做人人爽| 亚洲国产中文字幕在线视频| av视频免费观看在线观看| 国产片内射在线| 高清视频免费观看一区二区| 人人妻人人添人人爽欧美一区卜| 国产一区二区三区综合在线观看| 黄色毛片三级朝国网站| 国产免费现黄频在线看| 少妇猛男粗大的猛烈进出视频| 日韩熟女老妇一区二区性免费视频| 丝袜美足系列| 国产精品av久久久久免费| 咕卡用的链子| 欧美激情高清一区二区三区| 欧美性长视频在线观看| 午夜福利免费观看在线| 99国产精品99久久久久| 色网站视频免费| 五月天丁香电影| 欧美精品人与动牲交sv欧美| 久热爱精品视频在线9| 只有这里有精品99| 晚上一个人看的免费电影| 日本一区二区免费在线视频| www.自偷自拍.com| 亚洲免费av在线视频| 亚洲国产看品久久| 国产视频首页在线观看| 在线精品无人区一区二区三| 黄色毛片三级朝国网站| 香蕉丝袜av| 国产黄色视频一区二区在线观看| 欧美精品一区二区免费开放| 亚洲 国产 在线| 日韩中文字幕视频在线看片| 国产av一区二区精品久久| 国产精品免费大片| 午夜免费男女啪啪视频观看| 91老司机精品| 国产免费现黄频在线看| 中文字幕精品免费在线观看视频| 亚洲av日韩在线播放| 亚洲精品第二区| 最黄视频免费看| 日本91视频免费播放| 国产精品欧美亚洲77777| 黄片小视频在线播放| 91精品伊人久久大香线蕉| 日本a在线网址| 一二三四在线观看免费中文在| 一级片免费观看大全| 热re99久久精品国产66热6|