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

    基于Curvelet變換的隧道裂隙水GPR數(shù)據(jù)處理研究

    2014-05-25 00:30:29朱自強(qiáng)魯光銀王凡
    物探化探計(jì)算技術(shù) 2014年5期
    關(guān)鍵詞:波場探地雷達(dá)

    朱自強(qiáng),朱 賀*,魯光銀,王凡,譚 潔

    (1.中南大學(xué)地球科學(xué)與信息物理學(xué)院,長沙 410083;2.中國水電顧問集團(tuán)貴陽勘測設(shè)計(jì)研究院,貴陽 550081)

    基于Curvelet變換的隧道裂隙水GPR數(shù)據(jù)處理研究

    朱自強(qiáng)1,朱 賀1*,魯光銀1,王凡2,譚 潔1

    (1.中南大學(xué)地球科學(xué)與信息物理學(xué)院,長沙 410083;2.中國水電顧問集團(tuán)貴陽勘測設(shè)計(jì)研究院,貴陽 550081)

    應(yīng)用地質(zhì)雷達(dá)探測隧道內(nèi)部裂隙水的分布狀況時(shí),受隧道鋼筋網(wǎng)的影響以及噪聲干擾,導(dǎo)致探測信號(hào)中弱有效信號(hào)難以識(shí)別。Curvelet變換可以對(duì)二維信號(hào)從頻率、角度和空間位置實(shí)現(xiàn)有效反射波和干擾波的分離及降噪處理,應(yīng)用Curvelet變換對(duì)含鋼筋網(wǎng)干擾的裂隙水模型的探地雷達(dá)模擬數(shù)據(jù)進(jìn)行波場分離和降噪處理,在壓制直達(dá)波、去除鋼筋層反射信號(hào)及多次反射的基礎(chǔ)上,可以有效地提取出含裂隙水所引起的反射波信號(hào)。將該方法應(yīng)用到汝郴高速某隧道實(shí)測探地雷達(dá)數(shù)據(jù)中,較好地去除了鋼筋層反射信號(hào)的影響,準(zhǔn)確地判斷出了裂隙水的空間分布位置,證明了所提出方法的有效性。

    Curvelet變換;探地雷達(dá);隧道裂隙水;波場分離;降噪

    0 引言

    探地雷達(dá)法[1-2]是近些年來迅速發(fā)展起來的一種高分辨率、高效率的勘探方法,目前在隧道勘測中,特別是探測巖溶裂隙水方面應(yīng)用較多[3-6]。由于隧道周圍裂隙水被密集的鋼筋層覆蓋,用探地雷達(dá)對(duì)隧道裂隙水探測時(shí)會(huì)受到鋼筋層反射的干擾,因此,需要對(duì)雷達(dá)數(shù)據(jù)進(jìn)行處理,去除干擾和噪聲,從而提取有效的異常信號(hào)。劉斌等[7]將復(fù)信號(hào)技術(shù)應(yīng)用到探地雷達(dá)對(duì)地下裂隙水的探測中,對(duì)雷達(dá)數(shù)據(jù)中的異常分別從振幅、頻率和相位等各方面進(jìn)行綜合分析,提高了預(yù)報(bào)精度和準(zhǔn)確性,但沒能從根本上克服探地雷達(dá)法抗干擾能力弱的缺點(diǎn)。劉四新等[8]用鉆孔探地雷達(dá)對(duì)巖溶裂隙進(jìn)行勘探,勘探效果較好,但具有探測上方位不確定性和成本較高的缺點(diǎn)。鄒海林等[9]基于小波變換對(duì)探地雷達(dá)信號(hào)進(jìn)行處理,李才明[10]基于小波能譜分析對(duì)巖溶區(qū)探地雷達(dá)進(jìn)行目標(biāo)識(shí)別。小波變換很好地解決了時(shí)頻同步問題,而且對(duì)一維分段平穩(wěn)信號(hào)的表達(dá)具有最佳性能,但是對(duì)于描述二維圖像信息卻存在諸多缺點(diǎn)[11-13]。

    Candes等人[14-15]在1999年提出了第一代Curvelet變換,并于2002年提出了實(shí)現(xiàn)更簡單、運(yùn)行速度更快、更便于直觀理解的第二代Curvelet變換算法。在此之后,他們又于2005年提出了兩種基于第二代Curvelet 變換理論的快速離散實(shí)現(xiàn)方法[16],即USFFT(unequally-spaced fast Fourier transforms)和Wrapping兩種快速離散算法。2008年以來,該方法在地球物理資料處理尤其是地震信號(hào)處理中得到了迅速地研究和應(yīng)用[17-20]。在國內(nèi),仝中飛等人[21]提出了迭代閾值法壓制地震信號(hào)中的隨機(jī)噪聲,取得了很好的效果。潘雪輝[22]提出了一種用Curvelet變換壓制地震資料中線性噪聲的方法,為地震信號(hào)中線性噪聲的壓制問題提供了較好地解決方法。

    Curvelet變換建立在小波變換的基礎(chǔ)上,增加了一個(gè)方位參數(shù),解決了小波變化在處理二維信號(hào)時(shí)的不足,而探地雷達(dá)檢測數(shù)據(jù)主要為描述一定深度范圍內(nèi)剖面的二維信號(hào),這讓作者萌生思路,利用Curvelet變換對(duì)探地雷達(dá)數(shù)據(jù)進(jìn)行處理研究。在對(duì)Curvelet變換方法研究基礎(chǔ)上,將其應(yīng)用到隧道裂隙水的探地雷達(dá)檢測資料處理中。本研究首先模擬了在含有覆蓋鋼筋層和裂隙水的混凝土中的探地雷達(dá)波場,并通過Curvelet變換對(duì)模擬波場進(jìn)行去除直達(dá)波、提取水層信息、降噪等一系列數(shù)據(jù)處理流程,提取出有效的裂隙水異常信息;隨后將該方法應(yīng)用到汝郴高速某隧道實(shí)測地質(zhì)雷達(dá)數(shù)據(jù),獲取了隧道內(nèi)裂隙水異常的空間分布,通過鉆孔注漿驗(yàn)證,證明了此方法的有效性。

    1 Curvelet變換基本原理

    小波變換在信號(hào)處理中的應(yīng)用得到了很大的發(fā)展,遺憾的是,由一維小波所生成的可分離小波只具有有限的方向,因此在處理圖像邊緣時(shí)的效果不是很好。Curvelet變換除了與小波變換一樣具有尺度和位移參數(shù),還增加了一個(gè)方位參數(shù),因此具有更好的方位識(shí)別能力,可以精確表達(dá)出圖像中邊緣的方向信息。

    1.1 連續(xù)Curvelet變換

    Curvelet變換是利用基函數(shù)與信號(hào)(或函數(shù))的內(nèi)積形式來實(shí)現(xiàn)信號(hào)(或函數(shù))的稀疏表示,表示為:

    其中:φj,k,l表示Curvelet函數(shù);j、k、l分別表示尺度、位置和方向參數(shù)。

    Curvelet變換在頻域內(nèi)采用窗函數(shù)U來實(shí)現(xiàn)的。分別定義徑向窗函數(shù)W(r)、r∈(1/2,2)和角度窗函數(shù)V(t),t∈[-1,1],二者滿足式(2)與式(3)。

    對(duì)于每一個(gè)j≥j0,在頻域中定義窗函數(shù)Uj如式(4)。

    綜上所述,我們可以得到Curvelet變換公式:

    圖1 Curvelet變換示意圖Fig.1 Curvelet transform schematic diagram

    1.2 離散Curvelet變換

    在直角坐標(biāo)系下f[t1,t2],0≤t1,t2<n為輸入,Curvelet變換離散形式為:

    1.3 離散Curvelet變換快速實(shí)現(xiàn)方法

    本研究采用基于USFFT算法的快速離散Cur-velet變換方法,實(shí)現(xiàn)過程如下:

    1)對(duì)于給定的直角坐標(biāo)系下二維函數(shù)f[t1,t2],0≤t1,t2<ω進(jìn)行2DFFT,得到其二維頻率域表達(dá)式:

    2)在頻率域,對(duì)每一對(duì)(i,j),對(duì)函數(shù)^f[n1,n2]重新采樣,得到函數(shù)采樣值:

    其中:Pj={(n1,n2):n1,0≤n1≤n1,0+L1,j,n2,0≤n2≤n2,0+L2,j};L1,j表示窗函數(shù)支撐區(qū)間的長度;L2,j表示支撐區(qū)間的寬度。

    1.4 Curvelet變換的非線性逼近能力

    Curvelet變換之所以比小波變換更適合表示二維信號(hào),其主要原因在于Curvelet變換比小波在二維平面具有更好的稀疏性。假設(shè)f∈L2[0,1]為一條曲線的光滑部分,分別用小波變換和Curvelet變換對(duì)其分解,再分別選擇最佳的M個(gè)系數(shù)對(duì)函數(shù)進(jìn)行重構(gòu),得到三種變換下的的最佳逼近它們的逼近誤差分別為:

    從重構(gòu)誤差式中可以看出,Curvelet變換的逼近率要優(yōu)于小波變換。Curvelet變換之所在表示二維曲線時(shí)有更好的逼近率,主要在于它在小波變換的基礎(chǔ)上增加了方位參數(shù),使得它具有多方向性、多尺度性、波動(dòng)性和各向異性的特點(diǎn),這也使得用更少的Curvelet系數(shù)就能表示二維信號(hào)的主要特征。

    圖2 小波和Curvelet對(duì)二維曲線逼近示意圖Fig.2 Schematic of two-dimensional curve approximation with wavelet and Curvelet

    2 模擬數(shù)據(jù)處理分析

    針對(duì)表層存在鋼筋干擾下的含裂隙水混凝土模型進(jìn)行數(shù)值模擬,模型大小為5.2 m×1.28 m,混凝土介電常數(shù)ε0為9,電導(dǎo)率σ0為0.01 s/m。距離表層0.3 m處設(shè)置一排鋼筋,半徑0.01 m。鋼筋層下方設(shè)置一個(gè)水平和一個(gè)傾斜含裂隙水,橫截面大小均為0.8 m×0.08 m,介電常數(shù)ε1為81,電導(dǎo)率σ1為0.05 s/m,模擬天線頻率為400 MHz。模型剖面和雷達(dá)模擬剖面如圖3所示。

    圖3 雷達(dá)數(shù)值模擬Fig.3 GPR numerical simulation

    在圖3中,表層直達(dá)波強(qiáng)度較高,弱化了下面鋼筋層和裂隙水層的反射信號(hào)。將總波場變換到Curvelet域后,然后在角度窗函數(shù)下選擇變換域中θ=0(θ為方向因子)附近的Curvelet系數(shù)進(jìn)行重構(gòu),得到直達(dá)波的波場信號(hào),進(jìn)而分離得到的直達(dá)波波場信號(hào)和剩余波場信號(hào)(圖4)。

    圖4 波場分離Fig.4 Wave field separation

    當(dāng)直達(dá)波去除后,下面鋼筋層和裂隙水層的反射信息加強(qiáng),但是二次反射和干擾信息的存在干擾了異常信號(hào),通過Curvelet變換的降噪功能可以去掉二次反射及其他干擾信號(hào)。基于Curvelet變換的特點(diǎn),前人提出了各種關(guān)于降噪的閾值方法[23-24],本研究采取一種自適應(yīng)閾值方法[25],即隨著尺度的改變,每一層閾值可以相應(yīng)地改變。閾值函數(shù)表達(dá)式如下:其中 sigma為噪聲的方差估計(jì)值;N為Curvelet域系數(shù)矩陣大?。籹為Curvelet域尺度分解下的層數(shù)。

    圖5 波場降噪Fig.5 Wavefield denoising

    圖5為sigma系數(shù)設(shè)為700時(shí)降噪的效果圖,降噪后,鋼筋和裂隙水的異常信息已非常明顯。原始二維信號(hào)在經(jīng)過Curvelet變換后會(huì)得到的一系列Curvelet系數(shù)矩陣,將這些系數(shù)矩陣在MATLAB中生成時(shí)間-頻率圖像,通過對(duì)比原始圖像異常位置可以識(shí)別出包含有效信號(hào)和干擾信號(hào)的矩陣。如果我們只需要裂隙水的異常信息,只需對(duì)這些離散化的系數(shù)矩陣進(jìn)行人工識(shí)別,并挑選出含裂隙水信息的系數(shù)矩陣,進(jìn)行Curvelet反變換,得到只含裂隙水異常信息的二維數(shù)據(jù)。圖6為應(yīng)用Curvelet變換波場分離方法得到的裂隙水異常信息數(shù)據(jù)。

    圖6 提取裂隙水波場信息Fig.6 Wave signal extraction of fissure water

    處理結(jié)果顯示,用Curvelet變換對(duì)雷達(dá)模擬數(shù)據(jù)的處理效果良好,波場分離和降噪都能夠基本達(dá)到目的。兩個(gè)方向有差異的裂隙水的異常信息在處理過程中表現(xiàn)出不同的處理效果,這也說明了Curvelet變換在處理圖像信息時(shí)帶有很強(qiáng)的方向選擇性。由于水平的裂隙水在方向上與上面鋼筋層保持平行,用Curvelet變換離散化圖像信息時(shí),水平裂隙水與鋼筋層的信息在角度窗函數(shù)選擇下沒能夠完全分離,只能夠通過尺度來分離部分信息。這樣傾斜的裂隙水異常信息才會(huì)顯得比水平裂隙更明顯。

    3 實(shí)測數(shù)據(jù)處理分析

    在Curvelet變換方法研究和模擬處理的基礎(chǔ)上,對(duì)汝郴高速公路某隧道內(nèi)采集所得探地雷達(dá)數(shù)據(jù)進(jìn)行處理。探測目的是找出隧道邊墻和底部存在的裂隙水,確定這些裂隙水的位置及相關(guān)信息,為工程治理工作提供較好的參考信息。探測中使用的天線為400 MHz,采樣點(diǎn)數(shù)設(shè)為512,參考設(shè)計(jì)方案,鋼筋間距約為50 cm。

    圖7(a)為截取采集數(shù)據(jù)中某一異常位置信息,直達(dá)波和鋼筋產(chǎn)生的干擾信息強(qiáng)度較強(qiáng),同時(shí)數(shù)據(jù)中含有一些噪聲干擾。首先我們要將有效異常信息和鋼筋層干擾信息及直達(dá)波信息分離開,提取含有有效異常信息的數(shù)據(jù)。圖7(b)為分離出的鋼筋層信息及一些干擾信息,圖7(c)為提取得到的含有裂隙水層異常信息的波場。由于波場分離更多在于角度窗下的選擇,圖7(b)中鋼筋層干擾和直達(dá)波與圖7(c)中有效信息被基本分離開,但仍然有一些因?yàn)榕c裂隙水異常信息方向一致而未被分離的遺留信息,且存在一些噪聲。遺留下的干擾信息強(qiáng)度較弱,可以通過Curvelet降噪功能,應(yīng)用自適應(yīng)尺度閾值法,通過尺度窗的選擇進(jìn)行降噪,取sigma系數(shù)為48 000得到降噪后的異常圖7(d)。圖7(d)異常信息清晰,進(jìn)而可以判斷裂隙水的空間位置。

    根據(jù)探測和處理結(jié)果,對(duì)發(fā)現(xiàn)的異常位置進(jìn)行鉆孔驗(yàn)證如圖8所示。驗(yàn)證結(jié)果支持了該方法對(duì)于裂隙水空間位置的判斷,并在該異常位置進(jìn)行了注漿治理,總注漿量為10 T左右水泥和38 T左右化學(xué)漿。對(duì)于本次探測中發(fā)現(xiàn)的其他幾個(gè)異常位置,用文中所述方法進(jìn)行分析并打鉆注漿,在判斷得到的異常較大位置都能注入不同量的漿液,驗(yàn)證了該處理方法的準(zhǔn)確性與可行性。

    4 結(jié)論

    1)基于Curvelet變換,從波場分離的角度來提高探地雷達(dá)異常信號(hào)的辨識(shí)度,較好地解決了直達(dá)波和鋼筋層強(qiáng)反射對(duì)下層異常信息的干擾問題。

    2)角度參數(shù)的增加使得Curvelet變換在二維信號(hào)的分離上有了很強(qiáng)的優(yōu)勢,但是在干擾信息與有效信息方向性一致的情況下,分離效果并不理想。

    圖7 工程應(yīng)用實(shí)例Fig.7 Examples of application engineering

    圖8 現(xiàn)場鉆孔照片F(xiàn)ig.8 Photos of drilling

    3)閾值函數(shù)的選取直接關(guān)系到Curvelet降噪的效果,采取適應(yīng)不同尺度的自適應(yīng)閾值函數(shù),基本能夠達(dá)到消除噪聲的目的,但閾值函數(shù)thresh沒有從角度去適應(yīng)系數(shù)矩陣,找到同時(shí)適應(yīng)尺度和角度的閾值函數(shù)將是下一步急需解決的問題。

    [1] 李大心.探地雷達(dá)方法與應(yīng)用[M].北京:地質(zhì)出版社,1994.

    [2] 楊天春,呂紹林,王齊仁.探地雷達(dá)檢測道路厚度結(jié)構(gòu)的應(yīng)用現(xiàn)狀及進(jìn)展[J].物探與化探,2003,27(1):79-82.

    [3] 楊龍江,閻滿存.探地雷達(dá)在隱伏巖溶區(qū)工程地質(zhì)勘探中的應(yīng)用[J].物探化探計(jì)算技術(shù),1999,21(2):172-176.

    [4] 葛雙成,邵長云.巖溶勘察中的探地雷達(dá)技術(shù)及應(yīng)用[J].地球物理學(xué)進(jìn)展,2005,20(2):476-481.

    [5] 李孟娟,李川.探地雷達(dá)檢測隧道襯砌厚度的研究[J].物探化探計(jì)算技術(shù),2008,30(3):231-234.

    [6] 李仁海,楊磊,許新剛,等.地質(zhì)雷達(dá)探測技術(shù)在巖溶地形勘察中的應(yīng)用[J].物探化探計(jì)算技術(shù),2009,31(5):442-446.

    [7] 劉斌,李術(shù)才,李樹忱,等.復(fù)信號(hào)分析技術(shù)在地質(zhì)雷達(dá)預(yù)報(bào)巖溶裂隙水中的應(yīng)用研究[J].巖土力學(xué),2009,30(7):2191-2196.

    [8] 劉四新,曾昭發(fā),徐波.利用鉆孔雷達(dá)探測地下含水裂縫[J].地球物理學(xué)進(jìn)展,2006,21(2):620-624.

    [9] 鄒海林,寧書年,林捷.小波理論在探地雷達(dá)信號(hào)處理中的應(yīng)用[J].地球物理學(xué)進(jìn)展,2004,19(2):268-275.

    [10] 李才明,王良書,徐鳴潔,等.基于小波能譜分析的巖溶區(qū)探地雷達(dá)目標(biāo)識(shí)別[J].地球物理學(xué)報(bào),2006,49(5):1499-1504.

    [11]WALKER J S,CHEN Y.Image denoising using treebased wavelet subband correlations and shrinkage[J].Optical Engineering,2000,39(11):2900-2908.

    [12]楊鳳娟,羅省賢.小波變換在探地雷達(dá)檢測鋼筋中的應(yīng)用[J].物探化探計(jì)算技術(shù),2009,31(4):354-360.

    [13]黃敏,朱德兵,郭政學(xué),等.連續(xù)小波變換在探地雷達(dá)信號(hào)分析中的應(yīng)用研究[J].物探化探計(jì)算技術(shù),2012,34(5):593-598.

    [14]CANDES E J,DONOHO D L.Curvelets:A surprisingly effective nonadaptive representation for objects with edges[R].DTIC Document,2000.

    [15]STARCK J,CANDèS E J,DONOHO D L.The curvelet transform for image denoising[J].Image Processing,IEEE Transactions on.2002,11(6):670-684.

    [16]CANDES E,DEMANET L,DONOHO D,et al.Fast discrete curvelet transforms[J].Multiscale Modeling&Simulation,2006,5(3):861-899.

    [17]Neelamani R,Baumstein A I,Gillard D G,et al.Coherent and random noise attenuation using the curvelet transform[J].The Leading Edge.2008,27(2):240-248.

    [18]姜宇東,楊勤勇,何柯,等.基于曲波變換的地面微地震資料去噪方法研究[J].石油物探,2012,51(6):620-624.

    [19]吳愛弟,趙秀玲.基于曲波變換的地震信號(hào)去噪方法[J].中國石油大學(xué)學(xué)報(bào):自然科學(xué)版,2010,34(003):30-33.

    [20]張?jiān)茝?qiáng),張培林,王國德,等.基于曲波變換和色度模型的彩色圖像去噪[J].中國圖象圖形學(xué)報(bào),2012,17(012):1472-1477.

    [21]仝中飛,王德利,劉冰.基于Curvelet變換閾值法的地震數(shù)據(jù)去噪方法[J].吉林大學(xué)學(xué)報(bào):地球科學(xué)版,2008:48-52.

    [22]潘雪輝.曲波變換在地震信號(hào)去噪中的應(yīng)用[D].西安:西安理工大學(xué),2009.

    [23]HERRMANN F J,HENNENFENT G.Non-parametric seismic data recovery with curvelet frames[J].Geophysical Journal International,2008,173(1):233-248.

    [24]HERRMANN F J,MOGHADDAM P,STOLK C C.Sparsity-and continuity-promoting seismic image recovery with curvelet frames[J].Applied and Computational Harmonic Analysis,2008,24(2):150-173.

    [25]劉靜.地震圖像噪聲壓制在Curvelet域的閾值選取與噪聲定位[D].成都:西南交通大學(xué),2011.

    Processing of GPR data in tunnel fissure water based on Curvelet transform

    ZHU Zi-qiang1,ZHU He1*,LU Guang-yin1,WANG Fan2,TAN Jie1
    (1.School of Geosciences and Info-Physics,Central South University,Changsha 410083,China;2.HYDROChina Guiyang Engineering Corporation,Guiyang 550081,China)

    When we detect the distribution of fissure water in tunnel with ground penetrating radar,it is hard to identify the weak effective signal because of the influence of mesh reinforcement and noise.Curvelet transform can achieve effective wave field separation and denoising of two-dimensional signal from the frequency,angle and spatial location.The Curvelet transform is applied to wave field separation and signal denoise processing by using of GPR simulated data from fissured water which is disturbed seriously by the existence of mesh reinforcement.After the direct signal suppression,removal of the reflected signal and multiple reflections from mesh reinforcement,we can effectively extract the reflected wave signal which is caused by the water bearing structure crack.Our method is applied to real data measured in a tunnel of the Ruchen highway.The results show that after removal of the reflected signal from mesh reinforcement we can accurately determine the location of fissure water,which proves the effectiveness of the proposed method.

    Curvelet transform;ground-penetrating radar(GPR);fissure water in tunnel;wavefield separation;denoising

    P 631.4

    A

    10.3969/j.issn.1001-1749.2014.05.10

    1001-1749(2014)05-0571-06

    2014-03-33 改回日期:2014-06-22

    國家自然科學(xué)基金項(xiàng)目(41174061);中南大學(xué)自由探索計(jì)劃(2011QNZT011)

    朱自強(qiáng)(1964-),男,教授,博士生導(dǎo)師,主要從事地質(zhì)災(zāi)害探測與監(jiān)測工作,E-mail:13507319431@139.com。

    *通訊作者:朱賀(1988-),男,碩士,主要從事地球物理信號(hào)處理研究,E-mail:zhuhe24@sina.com。

    猜你喜歡
    波場探地雷達(dá)
    有雷達(dá)
    大自然探索(2023年7期)2023-08-15 00:48:21
    基于超表面的探地雷達(dá)增強(qiáng)探測研究
    全極化探地雷達(dá)系統(tǒng)
    雷達(dá)
    彈性波波場分離方法對(duì)比及其在逆時(shí)偏移成像中的應(yīng)用
    基于探地雷達(dá)法的地下管線探測頻譜分析
    交錯(cuò)網(wǎng)格與旋轉(zhuǎn)交錯(cuò)網(wǎng)格對(duì)VTI介質(zhì)波場分離的影響分析
    基于Hilbert變換的全波場分離逆時(shí)偏移成像
    基于空時(shí)二維隨機(jī)輻射場的彈載雷達(dá)前視成像
    現(xiàn)代“千里眼”——雷達(dá)
    欧美国产精品一级二级三级| 一区二区三区精品91| 亚洲婷婷狠狠爱综合网| 18禁裸乳无遮挡动漫免费视频| 亚洲天堂av无毛| 国产高清国产精品国产三级| 一区二区三区激情视频| 天美传媒精品一区二区| 极品人妻少妇av视频| 亚洲一码二码三码区别大吗| 大陆偷拍与自拍| 午夜福利,免费看| 日韩电影二区| 搡老岳熟女国产| 日日啪夜夜爽| 国产爽快片一区二区三区| 青春草国产在线视频| 中文字幕精品免费在线观看视频| 亚洲精品av麻豆狂野| 成人漫画全彩无遮挡| 少妇人妻久久综合中文| 哪个播放器可以免费观看大片| 母亲3免费完整高清在线观看| 欧美日韩综合久久久久久| 99re6热这里在线精品视频| 国产又色又爽无遮挡免| 久久久精品免费免费高清| 成年美女黄网站色视频大全免费| 日韩熟女老妇一区二区性免费视频| a级片在线免费高清观看视频| 免费黄频网站在线观看国产| 波多野结衣av一区二区av| 亚洲精品在线美女| 老司机在亚洲福利影院| 超碰成人久久| 中文字幕精品免费在线观看视频| 一级毛片我不卡| 咕卡用的链子| 咕卡用的链子| 大片免费播放器 马上看| 高清不卡的av网站| 国产99久久九九免费精品| 日本av手机在线免费观看| 男女高潮啪啪啪动态图| 免费久久久久久久精品成人欧美视频| e午夜精品久久久久久久| 别揉我奶头~嗯~啊~动态视频 | 久久久久精品国产欧美久久久 | av女优亚洲男人天堂| 午夜福利影视在线免费观看| www日本在线高清视频| 国产av国产精品国产| 午夜久久久在线观看| 亚洲国产中文字幕在线视频| 国产精品亚洲av一区麻豆 | 99久久人妻综合| 国产在线一区二区三区精| 国产亚洲精品第一综合不卡| 男人添女人高潮全过程视频| 国产免费又黄又爽又色| 1024香蕉在线观看| 国产精品香港三级国产av潘金莲 | www.自偷自拍.com| 国产爽快片一区二区三区| 91成人精品电影| 老熟女久久久| tube8黄色片| 国产淫语在线视频| 欧美黑人精品巨大| 人人妻人人爽人人添夜夜欢视频| 久久天堂一区二区三区四区| 一区二区日韩欧美中文字幕| 人妻 亚洲 视频| 欧美老熟妇乱子伦牲交| 亚洲精品久久成人aⅴ小说| 午夜福利乱码中文字幕| 大香蕉久久网| 国产有黄有色有爽视频| 制服丝袜香蕉在线| 老鸭窝网址在线观看| 亚洲一区二区三区欧美精品| 亚洲精品自拍成人| 制服诱惑二区| 午夜老司机福利片| 亚洲欧洲国产日韩| 成人国产麻豆网| av又黄又爽大尺度在线免费看| 中文精品一卡2卡3卡4更新| 久久久精品94久久精品| 制服人妻中文乱码| 久久热在线av| 精品少妇内射三级| 亚洲欧美清纯卡通| 亚洲四区av| 最近的中文字幕免费完整| 亚洲精品中文字幕在线视频| 亚洲国产成人一精品久久久| 精品人妻一区二区三区麻豆| avwww免费| 国产伦理片在线播放av一区| 人体艺术视频欧美日本| 大香蕉久久成人网| 亚洲美女视频黄频| 久久久精品94久久精品| 精品久久久精品久久久| 免费黄色在线免费观看| 最近的中文字幕免费完整| videos熟女内射| 国产又色又爽无遮挡免| 国产有黄有色有爽视频| 久久久久精品久久久久真实原创| 多毛熟女@视频| 啦啦啦在线免费观看视频4| 人妻人人澡人人爽人人| 又黄又粗又硬又大视频| 在线观看www视频免费| 国产淫语在线视频| 精品国产露脸久久av麻豆| 黑丝袜美女国产一区| 19禁男女啪啪无遮挡网站| 天天操日日干夜夜撸| 国产国语露脸激情在线看| 精品一区在线观看国产| av国产精品久久久久影院| 夜夜骑夜夜射夜夜干| 飞空精品影院首页| 中文字幕高清在线视频| 王馨瑶露胸无遮挡在线观看| 不卡av一区二区三区| 亚洲国产成人一精品久久久| 一区福利在线观看| 日韩熟女老妇一区二区性免费视频| 最近最新中文字幕大全免费视频 | 久久精品亚洲av国产电影网| 成人国产av品久久久| 久久久久久久久免费视频了| 王馨瑶露胸无遮挡在线观看| 色婷婷久久久亚洲欧美| 亚洲熟女毛片儿| 国产深夜福利视频在线观看| 日韩欧美一区视频在线观看| 国产又色又爽无遮挡免| 丝袜美腿诱惑在线| 亚洲精品aⅴ在线观看| 国产 一区精品| 国产精品女同一区二区软件| 咕卡用的链子| 飞空精品影院首页| 少妇被粗大的猛进出69影院| 日韩,欧美,国产一区二区三区| 中文字幕人妻丝袜一区二区 | 伦理电影免费视频| 好男人视频免费观看在线| 亚洲欧美一区二区三区国产| 欧美日韩精品网址| netflix在线观看网站| 亚洲国产欧美网| 另类精品久久| 亚洲第一av免费看| www.av在线官网国产| 国产精品 国内视频| 在线观看免费午夜福利视频| 午夜激情久久久久久久| 精品久久久精品久久久| 国产乱来视频区| 国产野战对白在线观看| 九草在线视频观看| 大话2 男鬼变身卡| 街头女战士在线观看网站| 国产国语露脸激情在线看| 欧美国产精品一级二级三级| 欧美国产精品va在线观看不卡| 国产精品久久久久久精品电影小说| 蜜桃国产av成人99| 久久久久久久国产电影| 午夜福利视频在线观看免费| 国产精品香港三级国产av潘金莲 | 精品少妇黑人巨大在线播放| 国产精品.久久久| 亚洲av国产av综合av卡| 操出白浆在线播放| 国产成人欧美在线观看 | 国产精品一国产av| 国产男女内射视频| 99久久综合免费| 欧美亚洲日本最大视频资源| 在线观看一区二区三区激情| 欧美老熟妇乱子伦牲交| 丝袜美腿诱惑在线| 亚洲一卡2卡3卡4卡5卡精品中文| 国产成人免费观看mmmm| 久久久精品国产亚洲av高清涩受| 黑人猛操日本美女一级片| 国产1区2区3区精品| 亚洲精品日本国产第一区| 蜜桃国产av成人99| 亚洲人成网站在线观看播放| 大话2 男鬼变身卡| 亚洲一级一片aⅴ在线观看| 毛片一级片免费看久久久久| 国产精品国产三级国产专区5o| 久热爱精品视频在线9| 黄色 视频免费看| 男女午夜视频在线观看| 免费看不卡的av| 久久国产精品大桥未久av| av在线观看视频网站免费| 街头女战士在线观看网站| 久久ye,这里只有精品| 热99国产精品久久久久久7| 激情视频va一区二区三区| 晚上一个人看的免费电影| 尾随美女入室| 亚洲一码二码三码区别大吗| 亚洲色图 男人天堂 中文字幕| 精品久久蜜臀av无| 国产免费一区二区三区四区乱码| 欧美日韩精品网址| 国产成人系列免费观看| 美女脱内裤让男人舔精品视频| 国产熟女午夜一区二区三区| svipshipincom国产片| 王馨瑶露胸无遮挡在线观看| 秋霞在线观看毛片| a级片在线免费高清观看视频| 最新的欧美精品一区二区| 精品久久久久久电影网| 久久青草综合色| 超碰成人久久| 久久人人97超碰香蕉20202| 尾随美女入室| 亚洲激情五月婷婷啪啪| 亚洲免费av在线视频| 国产黄频视频在线观看| 最黄视频免费看| 人人澡人人妻人| 国产黄色视频一区二区在线观看| 国产欧美日韩一区二区三区在线| 国产又色又爽无遮挡免| 天天操日日干夜夜撸| 又大又黄又爽视频免费| 中国国产av一级| 18禁动态无遮挡网站| tube8黄色片| 黄频高清免费视频| 视频区图区小说| 看免费av毛片| 女人被躁到高潮嗷嗷叫费观| 日韩视频在线欧美| 母亲3免费完整高清在线观看| 精品久久蜜臀av无| 超碰成人久久| 国产精品麻豆人妻色哟哟久久| 国产在线免费精品| 国产极品天堂在线| 又粗又硬又长又爽又黄的视频| 美国免费a级毛片| 久久久久久久国产电影| 亚洲久久久国产精品| 亚洲第一区二区三区不卡| 日本vs欧美在线观看视频| 欧美乱码精品一区二区三区| 亚洲av成人不卡在线观看播放网 | av在线播放精品| 中文字幕亚洲精品专区| 丝袜美足系列| 黄色视频在线播放观看不卡| 在线观看人妻少妇| 久久精品国产亚洲av涩爱| 亚洲精品日韩在线中文字幕| 一级毛片电影观看| 一级毛片我不卡| 亚洲欧美精品综合一区二区三区| 午夜福利视频在线观看免费| 成人午夜精彩视频在线观看| 男女免费视频国产| 精品久久蜜臀av无| 一边摸一边抽搐一进一出视频| 美女中出高潮动态图| 在线看a的网站| 少妇的丰满在线观看| 亚洲综合精品二区| 在线观看人妻少妇| 天堂中文最新版在线下载| 多毛熟女@视频| 麻豆精品久久久久久蜜桃| 一本一本久久a久久精品综合妖精| 十八禁人妻一区二区| 久久久久精品国产欧美久久久 | 亚洲欧美成人精品一区二区| 国产极品天堂在线| 搡老乐熟女国产| 在线观看www视频免费| 男女床上黄色一级片免费看| 成人手机av| 国产av国产精品国产| videos熟女内射| 亚洲人成电影观看| 国产人伦9x9x在线观看| 成年人午夜在线观看视频| 国产深夜福利视频在线观看| 欧美最新免费一区二区三区| 只有这里有精品99| 日日爽夜夜爽网站| 18禁观看日本| 国产福利在线免费观看视频| 大码成人一级视频| 丝袜脚勾引网站| 国产一卡二卡三卡精品 | 午夜激情久久久久久久| 午夜激情av网站| 国产午夜精品一二区理论片| 久久性视频一级片| 亚洲精品久久午夜乱码| 国产一级毛片在线| 90打野战视频偷拍视频| 精品一区二区三区四区五区乱码 | 精品国产乱码久久久久久小说| 久久久久久久精品精品| 久久久久久久大尺度免费视频| 看免费成人av毛片| 亚洲av成人精品一二三区| 美国免费a级毛片| 久久久精品免费免费高清| 1024视频免费在线观看| 丰满乱子伦码专区| 侵犯人妻中文字幕一二三四区| 久久久久精品人妻al黑| 国产精品秋霞免费鲁丝片| 国产精品熟女久久久久浪| 秋霞伦理黄片| 这个男人来自地球电影免费观看 | 午夜免费男女啪啪视频观看| 国产一区二区三区av在线| 女人爽到高潮嗷嗷叫在线视频| 亚洲国产精品国产精品| 视频在线观看一区二区三区| 91精品伊人久久大香线蕉| 美女午夜性视频免费| 中文欧美无线码| 精品一区在线观看国产| 亚洲图色成人| 伊人亚洲综合成人网| 久久 成人 亚洲| 巨乳人妻的诱惑在线观看| 欧美精品人与动牲交sv欧美| 久久影院123| 精品亚洲成国产av| 国产免费福利视频在线观看| 国产精品三级大全| 黄色 视频免费看| 国产女主播在线喷水免费视频网站| 精品第一国产精品| 一级a爱视频在线免费观看| 亚洲欧美一区二区三区久久| 久久婷婷青草| 国产av精品麻豆| 国产成人系列免费观看| 久久精品亚洲av国产电影网| 中文字幕精品免费在线观看视频| 看免费成人av毛片| 一级毛片电影观看| 久久人人97超碰香蕉20202| 久久久亚洲精品成人影院| 日韩欧美精品免费久久| 天天影视国产精品| 一区在线观看完整版| 久久久久精品性色| 午夜福利免费观看在线| 精品一区在线观看国产| 日本欧美视频一区| 色婷婷av一区二区三区视频| 婷婷成人精品国产| 1024视频免费在线观看| 五月天丁香电影| 人人妻人人澡人人爽人人夜夜| 久久综合国产亚洲精品| 美女中出高潮动态图| 久久这里只有精品19| 青春草国产在线视频| 亚洲精品aⅴ在线观看| 久久ye,这里只有精品| 国产成人欧美| 午夜福利免费观看在线| 精品一区二区免费观看| 别揉我奶头~嗯~啊~动态视频 | 亚洲第一av免费看| 中文天堂在线官网| 久久久久久久精品精品| 国产在线一区二区三区精| 婷婷色麻豆天堂久久| 人妻人人澡人人爽人人| 国产伦理片在线播放av一区| 国产成人欧美在线观看 | 性少妇av在线| 不卡av一区二区三区| 捣出白浆h1v1| 90打野战视频偷拍视频| 成年美女黄网站色视频大全免费| 午夜免费男女啪啪视频观看| 青春草国产在线视频| 亚洲五月色婷婷综合| 国产精品熟女久久久久浪| 亚洲欧洲日产国产| 啦啦啦在线免费观看视频4| av线在线观看网站| 亚洲av国产av综合av卡| 一区二区三区四区激情视频| 欧美变态另类bdsm刘玥| 极品人妻少妇av视频| 黄网站色视频无遮挡免费观看| 国产精品av久久久久免费| 色视频在线一区二区三区| 男女无遮挡免费网站观看| 欧美日韩福利视频一区二区| 欧美精品亚洲一区二区| 男女床上黄色一级片免费看| 免费不卡黄色视频| 国产 精品1| 伦理电影大哥的女人| 免费久久久久久久精品成人欧美视频| 久久久久久免费高清国产稀缺| 美女脱内裤让男人舔精品视频| 成人漫画全彩无遮挡| 极品少妇高潮喷水抽搐| 人成视频在线观看免费观看| 亚洲色图综合在线观看| 国产乱人偷精品视频| 亚洲精品日韩在线中文字幕| 日韩欧美精品免费久久| kizo精华| 亚洲免费av在线视频| 国产精品一区二区精品视频观看| 亚洲成av片中文字幕在线观看| 亚洲国产欧美在线一区| 精品国产超薄肉色丝袜足j| 久久精品国产亚洲av高清一级| 青青草视频在线视频观看| 精品一区二区三卡| 久久久精品国产亚洲av高清涩受| 一边摸一边抽搐一进一出视频| 欧美最新免费一区二区三区| 丰满乱子伦码专区| 久久久国产欧美日韩av| 男女边摸边吃奶| 丰满少妇做爰视频| 国产成人午夜福利电影在线观看| 18禁裸乳无遮挡动漫免费视频| 亚洲精品国产区一区二| 日韩人妻精品一区2区三区| 一本久久精品| 黄色视频不卡| 伦理电影免费视频| 日本91视频免费播放| 国产麻豆69| 悠悠久久av| 亚洲精品久久久久久婷婷小说| 婷婷色综合www| 日本黄色日本黄色录像| 精品一区在线观看国产| 国产亚洲精品第一综合不卡| 99热网站在线观看| 最新在线观看一区二区三区 | netflix在线观看网站| 亚洲成av片中文字幕在线观看| 国产精品免费大片| 乱人伦中国视频| 午夜av观看不卡| 丝袜美足系列| 国产精品.久久久| 久久97久久精品| 亚洲欧美成人精品一区二区| av网站在线播放免费| 久久久久久人妻| 中文欧美无线码| 男女边吃奶边做爰视频| 亚洲国产欧美日韩在线播放| 又大又爽又粗| 超碰成人久久| 欧美日韩亚洲综合一区二区三区_| 精品卡一卡二卡四卡免费| 欧美日本中文国产一区发布| 日本猛色少妇xxxxx猛交久久| 亚洲精品一区蜜桃| 少妇被粗大猛烈的视频| 亚洲三区欧美一区| 一级片免费观看大全| 日本欧美国产在线视频| 国产男女超爽视频在线观看| 国产一区二区 视频在线| 亚洲国产精品999| 黑人猛操日本美女一级片| 免费少妇av软件| 成年女人毛片免费观看观看9 | 日韩人妻精品一区2区三区| 美女国产高潮福利片在线看| 欧美激情高清一区二区三区 | 亚洲久久久国产精品| 中国国产av一级| 亚洲精华国产精华液的使用体验| 日韩av免费高清视频| 国产成人系列免费观看| 成人国产av品久久久| 免费看av在线观看网站| av在线观看视频网站免费| 国产精品免费视频内射| 永久免费av网站大全| 色婷婷久久久亚洲欧美| 丝袜美腿诱惑在线| 丰满少妇做爰视频| 国产免费现黄频在线看| 久久影院123| av视频免费观看在线观看| av又黄又爽大尺度在线免费看| tube8黄色片| 国产成人91sexporn| 97在线人人人人妻| 精品国产一区二区久久| 国产精品久久久久久精品古装| 十八禁人妻一区二区| 麻豆乱淫一区二区| 色精品久久人妻99蜜桃| 丝袜脚勾引网站| 交换朋友夫妻互换小说| 曰老女人黄片| 男女之事视频高清在线观看 | 大片电影免费在线观看免费| 午夜福利视频在线观看免费| e午夜精品久久久久久久| 韩国精品一区二区三区| 中文字幕人妻熟女乱码| 欧美亚洲日本最大视频资源| 日本av免费视频播放| 久久久久久免费高清国产稀缺| 十八禁人妻一区二区| 亚洲欧美日韩另类电影网站| 成年美女黄网站色视频大全免费| 日韩欧美精品免费久久| 国产av精品麻豆| 精品卡一卡二卡四卡免费| 亚洲人成77777在线视频| 久久精品国产综合久久久| 成人亚洲精品一区在线观看| 午夜日韩欧美国产| 黑丝袜美女国产一区| netflix在线观看网站| 欧美最新免费一区二区三区| 一级毛片我不卡| 中文字幕另类日韩欧美亚洲嫩草| 亚洲欧美清纯卡通| 天堂中文最新版在线下载| 在线 av 中文字幕| 男女边摸边吃奶| 欧美成人午夜精品| 国产黄频视频在线观看| 亚洲精品国产区一区二| 国产 一区精品| 女人被躁到高潮嗷嗷叫费观| 精品国产乱码久久久久久小说| 女人被躁到高潮嗷嗷叫费观| 日韩不卡一区二区三区视频在线| 亚洲伊人色综图| 亚洲成国产人片在线观看| 亚洲国产看品久久| 久久久久精品国产欧美久久久 | 可以免费在线观看a视频的电影网站 | 一二三四在线观看免费中文在| 母亲3免费完整高清在线观看| 老司机在亚洲福利影院| 又大又爽又粗| 另类精品久久| 欧美日韩视频精品一区| 亚洲国产成人一精品久久久| 美女主播在线视频| 熟妇人妻不卡中文字幕| 菩萨蛮人人尽说江南好唐韦庄| 亚洲国产成人一精品久久久| 看非洲黑人一级黄片| av在线app专区| 国产一区亚洲一区在线观看| 一个人免费看片子| av网站在线播放免费| h视频一区二区三区| www.熟女人妻精品国产| 国产不卡av网站在线观看| 精品一区二区三区av网在线观看 | 欧美av亚洲av综合av国产av | 国产日韩欧美亚洲二区| 中国三级夫妇交换| 日韩视频在线欧美| 多毛熟女@视频| 国产精品国产av在线观看| 久久综合国产亚洲精品| 在线天堂中文资源库| 国产成人a∨麻豆精品| 天天添夜夜摸| 中文字幕精品免费在线观看视频| 精品免费久久久久久久清纯 | 如日韩欧美国产精品一区二区三区| 飞空精品影院首页| 亚洲欧美成人综合另类久久久| 无遮挡黄片免费观看| av又黄又爽大尺度在线免费看| 午夜福利乱码中文字幕| √禁漫天堂资源中文www| 国产精品国产av在线观看| 日日爽夜夜爽网站| 肉色欧美久久久久久久蜜桃| 国产片特级美女逼逼视频| 汤姆久久久久久久影院中文字幕| 蜜桃国产av成人99| 久久久精品国产亚洲av高清涩受| www.自偷自拍.com| 成年女人毛片免费观看观看9 | 欧美日韩福利视频一区二区|