• <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á)
    久久久精品94久久精品| 欧美一级毛片孕妇| 麻豆成人av在线观看| 丰满迷人的少妇在线观看| 国产成人av教育| 国产成人精品在线电影| 午夜视频精品福利| 美女国产高潮福利片在线看| 亚洲精品久久午夜乱码| 夜夜夜夜夜久久久久| 黄频高清免费视频| 黄色视频不卡| 国产精品麻豆人妻色哟哟久久| 国产高清视频在线播放一区| 国产亚洲av高清不卡| 成年人午夜在线观看视频| 成人影院久久| 在线亚洲精品国产二区图片欧美| 免费av中文字幕在线| 性少妇av在线| 超色免费av| 美女午夜性视频免费| 五月天丁香电影| 精品久久久精品久久久| 亚洲av成人一区二区三| 亚洲伊人久久精品综合| 欧美激情极品国产一区二区三区| 欧美 亚洲 国产 日韩一| 超色免费av| 国产91精品成人一区二区三区 | 王馨瑶露胸无遮挡在线观看| 日韩一卡2卡3卡4卡2021年| 成人手机av| 国产免费福利视频在线观看| 精品乱码久久久久久99久播| 日韩成人在线观看一区二区三区| 在线观看66精品国产| 久久热在线av| 国产成人av教育| a在线观看视频网站| 日韩人妻精品一区2区三区| 在线观看www视频免费| 国产一区二区 视频在线| 免费av中文字幕在线| 亚洲一区二区三区欧美精品| 亚洲av第一区精品v没综合| 老汉色av国产亚洲站长工具| 亚洲第一av免费看| 美国免费a级毛片| 免费久久久久久久精品成人欧美视频| 在线av久久热| 国产精品99久久99久久久不卡| 亚洲视频免费观看视频| 国产视频一区二区在线看| 亚洲欧美一区二区三区久久| 国产一区二区 视频在线| 一本色道久久久久久精品综合| 热99久久久久精品小说推荐| 在线 av 中文字幕| 久久久久精品国产欧美久久久| 十八禁人妻一区二区| 国产xxxxx性猛交| 99精品在免费线老司机午夜| 90打野战视频偷拍视频| 老熟女久久久| 国产视频一区二区在线看| 成人永久免费在线观看视频 | 亚洲全国av大片| 国产一区二区三区视频了| 少妇猛男粗大的猛烈进出视频| 国产精品 欧美亚洲| 在线av久久热| 亚洲天堂av无毛| 丝袜人妻中文字幕| 亚洲一区二区三区欧美精品| 国产精品一区二区在线观看99| 1024香蕉在线观看| 1024香蕉在线观看| 久久中文看片网| 亚洲少妇的诱惑av| av免费在线观看网站| 亚洲av片天天在线观看| 亚洲专区字幕在线| netflix在线观看网站| 亚洲欧洲精品一区二区精品久久久| 亚洲avbb在线观看| 夫妻午夜视频| 国产精品 国内视频| 国产在线视频一区二区| 超碰97精品在线观看| 色综合欧美亚洲国产小说| 俄罗斯特黄特色一大片| 国产成人系列免费观看| 99九九在线精品视频| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲av第一区精品v没综合| 精品免费久久久久久久清纯 | 亚洲精品成人av观看孕妇| 成年人午夜在线观看视频| 免费在线观看视频国产中文字幕亚洲| 亚洲av片天天在线观看| 欧美中文综合在线视频| 丰满人妻熟妇乱又伦精品不卡| 国产高清国产精品国产三级| 欧美大码av| 男人操女人黄网站| 后天国语完整版免费观看| 美女主播在线视频| 后天国语完整版免费观看| 日本wwww免费看| 69av精品久久久久久 | 国产成人啪精品午夜网站| 在线av久久热| 美女扒开内裤让男人捅视频| 成人18禁在线播放| 久久ye,这里只有精品| cao死你这个sao货| 精品一区二区三卡| 色视频在线一区二区三区| 亚洲精品中文字幕一二三四区 | 亚洲中文av在线| tube8黄色片| 12—13女人毛片做爰片一| 亚洲专区中文字幕在线| 国产成人系列免费观看| 天堂动漫精品| 精品免费久久久久久久清纯 | 国产在线视频一区二区| 男女床上黄色一级片免费看| 国产精品九九99| 中文字幕人妻丝袜制服| 国产精品av久久久久免费| 久久热在线av| 欧美久久黑人一区二区| 丰满人妻熟妇乱又伦精品不卡| 大香蕉久久网| 激情在线观看视频在线高清 | 国产在线精品亚洲第一网站| 国产精品一区二区免费欧美| 久久青草综合色| 精品少妇久久久久久888优播| 9色porny在线观看| 一级毛片女人18水好多| av免费在线观看网站| 国产成+人综合+亚洲专区| 亚洲 欧美一区二区三区| 国产在视频线精品| 香蕉久久夜色| 曰老女人黄片| 午夜激情av网站| 成人精品一区二区免费| 9191精品国产免费久久| www日本在线高清视频| 国产精品免费视频内射| 亚洲精品成人av观看孕妇| 国产精品久久久久久精品电影小说| 岛国在线观看网站| 亚洲中文字幕日韩| 手机成人av网站| 亚洲精品一卡2卡三卡4卡5卡| 每晚都被弄得嗷嗷叫到高潮| 久久久久久人人人人人| 日韩中文字幕欧美一区二区| 免费在线观看影片大全网站| 日韩欧美国产一区二区入口| 国产高清激情床上av| 亚洲视频免费观看视频| 免费观看人在逋| 亚洲国产av新网站| 国产成人av激情在线播放| 日韩免费高清中文字幕av| 午夜免费成人在线视频| 精品久久久久久久毛片微露脸| 国产成人影院久久av| 国产精品成人在线| 女警被强在线播放| 天堂8中文在线网| netflix在线观看网站| 国产精品国产高清国产av | 色视频在线一区二区三区| 高潮久久久久久久久久久不卡| 99精品在免费线老司机午夜| 又黄又粗又硬又大视频| 一区二区三区激情视频| 精品亚洲成a人片在线观看| 亚洲成av片中文字幕在线观看| 99精国产麻豆久久婷婷| 大香蕉久久成人网| 一区二区av电影网| 国产激情久久老熟女| av天堂久久9| 欧美 亚洲 国产 日韩一| 亚洲专区中文字幕在线| 亚洲男人天堂网一区| 亚洲国产欧美一区二区综合| 91麻豆av在线| 国产免费视频播放在线视频| 久久国产精品影院| 精品亚洲成a人片在线观看| 91麻豆精品激情在线观看国产 | 国产成人av激情在线播放| 久久精品亚洲精品国产色婷小说| 不卡一级毛片| 久热这里只有精品99| 国产伦人伦偷精品视频| 三级毛片av免费| 久久精品国产99精品国产亚洲性色 | 丝袜美足系列| 欧美亚洲日本最大视频资源| 一个人免费看片子| 岛国毛片在线播放| 久久精品91无色码中文字幕| 黄色片一级片一级黄色片| 老司机午夜福利在线观看视频 | 黑人操中国人逼视频| 九色亚洲精品在线播放| 久久久久久免费高清国产稀缺| 丝瓜视频免费看黄片| 亚洲成人国产一区在线观看| 纵有疾风起免费观看全集完整版| 久久九九热精品免费| e午夜精品久久久久久久| 国产片内射在线| 精品少妇久久久久久888优播| 我要看黄色一级片免费的| 久久精品亚洲精品国产色婷小说| 手机成人av网站| 午夜免费成人在线视频| 欧美激情极品国产一区二区三区| 国产有黄有色有爽视频| 精品欧美一区二区三区在线| 成年人午夜在线观看视频| 80岁老熟妇乱子伦牲交| 色精品久久人妻99蜜桃| 国产麻豆69| 国产av国产精品国产| 亚洲精品成人av观看孕妇| 久久狼人影院| 精品卡一卡二卡四卡免费| 亚洲国产精品一区二区三区在线| 日韩大码丰满熟妇| 午夜成年电影在线免费观看| 一本综合久久免费| 国产日韩欧美视频二区| 亚洲国产毛片av蜜桃av| 欧美日韩中文字幕国产精品一区二区三区 | 999久久久精品免费观看国产| 成人18禁高潮啪啪吃奶动态图| 在线观看免费高清a一片| 9191精品国产免费久久| 亚洲精品久久午夜乱码| 成人18禁高潮啪啪吃奶动态图| √禁漫天堂资源中文www| 精品第一国产精品| 久久影院123| 午夜福利视频精品| av福利片在线| 欧美日韩一级在线毛片| 香蕉丝袜av| 精品国产乱码久久久久久男人| 欧美精品高潮呻吟av久久| 亚洲美女黄片视频| 一区福利在线观看| 国产成人系列免费观看| 色在线成人网| 国产野战对白在线观看| 免费少妇av软件| avwww免费| 宅男免费午夜| 久久天躁狠狠躁夜夜2o2o| 久久国产精品人妻蜜桃| 一区二区三区国产精品乱码| 麻豆成人av在线观看| 日韩制服丝袜自拍偷拍| 国产野战对白在线观看| videos熟女内射| avwww免费| 日本黄色日本黄色录像| 日韩欧美三级三区| 黄色视频,在线免费观看| 国产伦人伦偷精品视频| 18禁美女被吸乳视频| 久久久久久亚洲精品国产蜜桃av| www.熟女人妻精品国产| 高清av免费在线| 久久ye,这里只有精品| 国产精品 欧美亚洲| 日日摸夜夜添夜夜添小说| 国产免费视频播放在线视频| 国产精品国产高清国产av | 国产精品亚洲av一区麻豆| 亚洲第一欧美日韩一区二区三区 | 午夜91福利影院| 欧美精品啪啪一区二区三区| 午夜激情av网站| 在线观看www视频免费| 999精品在线视频| 丁香欧美五月| 亚洲国产看品久久| 国产伦理片在线播放av一区| 亚洲中文av在线| 免费观看av网站的网址| 中国美女看黄片| 69精品国产乱码久久久| 狂野欧美激情性xxxx| 精品熟女少妇八av免费久了| 久久久久国产一级毛片高清牌| 国产精品免费大片| 欧美 亚洲 国产 日韩一| 我要看黄色一级片免费的| 亚洲欧美激情在线| 人人妻,人人澡人人爽秒播| 欧美+亚洲+日韩+国产| 国产真人三级小视频在线观看| 1024香蕉在线观看| 精品久久久久久久毛片微露脸| 日本av免费视频播放| 亚洲欧洲日产国产| 女同久久另类99精品国产91| 在线观看人妻少妇| 亚洲av成人不卡在线观看播放网| 国产在线观看jvid| 亚洲av美国av| 日韩熟女老妇一区二区性免费视频| 叶爱在线成人免费视频播放| 久久人妻福利社区极品人妻图片| 欧美精品一区二区免费开放| 久久香蕉激情| 91老司机精品| 久久影院123| 国产在视频线精品| 高清欧美精品videossex| 黄色 视频免费看| 精品视频人人做人人爽| 亚洲美女黄片视频| 亚洲精品久久成人aⅴ小说| 一区福利在线观看| 人人妻人人爽人人添夜夜欢视频| videosex国产| 精品一品国产午夜福利视频| 国产成人系列免费观看| 免费在线观看完整版高清| 国产精品 国内视频| 麻豆av在线久日| av网站在线播放免费| 国产成+人综合+亚洲专区| 热re99久久精品国产66热6| 天天躁夜夜躁狠狠躁躁| 国产av一区二区精品久久| 亚洲精品美女久久av网站| 国产亚洲欧美在线一区二区| 叶爱在线成人免费视频播放| 日韩人妻精品一区2区三区| 人妻久久中文字幕网| 91老司机精品| 久久久精品区二区三区| 欧美亚洲日本最大视频资源| 国产91精品成人一区二区三区 | 91大片在线观看| 精品一区二区三区视频在线观看免费 | 香蕉久久夜色| 激情视频va一区二区三区| 搡老岳熟女国产| 亚洲男人天堂网一区| 日韩人妻精品一区2区三区| aaaaa片日本免费| 免费在线观看日本一区| 国产精品国产高清国产av | 在线观看免费视频日本深夜| 国产精品久久久人人做人人爽| 亚洲精华国产精华精| 欧美乱妇无乱码| 人妻久久中文字幕网| 中文字幕av电影在线播放| 欧美 亚洲 国产 日韩一| 欧美 日韩 精品 国产| 91麻豆精品激情在线观看国产 | 色播在线永久视频| 岛国在线观看网站| 午夜精品国产一区二区电影| 精品人妻在线不人妻| 国产成人精品在线电影| 无人区码免费观看不卡 | 男人舔女人的私密视频| 精品少妇黑人巨大在线播放| 免费观看a级毛片全部| 成人18禁高潮啪啪吃奶动态图| 亚洲午夜精品一区,二区,三区| 99国产极品粉嫩在线观看| 国产欧美亚洲国产| 交换朋友夫妻互换小说| 久久久久久亚洲精品国产蜜桃av| 黄色毛片三级朝国网站| 国产无遮挡羞羞视频在线观看| 好男人电影高清在线观看| 亚洲熟妇熟女久久| 色婷婷av一区二区三区视频| 一本久久精品| 巨乳人妻的诱惑在线观看| 久久精品国产亚洲av高清一级| 久久久精品国产亚洲av高清涩受| 怎么达到女性高潮| 天堂俺去俺来也www色官网| www.999成人在线观看| 免费人妻精品一区二区三区视频| 欧美黄色淫秽网站| av国产精品久久久久影院| 免费在线观看影片大全网站| 国产一卡二卡三卡精品| 亚洲精华国产精华精| 2018国产大陆天天弄谢| 咕卡用的链子| tube8黄色片| 大陆偷拍与自拍| 亚洲免费av在线视频| 国产成+人综合+亚洲专区| 丝瓜视频免费看黄片| 岛国毛片在线播放| 如日韩欧美国产精品一区二区三区| 天堂8中文在线网| 女人被躁到高潮嗷嗷叫费观| 亚洲视频免费观看视频| 一级片免费观看大全| 午夜成年电影在线免费观看| 欧美性长视频在线观看| 国产主播在线观看一区二区| 1024视频免费在线观看| 最近最新中文字幕大全电影3 | 女人高潮潮喷娇喘18禁视频| 久久精品人人爽人人爽视色| aaaaa片日本免费| 亚洲色图综合在线观看| 狠狠精品人妻久久久久久综合| 人人妻人人爽人人添夜夜欢视频| 麻豆av在线久日| 黄片播放在线免费| 欧美成人午夜精品| 国产一卡二卡三卡精品| tocl精华| 我的亚洲天堂| 久久人人97超碰香蕉20202| 母亲3免费完整高清在线观看| 久久国产精品人妻蜜桃| 国产高清视频在线播放一区| 国产精品一区二区在线观看99| 99久久国产精品久久久| 成人精品一区二区免费| 国产精品影院久久| 国产av精品麻豆| 大码成人一级视频| 国产av国产精品国产| 高潮久久久久久久久久久不卡| 亚洲精品在线美女| 国产一区二区三区综合在线观看| 大香蕉久久成人网| 黄色a级毛片大全视频| 纯流量卡能插随身wifi吗| 久久人妻福利社区极品人妻图片| av又黄又爽大尺度在线免费看| 少妇 在线观看| 99re6热这里在线精品视频| 国产成人av教育| 在线观看舔阴道视频| 亚洲国产毛片av蜜桃av| 69av精品久久久久久 | 激情视频va一区二区三区| 亚洲av美国av| 国产国语露脸激情在线看| 热99国产精品久久久久久7| a级片在线免费高清观看视频| 亚洲成国产人片在线观看| 老熟妇仑乱视频hdxx| 久久免费观看电影| 午夜免费鲁丝| 国产野战对白在线观看| 黄色片一级片一级黄色片| 国产亚洲av高清不卡| 男女之事视频高清在线观看| 久久国产亚洲av麻豆专区| 国产主播在线观看一区二区| 三级毛片av免费| 国产精品98久久久久久宅男小说| 老司机影院毛片| 亚洲精品久久午夜乱码| 精品一品国产午夜福利视频| 日本撒尿小便嘘嘘汇集6| 亚洲熟女毛片儿| 亚洲一码二码三码区别大吗| 超碰成人久久| 免费高清在线观看日韩| 91大片在线观看| 在线观看一区二区三区激情| 视频区图区小说| 精品国产超薄肉色丝袜足j| av网站免费在线观看视频| 中文字幕高清在线视频| 欧美老熟妇乱子伦牲交| 亚洲第一欧美日韩一区二区三区 | 天天躁狠狠躁夜夜躁狠狠躁| 汤姆久久久久久久影院中文字幕| 久久国产亚洲av麻豆专区| 一个人免费看片子| 久久香蕉激情| 国产亚洲精品一区二区www | 嫩草影视91久久| 后天国语完整版免费观看| 麻豆成人av在线观看| 中亚洲国语对白在线视频| 一本久久精品| 十八禁网站网址无遮挡| 国产精品久久电影中文字幕 | 午夜91福利影院| 啦啦啦 在线观看视频| 黑人猛操日本美女一级片| 国产野战对白在线观看| 亚洲五月婷婷丁香| 午夜福利视频精品| 亚洲视频免费观看视频| 欧美一级毛片孕妇| 老熟女久久久| 亚洲人成电影免费在线| 国产一区二区三区视频了| 新久久久久国产一级毛片| 久久天躁狠狠躁夜夜2o2o| 岛国在线观看网站| av有码第一页| 老司机福利观看| 极品少妇高潮喷水抽搐| 亚洲一卡2卡3卡4卡5卡精品中文| 欧美 日韩 精品 国产| 母亲3免费完整高清在线观看| 亚洲国产欧美网| 亚洲中文av在线| 精品少妇久久久久久888优播| 久久av网站| 国产成人av激情在线播放| 精品乱码久久久久久99久播| 欧美国产精品一级二级三级| 国产日韩欧美在线精品| 免费不卡黄色视频| 亚洲美女黄片视频| 一本色道久久久久久精品综合| 欧美日韩黄片免| 纵有疾风起免费观看全集完整版| 不卡av一区二区三区| videosex国产| 国精品久久久久久国模美| 19禁男女啪啪无遮挡网站| 中文字幕av电影在线播放| 19禁男女啪啪无遮挡网站| 久久99热这里只频精品6学生| 一进一出好大好爽视频| 色综合欧美亚洲国产小说| 国产深夜福利视频在线观看| 日韩有码中文字幕| 激情视频va一区二区三区| av超薄肉色丝袜交足视频| 亚洲国产欧美网| videosex国产| 久久国产亚洲av麻豆专区| 美女高潮喷水抽搐中文字幕| 美女扒开内裤让男人捅视频| 十八禁高潮呻吟视频| 国产欧美日韩综合在线一区二区| 国产精品98久久久久久宅男小说| 国产av精品麻豆| 99精品在免费线老司机午夜| 亚洲第一青青草原| 午夜久久久在线观看| 国产又爽黄色视频| 99热网站在线观看| 亚洲国产欧美日韩在线播放| 亚洲精品自拍成人| 老熟妇乱子伦视频在线观看| 国产在线免费精品| 久久久久精品国产欧美久久久| av在线播放免费不卡| 亚洲国产欧美在线一区| 国产成人系列免费观看| 精品国产一区二区久久| 亚洲国产成人一精品久久久| 青草久久国产| 999精品在线视频| 在线观看免费午夜福利视频| 欧美日韩亚洲高清精品| 两性夫妻黄色片| 午夜两性在线视频| 欧美亚洲 丝袜 人妻 在线| 黄网站色视频无遮挡免费观看| 亚洲精品自拍成人| 国产精品久久久久久精品电影小说| 男女下面插进去视频免费观看| a级毛片黄视频| 电影成人av| 女同久久另类99精品国产91| 久久精品国产亚洲av高清一级| 中文字幕精品免费在线观看视频| av电影中文网址| 免费高清在线观看日韩| 少妇猛男粗大的猛烈进出视频| 日韩熟女老妇一区二区性免费视频| 午夜福利视频精品| 日韩视频一区二区在线观看| 亚洲国产欧美网| 女性生殖器流出的白浆| 中文字幕高清在线视频| 99国产综合亚洲精品| 免费在线观看视频国产中文字幕亚洲| 中文字幕另类日韩欧美亚洲嫩草| 9191精品国产免费久久| 欧美激情高清一区二区三区| 五月天丁香电影| 一级,二级,三级黄色视频| 女同久久另类99精品国产91| 久久精品成人免费网站|