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

    場道鋼筋加固區(qū)的災(zāi)害目標(biāo)檢測

    2018-08-20 06:16:50曹蕓茜吳仁彪何煒琨
    信號處理 2018年2期
    關(guān)鍵詞:時變分段災(zāi)害

    曹蕓茜 吳仁彪 何煒琨

    (1. 中國民航大學(xué)智能信號與圖像處理天津市重點實驗室,天津 300300; 2. 中國民航大學(xué)基礎(chǔ)實驗中心,天津 300300)

    1 引言

    機場作為城市和航空運輸?shù)幕A(chǔ)設(shè)施,是交通運輸?shù)闹匾M成部分。機場道面是機場建設(shè)的核心,跑道質(zhì)量的好壞將會直接影響機場與飛機運行的安全,其中場道內(nèi)脫空和裂縫等隱性災(zāi)害的檢測和維護(hù)至關(guān)重要[1]。在不影響機場正常運行且不破壞道面結(jié)構(gòu)的基礎(chǔ)上,對道面進(jìn)行無損檢測是必然趨勢。探地雷達(dá)(Ground Penetrating Radar, GPR)是一種有效的無損檢測工具,其具有探測速度快、操作方便、分辨率高等特點[2-3]。在機場場道的建設(shè)施工中,為了加強其承載能力,需要鋪設(shè)鋼筋進(jìn)行加固[4]。鋼筋的強反射回波甚至多次反射后的回波都會淹沒災(zāi)害目標(biāo)的回波,所以有效的抑制鋼筋強反射回波對災(zāi)害目標(biāo)的檢測是十分關(guān)鍵和必要的。

    本文首先利用非一致性檢測對探地雷達(dá)接收數(shù)據(jù)進(jìn)行預(yù)處理,利用時變增益函數(shù)對鋼筋反射一次回波進(jìn)行粗抑制,然后對該數(shù)據(jù)在時間維進(jìn)行分段均衡,以濾除鋼筋在機場場道下各層分界面的多次反射回波,再利用時頻—時空域二重S變換方法對災(zāi)害目標(biāo)進(jìn)行檢測,并同時顯示其在場道下的空間位置。該方法能在災(zāi)害目標(biāo)回波信號被鋼筋多次強反射回波信號淹沒的情況下,有效抑制鋼筋對災(zāi)害目標(biāo)的影響,達(dá)到檢測道面下微小災(zāi)害目標(biāo)空間位置的目的,為機場養(yǎng)護(hù)施工提供直觀依據(jù)。

    2 信號模型

    探地雷達(dá)接收信號模型可以表示為:

    w(t)=o(t)+d(t)+s(t)

    (1)

    如圖1所示,o(t)表示收發(fā)天線直接耦合波,在收發(fā)天線間距確定后直接耦合波分量一般相對固定。d(t)表示各層媒質(zhì)表面反射回波,由于探測環(huán)境往往是變化的,所以d(t)的變化性和不可預(yù)測性較大。s(t)表示媒質(zhì)內(nèi)目標(biāo)反射回波,也就是我們需要的回波分量。

    圖1 GPR信號模型

    我們將目標(biāo)反射回波s(t)進(jìn)行采樣,那么二維B-Scan數(shù)據(jù)中,雷達(dá)測試線l=1,2,…,L的第i=1,2,…,T個采樣點的目標(biāo)反射回波信號可以表示為sl(i)。

    3 鋼筋回波抑制

    本節(jié)研究鋼筋的強反射回波及多次反射后回波的抑制方法。首先利用基于廣義內(nèi)積的非一致性檢測對探地雷達(dá)接收數(shù)據(jù)進(jìn)行預(yù)處理,然后利用時變增益函數(shù)對鋼筋的一次強反射回波進(jìn)行粗抑制,最后利用時間域的分段均衡算法抑制鋼筋的多次反射回波。

    3.1 預(yù)處理

    我們將探地雷達(dá)接收到的回波信號w(t)進(jìn)行預(yù)處理,去除直接耦合波o(t)和各層媒質(zhì)表面反射回波d(t),以提取出鋼筋及災(zāi)害目標(biāo)的反射回波s(t)。

    首先去除探地雷達(dá)接收到的回波信號中收發(fā)天線間的直接耦合波o(t)。我們可以將雷達(dá)天線對準(zhǔn)開闊無物體的空間,或?qū)⑻炀€置于暗室內(nèi),錄取對空信號,此信號可以近似為直接耦合信號,再將雷達(dá)接收回波減去此直接耦合信號即有效去除了直接耦合波。

    然后去除機場場道下各層媒質(zhì)表面反射回波d(t)。由于探測環(huán)境的變化,d(t)的變化性和不可預(yù)知性較大,所以,我們需要根據(jù)探測區(qū)域的數(shù)據(jù)對d(t)進(jìn)行有效的去除。這里,我們采用基于廣義內(nèi)積的非一致性檢測方法合理去除d(t)。

    基于廣義內(nèi)積的非一致性檢測算法,其輸出定義為:

    (2)

    其中l(wèi)=1,2,…,L為雷達(dá)測試線數(shù),xl為去除直接耦合波后的雷達(dá)采樣數(shù)據(jù)的矩陣表示,Rl為自相關(guān)矩陣,其可以表示為:

    (3)

    計算出雷達(dá)所有探測點的gl(GIP)后,給出閾值,大于此閾值則說明與周圍差異較大,則可認(rèn)為此部分探測區(qū)域存在目標(biāo)反射信號,我們將此類區(qū)域數(shù)據(jù)濾除,剩下的小于閾值部分的回波即為各層媒質(zhì)表面強反射回波,選取該回波信號作為背景回波,將雷達(dá)接收回波減去背景回波即可以有效去除各層媒質(zhì)表面反射回波而有效的保留了目標(biāo)反射回波。

    3.2 鋼筋回波抑制

    在機場跑道使用過程中,初期常見的隱性災(zāi)害主要有脫空和裂縫[5- 6]。脫空主要出現(xiàn)在面層和上基層的交界處,裂縫主要出現(xiàn)在面層,而加固的鋼筋主要位于面層的上部,如圖2所示,所以鋼筋的強反射回波會嚴(yán)重影響甚至淹沒脫空和裂縫的反射回波,并且鋼筋與場道各層介面的多次反射回波也會嚴(yán)重影響場道災(zāi)害目標(biāo)的檢測,所以有效地抑制鋼筋的強反射回波對災(zāi)害目標(biāo)進(jìn)行檢測是十分關(guān)鍵和必要的。

    圖2 鋼筋及災(zāi)害目標(biāo)分布模型

    由于鋼筋的一次反射回波遠(yuǎn)遠(yuǎn)強于多次反射回波,且多次反射回波往往會與災(zāi)害目標(biāo)反射回波混疊,不易處理,所以我們將鋼筋回波抑制分為兩步,第一步利用時變增益函數(shù)對鋼筋一次反射回波進(jìn)行粗抑制,第二步利用時域分段均衡方法對鋼筋的多次反射回波進(jìn)行細(xì)抑制,以便提取出有效的災(zāi)害目標(biāo)反射回波。

    3.2.1粗抑制—時變增益

    機場場道在進(jìn)行鋼筋加固時有嚴(yán)格的規(guī)范并且有先驗的鋼筋位置信息,我們可以根據(jù)鋼筋位置的先驗知識,構(gòu)造合適的時變增益函數(shù)Pl(t),其中l(wèi)=1,2,…,L為雷達(dá)測試線,將其采樣后記為Pl(i),那么時變增益后的雷達(dá)數(shù)據(jù)為

    ul(i)=Pl(i)·sl(i),i=1,2,…,T,l=1,2,…,L

    (4)

    時變增益函數(shù)可以抑制鋼筋一次強反射回波時間窗內(nèi)數(shù)據(jù),增強需要觀察的機場場道面層及基層時間窗內(nèi)數(shù)據(jù),抑制不需要觀察的場道底部數(shù)據(jù),即選擇需要觀察的數(shù)據(jù)予以增強,不需要觀察的數(shù)據(jù)予以抑制,對鋼筋的一次反射回波進(jìn)行粗抑制,此方法可以有效抑制鋼筋一次反射回波,而增強了災(zāi)害目標(biāo)反射回波,但是鋼筋的多次反射回波也同災(zāi)害目標(biāo)反射回波一同被增強,所以我們還要對鋼筋的多次反射回波進(jìn)一步抑制。

    3.2.2細(xì)抑制—分段均衡

    如果某一深度存在目標(biāo),則此深度目標(biāo)所在位置必與同一深度其他位置的回波有差異且鋼筋對各層表面多次反射回波呈現(xiàn)較強的規(guī)律性,所以本文提出分段均衡的方法,對B-Scan數(shù)據(jù)進(jìn)行時間維分段,再利用數(shù)據(jù)信息求出分段矢量,對每一段數(shù)據(jù)進(jìn)行均衡,這樣可以對相對規(guī)律的數(shù)據(jù)進(jìn)行均衡濾波,對存在局部差異的數(shù)據(jù)起到增強的效果,有利于細(xì)節(jié)的辨認(rèn)。本文利用分段均衡的方法來抑制鋼筋的多次反射回波。

    首先計算鋼筋粗抑制后的數(shù)據(jù)中所有雷達(dá)測試線l=1,2,…,L在第i個采樣點的平均值:

    (5)

    根據(jù)災(zāi)害目標(biāo)的厚度合理選取分段數(shù)據(jù)點數(shù)m=T/K,其中K為分段數(shù)目,T為總時間采樣點數(shù)。局部搜索第k=1,2,…,K段的最大值maxk及最小值mink,定義差值矢量為

    d=[max1-min1max2-min2…maxK-minK]

    (6)

    為了方便將上式記為

    d=[d1d2…dk]

    (7)

    設(shè)dmax為d中最大值,定義均衡放大矢量為

    a=[d1/dmaxd2/dmax…dK/dmax]

    (8)

    記為a=[a1a2…ak]。

    鋼筋粗抑制后的B-Scan接收數(shù)據(jù)用矩陣可以表示為u。根據(jù)上述分段方法,按照相同原理對u在時間維進(jìn)行同樣的分段。設(shè)uk為對u進(jìn)行分段的第k段數(shù)據(jù),則uk可以表示為:

    uk={ul(i),i=(k-1)m+1,(k-1)m+2,…,km;

    l=1,2,…,L;k=1,2,…,K}

    (9)

    利用每一段數(shù)據(jù)的均衡放大系數(shù)對相應(yīng)的該段數(shù)據(jù)進(jìn)行加權(quán)均衡,得到第k段數(shù)據(jù)的均衡處理結(jié)果:

    yk=akuk

    (10)

    其中,k=1,2,…,K。那么,B-Scan接收數(shù)據(jù)的均衡處理結(jié)果可以表示為:

    Y=[y1y2…yK]T

    (11)

    4 災(zāi)害目標(biāo)檢測

    鋼筋回波抑制后的數(shù)據(jù)在災(zāi)害目標(biāo)附近可能還會存在少許干擾回波,這也會干擾體積或厚度較小的災(zāi)害目標(biāo),譬如小裂縫、脫空薄層的檢測。本文提出了時頻—時空域二重S變換災(zāi)害目標(biāo)檢測方法,可以有效的檢測體積及厚度較小的災(zāi)害目標(biāo)。

    S變換[7-9]是基于短時窗傅里葉變換和小波變換提出來的一種時頻分析方法,該方法吸收了兩者的優(yōu)點,同時克服了它們的缺點。

    設(shè)函數(shù)h(t)∈L2(R)(L2(R)表示能量有限的函數(shù)空間)的S變換的定義如下:

    (12)

    S變換的基本小波定義為:

    (13)

    (14)

    由式(14)可以看出,S變換采用寬度隨頻率呈反比變化的高斯窗函數(shù),低頻段的時窗較寬,可以獲得較高的頻率分辨率,高頻段的時窗較窄,可獲得較高的時間分辨率[10]。利用高頻段時間分辨率較高的特性可以有效的檢測出體積及厚度較小的災(zāi)害目標(biāo)。

    首先選取均衡后數(shù)據(jù)Y中的A-Scan數(shù)據(jù)yl(t),雷達(dá)測試線l∈[1,L]。對yl(t)進(jìn)行時間—頻率維S變換:

    (15)

    根據(jù)S變換結(jié)果選取時間分辨率較高的頻率段f0…fN,固定該頻率段中的一個頻率點fn(n∈[0,N]),在該頻率點下對yl(t)再次進(jìn)行S變換:

    (16)

    對B-Scan數(shù)據(jù)Y中的所有測試線l=1,2,…,L的數(shù)據(jù)yl(t)分別進(jìn)行S變換,即可得到一個固定頻點fn的時間—距離維S變換結(jié)果:

    Sfn(τ,l)={Sl(τ,fn),l=1,…,L}

    (17)

    將選取的頻率段f0…fN中所有頻點的時間—距離維S變換數(shù)據(jù)疊加,進(jìn)行頻域均衡濾波:

    (18)

    根據(jù)疊加結(jié)果不但可以有效檢測出體積及厚度較小的災(zāi)害目標(biāo),還可以得到災(zāi)害目標(biāo)的空間位置。

    綜上所述,機場場道鋼筋加固區(qū)的災(zāi)害目標(biāo)檢測方法的流程如圖3所示。

    圖3 災(zāi)害目標(biāo)檢測方法流程

    5 仿真結(jié)果與分析

    為了驗證本文方法的有效性,實驗中所有數(shù)據(jù)均由GPRMAX2.0[11]高保真仿真軟件產(chǎn)生,采用雙基地探地雷達(dá),選取天線間距為0.1 m,距地面高度為0.04 m,發(fā)射信號采用ricker波脈沖形式,中心頻率為900 MHz,采集數(shù)據(jù)時間窗為25 ns,采集測試線間距為0.02 m,共采集118道數(shù)據(jù)。仿真模型為一般機場場道基本四層模型,層厚從上到下依次為0.34 m,0.20 m,0.20 m,0.26 m,介電常數(shù)依次為9,12,15,22。

    圖4為本實驗中所采用的機場道面仿真數(shù)據(jù)原始幾何模型圖,該模型中含有五根鋼筋,一條裂縫和一個脫空,鋼筋均勻分布在面層中距地表面17 cm,每條鋼筋間距為50 cm,裂縫長18 cm,寬5 mm,脫空寬20 cm,厚7 mm。

    圖4 原始數(shù)據(jù)幾何模型

    圖5為探地雷達(dá)接收到的二維B-Scan回波數(shù)據(jù)。由于地表反射回波很強,淹沒了地下目標(biāo)的回波信息,所以首先必須有效的去除地表反射回波,提取需要的回波信息。圖6為利用非一致性檢測濾除背景回波后的目標(biāo)反射回波。背景的四層媒質(zhì)表面反射回波被有效濾除,凸顯出地下目標(biāo)的反射回波,但是由于場道淺層鋼筋反射回波很強,且存在鋼筋多次反射回波,嚴(yán)重影響甚至淹沒鋼筋下災(zāi)害目標(biāo)的回波,無法有效地檢測出災(zāi)害目標(biāo),尤其是細(xì)小的裂縫。

    圖6 目標(biāo)反射回波數(shù)據(jù)

    根據(jù)圖6所示數(shù)據(jù),選擇合適的時變增益時窗,對強金屬回波及不需要觀察的場道深層處數(shù)據(jù)進(jìn)行抑制,對需要觀察的數(shù)據(jù)段進(jìn)行增強,如圖7所示。圖8為鋼筋回波粗抑制后B-Scan數(shù)據(jù),脫空的反射回波被增強,可以觀察到,但是由于裂縫上下表面比較窄,反射回波很弱,并且鋼筋與地下各層媒質(zhì)表面的多次反射回波也比較強,所以裂縫還是不能有效觀測。圖9為分段均衡后B-Scan數(shù)據(jù),鋼筋在各層媒質(zhì)間的多次反射回波被有效均衡,突出了脫空的反射回波,但仍無法有效檢測出災(zāi)害目標(biāo)。

    圖7 時變增益

    圖8 鋼筋回波粗抑制

    圖9 鋼筋回波細(xì)抑制

    圖10為均衡濾波后采集數(shù)據(jù)進(jìn)行時間—頻率維S變換的時頻分析結(jié)果,我們選取其中時間分辨率較高的高頻段2500~2900 MHz數(shù)據(jù),先確定該頻段內(nèi)的一個頻點,對均衡后B-Scan數(shù)據(jù)進(jìn)行時間—距離維S變換,將所有頻段內(nèi)的S變換結(jié)果進(jìn)行頻域均衡濾波,得到如圖11所示結(jié)果,可以明顯檢測出脫空及裂縫的存在。圖12為不進(jìn)行均衡濾波而直接用相同的方法進(jìn)行S變換的結(jié)果,雖然能夠檢測到脫空,但是不能有效地檢測到裂縫的存在,且鋼筋的多次反射回波嚴(yán)重影響了災(zāi)害目標(biāo)的檢測。

    圖10 S變換時頻分析

    圖11 災(zāi)害目標(biāo)檢測結(jié)果

    圖12 未分段均衡檢測結(jié)果

    6 結(jié)論

    考慮到鋼筋的強反射回波對場道鋼筋加固區(qū)災(zāi)害目標(biāo)檢測的影響,本文首先提出時變增益函數(shù)和分段均衡的方法對鋼筋的多次反射回波進(jìn)行抑制,然后提出時頻—時空域二重S變換方法對災(zāi)害目標(biāo)進(jìn)行檢測。仿真結(jié)果表明本文算法可以有效去除鋼筋的多次反射回波對厚度及體積較小的常見災(zāi)害目標(biāo)回波的影響,對災(zāi)害目標(biāo)進(jìn)行檢測并顯示其在場道下的空間位置,為機場場道的養(yǎng)護(hù)施工提供直觀的依據(jù)。

    [1] 翁興中,蔡良才.機場道面設(shè)計[M].北京:人民交通出版社,2007.

    Weng Xingzhong, Cai Liangcai. Design of Airport Runway[M].Beijing: China Communications Press, 2007.(in Chinese)

    [2] 何煒琨, 吳仁彪, 劉家學(xué). 基于探地雷達(dá)的機場場道脫空層檢測及厚度估計[J]. 信號處理, 2011, 27(10):1547-1551.

    He Weikun, Wu Renbiao, Liu Jiaxue. Void-layer Detection and Depth Determination in Runways based on GPR[J]. Signal Processing, 2011, 27(10):1547-1551.(in Chinese)

    [3] Hoarau Q, Ginolhac G, Atto A M, et al. Robust Adaptive Detection of Buried Pipes using GPR[C]∥24th European Signal Processing Conference, 2016:533-537.

    [4] Harry M J. Ground penetrating radar theory and applications[M]. London: Elsevier, 2009.

    [5] 何煒琨. 基于探地雷達(dá)的機場場道質(zhì)量監(jiān)測關(guān)鍵技術(shù)研究[D]. 天津:天津大學(xué), 2012.

    He Weikun. Enabling Techniques for Runways Quality Surveillance via Ground Penetrating Radar[D]. Tianjin: The University of Tianjin, 2012.(in Chinese)

    [6] 龐希斌, 徐進(jìn), 盧小賓, 等. 地質(zhì)雷達(dá)在機場跑道缺陷檢測中的應(yīng)用——以“5·12”汶川地震后九寨黃龍機場檢測為例[J]. 西南民族大學(xué)學(xué)報, 2008, 34(6):1096-1100.

    Pang Xibin, Xu Jin, Lu Xiaobin, et al. Application of the ground penetrating radar to the inspection of the airfield runway defect——Taking the detection of Jiuhuang Airport after the 5.12 Wenchuan earthquake for example[J]. Journal of Southwest University for Nationalities Natural Science Edition, 2008, 34(6):1096-1100.(in Chinese)

    [7] Yan Jin, Liu Jie. Parameter estimation of frequency hopping signals based on the Robust S-transform algorithms in alpha stable noise environment[J]. International Journal of Electronics & Communications, 2016, 70(5): 611- 616.

    [8] Kazemi K, Amirian M, Dehghani M J. The S-transform using a new window to improve frequency and time resolutions[J]. Signal, Image and Video Processing, 2014, 8(3): 533-541.

    [9] 劉志成, 王殿偉. 非線性調(diào)頻信號的自適應(yīng)時頻濾波算法[J]. 信號處理, 2015, 31(3):356-363.

    Liu Zhicheng, Wang Dianwei.Adaptive Time-Frequency Filter Method of Nonlinear Frequency Modulation Signal[J]. Journal of Signal Processing, 2015, 31(3):356-363.(in Chinese)

    [10]Huang Nantian, Peng Hua, Cai Guowei, et al. Power Quality Disturbances Feature Selection and Recognition Using Optimal Multi-Resolution Fast S-Transform and CART Algorithm[J]. Energies, 2016, 9(11): 1-21.

    [11]GprMax2.0[OL]. http:∥www. gprmax.org.

    猜你喜歡
    時變分段災(zāi)害
    河南鄭州“7·20”特大暴雨災(zāi)害的警示及應(yīng)對
    一類連續(xù)和不連續(xù)分段線性系統(tǒng)的周期解研究
    我省汛期常見氣象災(zāi)害及防御
    推動災(zāi)害防治工作實現(xiàn)新跨越
    分段計算時間
    基于時變Copula的股票市場相關(guān)性分析
    智富時代(2017年4期)2017-04-27 17:08:47
    3米2分段大力士“大”在哪兒?
    太空探索(2016年9期)2016-07-12 10:00:04
    煙氣輪機復(fù)合故障時變退化特征提取
    基于MEP法的在役橋梁時變可靠度研究
    關(guān)于年齡分段的描述
    91av网站免费观看| 久久久久久久久中文| 国产精品影院久久| 别揉我奶头~嗯~啊~动态视频| 免费在线观看影片大全网站| 三级毛片av免费| 人人妻人人爽人人添夜夜欢视频| 久久久久久亚洲精品国产蜜桃av| 男女下面进入的视频免费午夜 | 国产成人系列免费观看| 热re99久久国产66热| 真人一进一出gif抽搐免费| 久久久久久久午夜电影 | 国产成人精品在线电影| 999久久久精品免费观看国产| 窝窝影院91人妻| 女警被强在线播放| 亚洲精品中文字幕一二三四区| 99re在线观看精品视频| 成人永久免费在线观看视频| 免费在线观看亚洲国产| 日韩视频一区二区在线观看| 精品国产一区二区三区四区第35| 亚洲色图av天堂| 手机成人av网站| 欧美一级毛片孕妇| 嫁个100分男人电影在线观看| 午夜福利在线免费观看网站| 丰满的人妻完整版| 精品电影一区二区在线| 手机成人av网站| 99精品在免费线老司机午夜| 国产精品 欧美亚洲| 色婷婷av一区二区三区视频| 黑丝袜美女国产一区| 成人三级黄色视频| 亚洲中文字幕日韩| 91在线观看av| 国产高清国产精品国产三级| 无人区码免费观看不卡| 99国产精品一区二区蜜桃av| 国产一区二区在线av高清观看| 国产精品久久久人人做人人爽| av天堂在线播放| 国产精品98久久久久久宅男小说| 久久久久久久久中文| 国产在线精品亚洲第一网站| 女警被强在线播放| 丝袜人妻中文字幕| 88av欧美| 91字幕亚洲| 9热在线视频观看99| 高清在线国产一区| 欧美人与性动交α欧美精品济南到| 久久久水蜜桃国产精品网| 国产野战对白在线观看| 亚洲av五月六月丁香网| av网站在线播放免费| 久久精品成人免费网站| 国产一区二区三区综合在线观看| 三上悠亚av全集在线观看| 午夜免费观看网址| 在线十欧美十亚洲十日本专区| 熟女少妇亚洲综合色aaa.| 无人区码免费观看不卡| 久久九九热精品免费| 人成视频在线观看免费观看| 美女高潮喷水抽搐中文字幕| 涩涩av久久男人的天堂| 亚洲精品在线美女| av国产精品久久久久影院| 国产单亲对白刺激| 精品无人区乱码1区二区| 51午夜福利影视在线观看| 精品一品国产午夜福利视频| 男女午夜视频在线观看| 午夜福利影视在线免费观看| 午夜成年电影在线免费观看| 成人国语在线视频| 制服诱惑二区| 天天添夜夜摸| 亚洲av美国av| 搡老岳熟女国产| 好男人电影高清在线观看| 757午夜福利合集在线观看| 精品国产乱子伦一区二区三区| 日韩欧美在线二视频| 免费看十八禁软件| 日韩人妻精品一区2区三区| 黑丝袜美女国产一区| 亚洲,欧美精品.| 久热这里只有精品99| 国产男靠女视频免费网站| 日本黄色视频三级网站网址| 变态另类成人亚洲欧美熟女 | 在线观看免费日韩欧美大片| 日韩视频一区二区在线观看| 国产一区二区在线av高清观看| 男女做爰动态图高潮gif福利片 | 国产成人一区二区三区免费视频网站| 夜夜爽天天搞| 淫秽高清视频在线观看| 黄色女人牲交| 久久人人精品亚洲av| 亚洲成人免费电影在线观看| 免费在线观看视频国产中文字幕亚洲| 最近最新中文字幕大全电影3 | 日日摸夜夜添夜夜添小说| 韩国av一区二区三区四区| 法律面前人人平等表现在哪些方面| 人成视频在线观看免费观看| 午夜福利免费观看在线| 亚洲久久久国产精品| 欧美黑人精品巨大| 日韩精品免费视频一区二区三区| 日本五十路高清| 亚洲专区中文字幕在线| 欧美乱色亚洲激情| 亚洲男人的天堂狠狠| 中出人妻视频一区二区| 久久香蕉国产精品| 久久久久亚洲av毛片大全| 国产精品久久久久成人av| 看黄色毛片网站| 日日夜夜操网爽| 国产精品 欧美亚洲| 欧美成人免费av一区二区三区| 欧美黄色片欧美黄色片| 日本黄色视频三级网站网址| 无限看片的www在线观看| 国产av一区在线观看免费| 夜夜夜夜夜久久久久| 精品无人区乱码1区二区| 男人操女人黄网站| 欧美激情高清一区二区三区| 80岁老熟妇乱子伦牲交| 久久久久亚洲av毛片大全| 黄色成人免费大全| 亚洲专区字幕在线| 日日夜夜操网爽| 国产精品久久电影中文字幕| 嫁个100分男人电影在线观看| 日本a在线网址| 亚洲精品粉嫩美女一区| 久久国产亚洲av麻豆专区| 国产区一区二久久| 香蕉国产在线看| 亚洲七黄色美女视频| 国产精品一区二区三区四区久久 | 国产精品影院久久| 岛国在线观看网站| netflix在线观看网站| 久久影院123| av网站在线播放免费| 日本一区二区免费在线视频| 真人一进一出gif抽搐免费| 老鸭窝网址在线观看| 亚洲av熟女| 久久人妻福利社区极品人妻图片| 亚洲精品国产精品久久久不卡| 黄色女人牲交| 成人永久免费在线观看视频| netflix在线观看网站| 久久草成人影院| 性少妇av在线| 国产精品av久久久久免费| 亚洲男人天堂网一区| 欧美午夜高清在线| 黄色视频不卡| 女生性感内裤真人,穿戴方法视频| 亚洲男人天堂网一区| 国产精品国产av在线观看| 亚洲少妇的诱惑av| 在线视频色国产色| 亚洲精品在线观看二区| 一进一出抽搐gif免费好疼 | 性色av乱码一区二区三区2| 黄频高清免费视频| 亚洲一码二码三码区别大吗| 国产亚洲精品久久久久5区| 国产av精品麻豆| 亚洲性夜色夜夜综合| 五月开心婷婷网| 亚洲人成77777在线视频| 丝袜美腿诱惑在线| 满18在线观看网站| 国产精华一区二区三区| 精品熟女少妇八av免费久了| 久99久视频精品免费| 国产成人一区二区三区免费视频网站| 一进一出抽搐动态| 久久久国产成人精品二区 | 99re在线观看精品视频| 天天影视国产精品| 亚洲国产精品sss在线观看 | 男人舔女人下体高潮全视频| 国产精品久久久人人做人人爽| 欧美激情 高清一区二区三区| 一区二区日韩欧美中文字幕| 国产精品国产av在线观看| avwww免费| 一a级毛片在线观看| 国产高清激情床上av| 少妇被粗大的猛进出69影院| 精品国产美女av久久久久小说| 精品一品国产午夜福利视频| avwww免费| 少妇被粗大的猛进出69影院| 1024香蕉在线观看| 久久中文看片网| 午夜福利一区二区在线看| 国产真人三级小视频在线观看| 最好的美女福利视频网| 国产男靠女视频免费网站| 精品一区二区三卡| 另类亚洲欧美激情| 亚洲,欧美精品.| 91精品三级在线观看| 国内毛片毛片毛片毛片毛片| 精品久久久久久电影网| 中文字幕人妻熟女乱码| 老司机午夜十八禁免费视频| 久久精品国产99精品国产亚洲性色 | 国产精品偷伦视频观看了| 一边摸一边抽搐一进一小说| 国产成人精品无人区| 亚洲九九香蕉| 国产午夜精品久久久久久| 午夜福利在线免费观看网站| av在线天堂中文字幕 | 国产乱人伦免费视频| 久久精品亚洲精品国产色婷小说| 色老头精品视频在线观看| 亚洲成人免费电影在线观看| www.自偷自拍.com| 亚洲欧美日韩无卡精品| 啦啦啦免费观看视频1| 91精品三级在线观看| 美女高潮到喷水免费观看| av在线天堂中文字幕 | 视频在线观看一区二区三区| 最近最新免费中文字幕在线| a在线观看视频网站| 久久久久国产一级毛片高清牌| 精品久久久久久久久久免费视频 | 变态另类成人亚洲欧美熟女 | 国产精品1区2区在线观看.| 精品电影一区二区在线| 国产亚洲精品第一综合不卡| 黄色女人牲交| 亚洲成人免费av在线播放| 日韩欧美一区二区三区在线观看| 成年人免费黄色播放视频| 叶爱在线成人免费视频播放| 午夜免费激情av| 99久久综合精品五月天人人| 脱女人内裤的视频| 又紧又爽又黄一区二区| 精品久久久久久,| 在线免费观看的www视频| 久久久国产精品麻豆| 中文字幕另类日韩欧美亚洲嫩草| 一级片免费观看大全| 国产极品粉嫩免费观看在线| 首页视频小说图片口味搜索| 级片在线观看| 亚洲中文日韩欧美视频| 国产亚洲av高清不卡| 12—13女人毛片做爰片一| 亚洲第一av免费看| 午夜福利一区二区在线看| 俄罗斯特黄特色一大片| 久久精品91无色码中文字幕| 久久午夜亚洲精品久久| 日韩一卡2卡3卡4卡2021年| 久久精品国产99精品国产亚洲性色 | 色婷婷久久久亚洲欧美| 一边摸一边做爽爽视频免费| 午夜免费成人在线视频| 精品国产美女av久久久久小说| 久久香蕉国产精品| 精品免费久久久久久久清纯| 别揉我奶头~嗯~啊~动态视频| 精品久久久久久成人av| netflix在线观看网站| 欧美日韩精品网址| 久久人人97超碰香蕉20202| videosex国产| 日韩有码中文字幕| 少妇的丰满在线观看| 高潮久久久久久久久久久不卡| 一级毛片高清免费大全| 亚洲专区国产一区二区| 久久久国产欧美日韩av| 亚洲精品粉嫩美女一区| 性色av乱码一区二区三区2| 别揉我奶头~嗯~啊~动态视频| 多毛熟女@视频| 香蕉久久夜色| 身体一侧抽搐| 亚洲视频免费观看视频| av在线天堂中文字幕 | 级片在线观看| 91在线观看av| 一区二区三区激情视频| 日韩 欧美 亚洲 中文字幕| 国产一区二区在线av高清观看| 中文字幕人妻丝袜制服| 亚洲av电影在线进入| 国产精品九九99| 久久天躁狠狠躁夜夜2o2o| 中国美女看黄片| 最近最新免费中文字幕在线| 欧美日韩视频精品一区| 如日韩欧美国产精品一区二区三区| 婷婷丁香在线五月| 美女福利国产在线| 在线观看www视频免费| av在线播放免费不卡| 国产伦人伦偷精品视频| 国产国语露脸激情在线看| 50天的宝宝边吃奶边哭怎么回事| 欧美日韩瑟瑟在线播放| 国产深夜福利视频在线观看| 男女床上黄色一级片免费看| 国产成人精品无人区| 欧美黄色淫秽网站| 窝窝影院91人妻| 丰满饥渴人妻一区二区三| 国产在线观看jvid| 少妇裸体淫交视频免费看高清 | 久久欧美精品欧美久久欧美| 琪琪午夜伦伦电影理论片6080| 在线观看日韩欧美| 亚洲 欧美一区二区三区| 日日爽夜夜爽网站| 亚洲熟女毛片儿| 另类亚洲欧美激情| 午夜福利在线观看吧| 老司机福利观看| 一本综合久久免费| 涩涩av久久男人的天堂| 久久久久久大精品| 99久久国产精品久久久| 久久青草综合色| 亚洲午夜精品一区,二区,三区| 国产成人系列免费观看| 丝袜美腿诱惑在线| 女人精品久久久久毛片| 一本大道久久a久久精品| 午夜影院日韩av| 亚洲性夜色夜夜综合| 99久久精品国产亚洲精品| 精品国产国语对白av| 亚洲少妇的诱惑av| 国产av精品麻豆| 女人精品久久久久毛片| 咕卡用的链子| 亚洲全国av大片| 欧美激情极品国产一区二区三区| 在线观看日韩欧美| 一级毛片高清免费大全| 久久人人97超碰香蕉20202| 久久久国产精品麻豆| 99久久国产精品久久久| av国产精品久久久久影院| 老汉色∧v一级毛片| 亚洲熟妇中文字幕五十中出 | 亚洲欧美一区二区三区久久| 国产xxxxx性猛交| 成人亚洲精品一区在线观看| 成人手机av| 精品久久久久久成人av| 免费在线观看视频国产中文字幕亚洲| av天堂在线播放| 亚洲狠狠婷婷综合久久图片| 欧美日韩一级在线毛片| 精品久久久久久久久久免费视频 | 悠悠久久av| 美女高潮喷水抽搐中文字幕| 在线观看午夜福利视频| 免费久久久久久久精品成人欧美视频| 99国产精品一区二区三区| av天堂在线播放| 淫秽高清视频在线观看| 热99re8久久精品国产| 亚洲国产欧美网| 日韩精品青青久久久久久| 日韩欧美免费精品| 亚洲第一欧美日韩一区二区三区| 老司机亚洲免费影院| 精品人妻在线不人妻| 久久久久国内视频| 搡老乐熟女国产| 亚洲人成网站在线播放欧美日韩| 欧美黄色片欧美黄色片| 99精品久久久久人妻精品| 亚洲 欧美一区二区三区| 亚洲全国av大片| 淫秽高清视频在线观看| 精品熟女少妇八av免费久了| 亚洲情色 制服丝袜| 国产亚洲av高清不卡| 欧美激情 高清一区二区三区| 国产精华一区二区三区| 亚洲情色 制服丝袜| 午夜精品国产一区二区电影| 欧美激情极品国产一区二区三区| 自拍欧美九色日韩亚洲蝌蚪91| 成人特级黄色片久久久久久久| 亚洲欧美激情在线| 亚洲自拍偷在线| 神马国产精品三级电影在线观看 | 久久中文字幕人妻熟女| 啦啦啦 在线观看视频| 黑丝袜美女国产一区| 日本一区二区免费在线视频| 亚洲欧美日韩高清在线视频| 18禁观看日本| 一夜夜www| √禁漫天堂资源中文www| 精品一区二区三区视频在线观看免费 | 精品国产超薄肉色丝袜足j| 亚洲熟妇熟女久久| 国产区一区二久久| 久久精品亚洲熟妇少妇任你| 精品一品国产午夜福利视频| 麻豆av在线久日| 久久人妻av系列| 纯流量卡能插随身wifi吗| 国内毛片毛片毛片毛片毛片| 欧美日韩亚洲国产一区二区在线观看| 亚洲五月色婷婷综合| 欧美激情极品国产一区二区三区| 最新在线观看一区二区三区| 午夜激情av网站| 午夜福利,免费看| 久久久精品国产亚洲av高清涩受| 一本大道久久a久久精品| 欧美午夜高清在线| 国产欧美日韩综合在线一区二区| 国产成+人综合+亚洲专区| 久久天堂一区二区三区四区| 99香蕉大伊视频| 在线十欧美十亚洲十日本专区| 欧美午夜高清在线| 日本wwww免费看| 老熟妇乱子伦视频在线观看| 韩国av一区二区三区四区| 一区二区三区激情视频| 国产精品爽爽va在线观看网站 | 色婷婷久久久亚洲欧美| 美女高潮喷水抽搐中文字幕| 欧美中文日本在线观看视频| 亚洲三区欧美一区| 免费高清在线观看日韩| 法律面前人人平等表现在哪些方面| 黄色 视频免费看| 欧美成人性av电影在线观看| 人人妻人人添人人爽欧美一区卜| 波多野结衣av一区二区av| 一边摸一边抽搐一进一小说| 日韩大尺度精品在线看网址 | 亚洲成人久久性| 高清黄色对白视频在线免费看| 脱女人内裤的视频| 午夜a级毛片| 久久久精品国产亚洲av高清涩受| 老司机在亚洲福利影院| 欧美激情极品国产一区二区三区| 熟女少妇亚洲综合色aaa.| 亚洲成人久久性| 夜夜爽天天搞| 看黄色毛片网站| 两个人免费观看高清视频| 免费在线观看黄色视频的| 一二三四在线观看免费中文在| 视频区欧美日本亚洲| 精品第一国产精品| 美女午夜性视频免费| 国产av又大| 亚洲色图av天堂| 亚洲人成77777在线视频| 丰满迷人的少妇在线观看| 久久久国产欧美日韩av| 又大又爽又粗| 中文字幕精品免费在线观看视频| 国产一区二区三区综合在线观看| 老司机福利观看| 一进一出抽搐gif免费好疼 | 精品久久久精品久久久| 国产精品香港三级国产av潘金莲| 99久久国产精品久久久| 国产乱人伦免费视频| 久久久久久大精品| 久久午夜综合久久蜜桃| 久久国产乱子伦精品免费另类| 欧美色视频一区免费| 欧美不卡视频在线免费观看 | 国产成年人精品一区二区 | 亚洲全国av大片| 久久精品91无色码中文字幕| 美国免费a级毛片| 国产精品一区二区三区四区久久 | 欧美人与性动交α欧美软件| 在线天堂中文资源库| 精品久久久久久成人av| 婷婷精品国产亚洲av在线| 在线观看免费日韩欧美大片| 国产成+人综合+亚洲专区| 99在线人妻在线中文字幕| 男女高潮啪啪啪动态图| 亚洲欧美日韩无卡精品| 欧美日韩福利视频一区二区| 成年版毛片免费区| 中文字幕另类日韩欧美亚洲嫩草| 人人澡人人妻人| 亚洲精品久久成人aⅴ小说| 国产成人一区二区三区免费视频网站| 女同久久另类99精品国产91| 中文字幕人妻丝袜制服| 丁香六月欧美| 欧美老熟妇乱子伦牲交| 成人三级做爰电影| 欧美+亚洲+日韩+国产| 成人手机av| 欧美成人午夜精品| 1024视频免费在线观看| 亚洲国产精品合色在线| 国产免费现黄频在线看| 黄片播放在线免费| 欧美黑人精品巨大| 国产精品亚洲av一区麻豆| 好男人电影高清在线观看| 高清毛片免费观看视频网站 | 亚洲少妇的诱惑av| 国产色视频综合| 又紧又爽又黄一区二区| 一夜夜www| 美国免费a级毛片| 久久精品人人爽人人爽视色| 亚洲在线自拍视频| 午夜福利一区二区在线看| 久久草成人影院| 免费看a级黄色片| 热99re8久久精品国产| 麻豆久久精品国产亚洲av | 欧美性长视频在线观看| 午夜福利免费观看在线| 免费观看精品视频网站| 日韩欧美三级三区| svipshipincom国产片| 国产真人三级小视频在线观看| 性少妇av在线| 久久人人97超碰香蕉20202| 后天国语完整版免费观看| 亚洲精品成人av观看孕妇| 免费搜索国产男女视频| 欧美日韩瑟瑟在线播放| 亚洲欧美一区二区三区黑人| 搡老熟女国产l中国老女人| 黄色视频,在线免费观看| 欧美日韩国产mv在线观看视频| 在线天堂中文资源库| av国产精品久久久久影院| 老汉色av国产亚洲站长工具| 黄色怎么调成土黄色| 亚洲久久久国产精品| 黄色a级毛片大全视频| 日本精品一区二区三区蜜桃| 久久热在线av| 黄片小视频在线播放| 久久午夜综合久久蜜桃| 亚洲精品一二三| 久久人妻av系列| 99热只有精品国产| 亚洲情色 制服丝袜| 深夜精品福利| 人人妻人人添人人爽欧美一区卜| 男女床上黄色一级片免费看| 淫妇啪啪啪对白视频| 国产精品一区二区三区四区久久 | 一边摸一边做爽爽视频免费| 亚洲黑人精品在线| 国产精品免费视频内射| 久久天躁狠狠躁夜夜2o2o| 一级黄色大片毛片| 桃红色精品国产亚洲av| 在线观看www视频免费| 亚洲色图 男人天堂 中文字幕| 欧美亚洲日本最大视频资源| 欧美人与性动交α欧美精品济南到| 欧美日韩亚洲高清精品| 国产国语露脸激情在线看| 乱人伦中国视频| 亚洲免费av在线视频| av中文乱码字幕在线| 女人被狂操c到高潮| 国产精品久久久av美女十八| 欧美成人免费av一区二区三区| 在线观看免费视频网站a站| 男女高潮啪啪啪动态图| 黄片小视频在线播放| 9色porny在线观看| 中文字幕高清在线视频| 免费高清在线观看日韩| 色婷婷久久久亚洲欧美| 成人手机av| 亚洲自拍偷在线| 日日爽夜夜爽网站| a在线观看视频网站| 美女午夜性视频免费| av国产精品久久久久影院| 夜夜爽天天搞| 精品久久久久久,|