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

    基于GprMax2D的地下管線正演模擬

    2020-12-02 07:48趙慧冰范文娟郭超超呂高陳紅莉
    科技創(chuàng)新與應(yīng)用 2020年34期
    關(guān)鍵詞:管徑時(shí)域有限元

    趙慧冰 范文娟 郭超超 呂高 陳紅莉

    摘? 要:地質(zhì)雷達(dá)常用于地下管線的檢測(cè)。文章利用GprMax2D和Matlab軟件對(duì)地下管線進(jìn)行時(shí)域有限元差分(FDTD)法正演模擬,分別針對(duì)不同管徑、不同材質(zhì)、不同天線頻率和管線上不同覆土層厚度進(jìn)行模擬,分析對(duì)比圖像變化特征。結(jié)果表明:隨著管徑不斷增大,雷達(dá)成像出現(xiàn)管線上下界面;當(dāng)管線材質(zhì)不同時(shí),金屬管僅能出現(xiàn)管線上界面,PVC管和混凝土管成像效果類似;檢測(cè)天線頻率降低時(shí),管線圖像的分辨率逐漸降低,但檢測(cè)的深度不斷加大;同一管徑隨著覆土層厚度加大,管線上界面反射弧可清晰判斷,但下界面的反射弧逐漸模糊。通過正演模擬可為地質(zhì)雷達(dá)探測(cè)地下埋管的判讀提供一定的科學(xué)依據(jù)和指導(dǎo)。

    關(guān)鍵詞:地下埋管;地質(zhì)雷達(dá);FDTD正演模擬

    中圖分類號(hào):P631 文獻(xiàn)標(biāo)志碼:A 文章編號(hào):2095-2945(2020)34-0011-04

    Abstract: Geological radar is commonly used in the detection of underground pipelines. In this paper, GprMax2D software and Matlab software are used to simulate the forward simulation of the underground pipeline with the Finite-Difference Time-Domain (FDTD) method, respectively for different pipe diameters, materials, antenna frequencies, and thickness of the overburden layer on the pipeline. The results show that: with the increasing diameter of the pipeline, the upper and lower interfaces of the pipeline appear in the radar imaging; when the pipeline materials are different, the metallic pipeline can only appear on the pipeline interface, the imaging effects of the PVC pipe and the concrete pipe are similar; when the frequency of the detection antenna decreases, the resolution of the pipeline image decreases gradually, but the depth of the detection continues to increase; with the increase of the thickness of the soil layer in the same pipe diameter, the reflection arc of the upper interface of the pipeline can be clearly judged, but the reflection arc of the lower interface is gradually blurred. The forward simulation can provide some scientific basis and guidance for the interpretation of underground buried pipes detected by ground-penetrating radar (GPR).

    Keywords: underground pipe; geological radar; FDTD forward simulation

    隨著城市化進(jìn)程的加快,城市人口不斷增加,導(dǎo)致了城市供水問題日益突出,而城市建設(shè)的發(fā)展還帶動(dòng)了地下埋管建設(shè)的興起,這就使得用水問題更加緊張[1]。在運(yùn)行過程中,經(jīng)常會(huì)存在滲漏,造成水源的浪費(fèi),如果不檢測(cè)直接進(jìn)行開挖,可能會(huì)造成管線的破壞,檢測(cè)成本增加,因此出于施工安全以及管線保護(hù)的需要,地下管線的探測(cè)成為一項(xiàng)必不可少的工作。

    目前地下管線探測(cè)方法有許多種,例如管線探測(cè)儀、高密度電阻率和地質(zhì)雷達(dá)等等[2],這些方法探測(cè)的原理不同,但現(xiàn)已均建立了一套完整和成熟的體系,各種方法均具有不同的優(yōu)缺點(diǎn)和不同的使用條件,其中地質(zhì)雷達(dá)具有高精度、高分辨率和無破壞性的特點(diǎn),使用較為廣泛。地質(zhì)雷達(dá)進(jìn)行檢測(cè)主要是依據(jù)管線與周圍介質(zhì)材料在介電常數(shù)上存在差異,據(jù)此進(jìn)行物探方面的檢測(cè)。其中利用地質(zhì)雷達(dá)和GprMax2D進(jìn)行管線檢測(cè)和模擬研究的很多,如張鵬等[3]利用二維FDTD正演模擬,研究地質(zhì)雷達(dá)在管線檢測(cè)和模擬中可能遇到的工況,為實(shí)際檢測(cè)時(shí)圖像判斷提供理論的依據(jù);姚顯春等[4]利用地質(zhì)雷達(dá)正演模擬,根據(jù)補(bǔ)貼交界面的相對(duì)波幅值判斷管線的材質(zhì)和成像特征;黃樂藝[5]主要研究地質(zhì)雷達(dá)在市政管線滲漏檢測(cè)方面的模擬,得到了不同工況下的管線滲漏病害雷達(dá)圖譜,總結(jié)了地質(zhì)雷達(dá)圖像判識(shí)過程中滲漏區(qū)圖譜的典型特征;梁小強(qiáng)[6]應(yīng)用地質(zhì)雷達(dá)正演模擬出不同管線埋深、間距、管內(nèi)物和管材等影響因素,分析其正演剖面特征,了解管線探測(cè)與各影響參數(shù)之間的關(guān)系,對(duì)實(shí)際檢測(cè)起到指導(dǎo)作用。

    在前人的工作基礎(chǔ)之上,本文主要利用FDTD正演模擬的方法,對(duì)土體內(nèi)部存在的埋管進(jìn)行檢測(cè),由于電磁波在不同介電常數(shù)的介質(zhì)中傳播的速度會(huì)不同,以此為理論依據(jù),結(jié)合時(shí)域有限元差分法對(duì)地下埋管的情況進(jìn)行數(shù)值模擬分析,建立地下埋管模型作為研究對(duì)象,對(duì)管線的不同管徑大小、同一埋深和管徑下不同管材、同一埋深和管徑下不同檢測(cè)頻率以及同一管徑不同埋深下進(jìn)行模擬分析,比較圖像,據(jù)此獲得相應(yīng)的地質(zhì)雷達(dá)圖像和波形信息,經(jīng)過相應(yīng)的對(duì)比分析以期達(dá)到識(shí)別管線直徑、簡(jiǎn)單判斷管材和選擇合適的檢測(cè)天線頻率,對(duì)以后類似工程的檢測(cè)提供參考。

    1 GprMax正演模擬原理

    GprMax地質(zhì)雷達(dá)正演模擬是基于時(shí)域有限差分法(FDTD法)的地質(zhì)雷達(dá)信號(hào)正演模擬有限元軟件。時(shí)域有限元差分法(finite difference time domain method)是將麥克斯韋時(shí)域場(chǎng)旋度方程中的微分式用有限差分式代替,這樣使得連續(xù)性問題變?yōu)殡x散型問題,得到關(guān)于場(chǎng)分量的差分方程組,之后用相同電參量的空間網(wǎng)格來模擬被研究的物體[7]。在二維直角坐標(biāo)系中,橫電磁波(TM波)即當(dāng)磁場(chǎng)和電場(chǎng)與z的方向均無關(guān)時(shí),方程組可表示成兩個(gè)相互獨(dú)立的方程,此時(shí)的麥克斯韋方程組可寫為:

    式(1)中:Ez為z方向的電場(chǎng)強(qiáng)度;Hx和Hy分別為x和y方向的磁場(chǎng)強(qiáng)度;?滓為電導(dǎo)率;?滓m為導(dǎo)磁率;?著為介電常數(shù);?滋為磁導(dǎo)率。

    時(shí)域有限元差分法是把Yee元胞作為最小的單位差分網(wǎng)格,電場(chǎng)與磁場(chǎng)在時(shí)間上交替出現(xiàn)且彼此之間相差半個(gè)時(shí)間步長(zhǎng),可得到二維TM波的FDTD公式[8]:

    (2)

    (3)

    (4)

    其中:

    式中:?駐x和?駐y分別是Yee元胞x和y方向上的空間步長(zhǎng);?駐t為時(shí)間步長(zhǎng);n為時(shí)間步;(i,j)為節(jié)點(diǎn)坐標(biāo),其中i和j取值受模型差分情況的影響。

    2 地下埋管數(shù)值模擬方案

    本模型基于以下假設(shè)而建立:

    (1)地埋管土層、PVC管材和混凝土管材均為半無限空間連續(xù)均勻各向同性的完美介電體。

    (2)反射和折射過程均在二維平面中進(jìn)行,反射界面平滑,符合Snell原理。

    土壤結(jié)構(gòu)面較為連續(xù),忽略內(nèi)部某些非連續(xù)的情況,可將其探測(cè)范圍以內(nèi)的部分看作是連續(xù)的介電體。根據(jù)上述假設(shè),建立兩維地電模型,如圖1所示。

    圖1中,“T”“R”分別為地質(zhì)雷達(dá)的發(fā)射和接收裝置;此模型為土壤中埋管,其模型長(zhǎng)度為1m,深度為1.5m,相應(yīng)材料的屬性參數(shù)如表1所示。

    3 埋管正演模擬圖像分析

    3.1 不同管徑的地下埋管檢測(cè)圖

    針對(duì)于上述埋管模型,當(dāng)上方土層深度h均為0.7m時(shí),分別對(duì)管徑為100mm、200mm和300mm的PVC管,使用900M天線進(jìn)行掃描,通過對(duì)模型時(shí)域有限元差分法計(jì)算,得到如圖2所示的掃描結(jié)果。

    就管線管徑變化而言:當(dāng)管徑為100mm時(shí),由圖2(a)僅能判斷出管線的上界面,這是由于900M天線檢測(cè)分辨率導(dǎo)致無法判斷,但當(dāng)管徑增到200mm和300mm時(shí),由圖2(b)和圖2(c)可知,管線的上下界面均清晰可見,并且隨著管徑的不斷增大,管線上界面和下界面之間的距離也在不斷的加大,通過選擇合適上下界面的時(shí)間節(jié)點(diǎn),可以大致計(jì)算出管徑值。

    3.2 不同管材的地下埋管檢測(cè)圖

    針對(duì)于埋管模型,當(dāng)上方土層深度h均為0.7m時(shí),對(duì)直徑為200mm的PVC管、金屬管和混凝土管,使用900M天線進(jìn)行掃描,通過對(duì)模型時(shí)域有限元差分法計(jì)算,得到如圖3所示的掃描結(jié)果。

    就管線材質(zhì)變化而言:由圖2(b)可得到PVC管的成像信息,由圖3(a)和圖3(b)可知,金屬管和混凝土管的成像情況,對(duì)比分析可得,PVC管與混凝土管圖像信息較為相似,可掃描觀測(cè)出管線上界面和管線下界面,這是由于PVC與混凝土介電常數(shù)相差不大,所以造成電磁波的反射與折射相類似。但金屬管線則不同,僅能觀測(cè)到管線上界面,并且成像效果較其他兩者更為清晰,這是由于電磁波在金屬管線表面全反射造成的。當(dāng)管線為非金屬管線時(shí),可通過上下界面的時(shí)間節(jié)點(diǎn)計(jì)算管徑,但金屬管無法計(jì)算出其對(duì)應(yīng)的管徑信息。

    3.3 不同天線頻率下的地下埋管檢測(cè)圖

    針對(duì)于埋管模型,當(dāng)上方土層深度h均為0.7m時(shí),對(duì)直徑為200mm的PVC管分別使用900M天線、700M天線和400M天線進(jìn)行掃描,通過時(shí)域有限元差分法計(jì)算,得到如圖4所示的掃描結(jié)果。

    就檢測(cè)天線頻率變化而言:由圖2(b)可得到PVC管的成像信息,由圖4(a)和圖4(b)分別為700M天線掃描結(jié)果和400M天線掃描結(jié)果。通過對(duì)比可知,隨著天線檢測(cè)頻率的降低,圖像的分辨率也在降低,在900M天線檢測(cè)的情況下,管線上下界面可清晰判斷,在700M天線檢測(cè)的情況下,管線上界面可清晰判斷,但下界面成像不夠清晰,在400M天線檢測(cè)的情況下,管線上界面已變得模糊,下界面已無法判斷。因此,在檢測(cè)之前需要先進(jìn)行管線埋深的預(yù)判,由預(yù)判結(jié)果選擇合理的天線進(jìn)行檢測(cè)。

    3.4 不同覆土厚度的地下埋管檢測(cè)圖

    針對(duì)于埋管模型,埋管采用直徑為200mm的PVC管,管線上方土層厚度h分別為0.7m、0.8m和0.9m變化時(shí),使用900M天線進(jìn)行掃描,通過時(shí)域有限元差分法計(jì)算,得到如圖5所示的掃描結(jié)果。

    就管線上方土層厚度變化而言:由圖2(b)可得到h=0.7m時(shí)PVC管的成像信息,由圖5(a)和圖5(b)分別為h=0.8m和h=0.9m的圖像信息。通過對(duì)比可知,隨著管線上方土層厚度的不斷增大,均能判斷出管線的上下界面,但是在管線上方土層為0.9m時(shí),管線下界面成像較為模糊,這是由于隨著埋深的增加,電磁波在不斷的進(jìn)行衰減,埋深超過一定的深度后,900M天線檢測(cè)將無法判斷出其管線下界面,這時(shí)需改用低頻率天線進(jìn)行檢測(cè),但分辨率亦會(huì)有所降低。

    4 結(jié)論

    (1)檢測(cè)地下埋管時(shí),若雷達(dá)成像只顯示管線上界面,分析其原因可能為管徑過小、金屬管線、檢測(cè)使用天線頻率較低和埋深較大。

    (2)在進(jìn)行地下埋管檢測(cè)時(shí),應(yīng)首先進(jìn)行預(yù)判管線的埋深,這樣可以根據(jù)其深度選擇合理的天線頻率,保證在檢測(cè)到埋管的前提下,確保其清晰度。

    (3)同一頻率天線檢測(cè),埋深相同時(shí),管線管徑越大,反射弧的曲率半徑越大,成像效果越清晰,結(jié)果越容易判斷和分析,且管頂與管底部反射時(shí)間相差越大。

    (4)通過檢測(cè)圖像中管線上下界面選取合適的時(shí)間節(jié)點(diǎn),可以大致進(jìn)行管徑的計(jì)算。

    參考文獻(xiàn):

    [1]盧石坤,吳獻(xiàn)文.城市地下管線探測(cè)及其信息管理系統(tǒng)應(yīng)用分析[J].測(cè)繪通報(bào),2016(S1):156-159+164.

    [2]張軍偉,劉秉峰,李雪,等.基于GPRMax2D的地下管線精細(xì)化探測(cè)方法[J].物探與化探,2019,43(02):422-427.

    [3]張鵬,王旭東,王曉文,等.基于GPR的地下管線圖譜特征的正演研究[J].地下空間與工程學(xué)報(bào),2014,10(2):304-310.

    [4]姚顯春,閆茂,呂高,等.地質(zhì)雷達(dá)探測(cè)地下管線分類判別方法研究[J].地球物理學(xué)進(jìn)展,2018,33(04):1740-1747.

    [5]黃樂藝.城市地下供水管線滲漏探地雷達(dá)正演模擬與解釋方法應(yīng)用研究[D].北京交通大學(xué),2016.

    [6]梁小強(qiáng),楊道學(xué),張可能,等.FDTD數(shù)值模擬在GPR管線探測(cè)中的應(yīng)用[J].地球物理學(xué)進(jìn)展,2017,32(04):1803-1807.

    [7]趙慧冰,楊杰,程琳,等.混凝土內(nèi)裂隙空間位置及填充物FDTD正演分析[J].西安理工大學(xué)學(xué)報(bào),2017,33(04):450-456.

    [8]葛德彪,閆玉波.電磁波時(shí)域有限差分法[M].西安:西安電子科技大學(xué)出版社,2011.

    猜你喜歡
    管徑時(shí)域有限元
    有限元基礎(chǔ)與應(yīng)用課程專業(yè)賦能改革與實(shí)踐
    基于有限元的Q345E鋼補(bǔ)焊焊接殘余應(yīng)力的數(shù)值模擬
    將有限元分析引入材料力學(xué)組合變形的教學(xué)探索
    晚霞淺淡少年糖
    市政室外給水管道工程的設(shè)計(jì)研究
    基于MATLAB 的信號(hào)時(shí)域采樣及頻率混疊現(xiàn)象分析
    城市污水管網(wǎng)優(yōu)化設(shè)計(jì)及計(jì)算方法
    兩種常用漂浮式風(fēng)力機(jī)平臺(tái)動(dòng)態(tài)特性分析
    不同入射角風(fēng)波流海上漂浮式風(fēng)力機(jī)頻域與時(shí)域動(dòng)態(tài)特性
    城市供水方案的應(yīng)用分析
    国产亚洲精品第一综合不卡| av线在线观看网站| 国产av国产精品国产| 午夜成年电影在线免费观看| 国产熟女午夜一区二区三区| 精品国产一区二区三区久久久樱花| 亚洲五月婷婷丁香| 丰满人妻熟妇乱又伦精品不卡| 香蕉久久夜色| 少妇精品久久久久久久| 久久人妻福利社区极品人妻图片| 欧美黄色淫秽网站| 在线十欧美十亚洲十日本专区| av线在线观看网站| 精品国产乱码久久久久久男人| 18禁黄网站禁片午夜丰满| 精品欧美一区二区三区在线| 黑人猛操日本美女一级片| 大香蕉久久成人网| 老汉色∧v一级毛片| 国产免费视频播放在线视频| 三级毛片av免费| 午夜老司机福利片| 美女福利国产在线| 久热爱精品视频在线9| 久久久精品区二区三区| 午夜免费鲁丝| 久久天躁狠狠躁夜夜2o2o| 十八禁网站网址无遮挡| 另类亚洲欧美激情| 丰满饥渴人妻一区二区三| 新久久久久国产一级毛片| 性高湖久久久久久久久免费观看| 日本撒尿小便嘘嘘汇集6| 少妇猛男粗大的猛烈进出视频| 另类精品久久| 久久毛片免费看一区二区三区| 中文字幕人妻丝袜一区二区| 欧美黄色片欧美黄色片| 中文字幕av电影在线播放| 老司机福利观看| 精品国产乱码久久久久久小说| 啪啪无遮挡十八禁网站| 国产主播在线观看一区二区| 国产精品.久久久| 亚洲人成电影观看| 麻豆国产av国片精品| 欧美在线一区亚洲| 国产无遮挡羞羞视频在线观看| 国产精品欧美亚洲77777| 亚洲av日韩在线播放| 国产不卡一卡二| 日韩欧美一区视频在线观看| 色精品久久人妻99蜜桃| 最近最新中文字幕大全电影3 | 亚洲综合色网址| 色视频在线一区二区三区| 久久国产精品人妻蜜桃| 久久人妻熟女aⅴ| 1024香蕉在线观看| 久久久国产欧美日韩av| 亚洲美女黄片视频| 国产伦人伦偷精品视频| 国产精品免费大片| 亚洲中文字幕日韩| 久久天堂一区二区三区四区| 日韩视频一区二区在线观看| 国产成人精品无人区| 岛国毛片在线播放| 在线观看免费视频网站a站| av线在线观看网站| 在线观看www视频免费| av天堂在线播放| 国产精品.久久久| 国产成人一区二区三区免费视频网站| 丁香六月欧美| 亚洲欧美日韩另类电影网站| a级毛片黄视频| 久久精品国产99精品国产亚洲性色 | 久久99一区二区三区| 三级毛片av免费| 亚洲av成人一区二区三| 变态另类成人亚洲欧美熟女 | 大香蕉久久成人网| 精品人妻熟女毛片av久久网站| 嫩草影视91久久| 人人妻,人人澡人人爽秒播| 亚洲av日韩精品久久久久久密| 精品熟女少妇八av免费久了| xxxhd国产人妻xxx| 欧美 亚洲 国产 日韩一| 窝窝影院91人妻| 国产免费福利视频在线观看| 十八禁网站网址无遮挡| 精品乱码久久久久久99久播| 久久久久精品人妻al黑| 18禁国产床啪视频网站| 成人亚洲精品一区在线观看| 亚洲avbb在线观看| 国产亚洲一区二区精品| 老司机影院毛片| 成人国产av品久久久| 一边摸一边抽搐一进一出视频| 亚洲色图综合在线观看| 亚洲伊人久久精品综合| 日本av免费视频播放| 国产精品成人在线| 免费观看a级毛片全部| xxxhd国产人妻xxx| av线在线观看网站| 日韩人妻精品一区2区三区| 国产高清videossex| av免费在线观看网站| 色在线成人网| 日韩欧美一区视频在线观看| 少妇精品久久久久久久| 色在线成人网| 涩涩av久久男人的天堂| 午夜福利一区二区在线看| 成年人免费黄色播放视频| 国产av又大| 国产99久久九九免费精品| 婷婷成人精品国产| 考比视频在线观看| 日韩三级视频一区二区三区| 国产亚洲精品第一综合不卡| 丝瓜视频免费看黄片| 在线亚洲精品国产二区图片欧美| 亚洲欧美激情在线| 日本黄色视频三级网站网址 | 国产精品二区激情视频| 欧美黄色淫秽网站| 日韩中文字幕欧美一区二区| 久久亚洲精品不卡| 女性生殖器流出的白浆| 日韩三级视频一区二区三区| 亚洲国产成人一精品久久久| 777久久人妻少妇嫩草av网站| 日韩一区二区三区影片| 亚洲精品在线美女| 91大片在线观看| 午夜福利视频精品| 日韩中文字幕欧美一区二区| av欧美777| 免费观看人在逋| 欧美日韩福利视频一区二区| 十分钟在线观看高清视频www| 久久精品国产99精品国产亚洲性色 | www.精华液| 大片电影免费在线观看免费| 国产亚洲一区二区精品| 男女之事视频高清在线观看| 成人影院久久| 午夜激情av网站| 激情视频va一区二区三区| 亚洲欧美精品综合一区二区三区| 国产人伦9x9x在线观看| 黄色丝袜av网址大全| av电影中文网址| 国产日韩欧美视频二区| 亚洲av成人不卡在线观看播放网| 亚洲七黄色美女视频| 午夜成年电影在线免费观看| 91字幕亚洲| 中文字幕人妻丝袜制服| 成人国语在线视频| 亚洲色图综合在线观看| 两性午夜刺激爽爽歪歪视频在线观看 | 69精品国产乱码久久久| 精品福利永久在线观看| 18禁黄网站禁片午夜丰满| 国产成人免费观看mmmm| 50天的宝宝边吃奶边哭怎么回事| 免费日韩欧美在线观看| 99re6热这里在线精品视频| 精品一区二区三区四区五区乱码| 亚洲午夜精品一区,二区,三区| 国产野战对白在线观看| 欧美日本中文国产一区发布| 女人被躁到高潮嗷嗷叫费观| 国产精品av久久久久免费| 超色免费av| 人人澡人人妻人| 岛国在线观看网站| 国产一区有黄有色的免费视频| 日韩欧美免费精品| av又黄又爽大尺度在线免费看| 在线十欧美十亚洲十日本专区| 精品一区二区三卡| 香蕉丝袜av| 色综合欧美亚洲国产小说| 成年版毛片免费区| 另类亚洲欧美激情| 久久九九热精品免费| 天堂8中文在线网| 亚洲三区欧美一区| 免费在线观看完整版高清| 亚洲av国产av综合av卡| 精品国产乱子伦一区二区三区| 妹子高潮喷水视频| 老汉色av国产亚洲站长工具| 制服诱惑二区| 久久人人爽av亚洲精品天堂| 日本欧美视频一区| 亚洲人成电影免费在线| 在线观看66精品国产| 十八禁网站网址无遮挡| 国产区一区二久久| 久久久久久久大尺度免费视频| 性色av乱码一区二区三区2| 欧美精品av麻豆av| 亚洲国产毛片av蜜桃av| 人人妻人人澡人人看| 国产精品99久久99久久久不卡| 日韩三级视频一区二区三区| 国产亚洲欧美精品永久| 在线亚洲精品国产二区图片欧美| 高清黄色对白视频在线免费看| 久久99热这里只频精品6学生| 亚洲国产毛片av蜜桃av| 中文欧美无线码| 香蕉久久夜色| 国产三级黄色录像| 丰满迷人的少妇在线观看| 无遮挡黄片免费观看| 成人手机av| 日韩免费高清中文字幕av| 蜜桃在线观看..| 99热网站在线观看| 国产欧美日韩一区二区三区在线| 久久精品国产亚洲av高清一级| 日本av手机在线免费观看| 亚洲欧美色中文字幕在线| 欧美精品一区二区免费开放| 欧美日韩av久久| 女性生殖器流出的白浆| 我的亚洲天堂| 欧美黑人精品巨大| 菩萨蛮人人尽说江南好唐韦庄| 亚洲精品成人av观看孕妇| 丝袜美腿诱惑在线| 国产av国产精品国产| 午夜福利乱码中文字幕| 精品久久久久久电影网| 蜜桃国产av成人99| 久久精品成人免费网站| 热99国产精品久久久久久7| 日韩欧美三级三区| 一区二区av电影网| 超碰97精品在线观看| 日本a在线网址| 久久99热这里只频精品6学生| 男女无遮挡免费网站观看| 精品国产乱码久久久久久小说| 九色亚洲精品在线播放| 精品免费久久久久久久清纯 | 中文字幕高清在线视频| 日本黄色日本黄色录像| 国产在线精品亚洲第一网站| 国产精品久久久人人做人人爽| 三级毛片av免费| 日韩中文字幕欧美一区二区| 麻豆av在线久日| 成人永久免费在线观看视频 | 国产成+人综合+亚洲专区| 高清视频免费观看一区二区| 成人影院久久| av视频免费观看在线观看| aaaaa片日本免费| 国产精品久久久久久精品古装| 在线看a的网站| 在线观看免费高清a一片| 宅男免费午夜| 伦理电影免费视频| 日韩成人在线观看一区二区三区| 性少妇av在线| 成人国产一区最新在线观看| 国产黄色免费在线视频| 最近最新中文字幕大全电影3 | 欧美日韩精品网址| 亚洲七黄色美女视频| 美女福利国产在线| 国产在线免费精品| 久久人人爽av亚洲精品天堂| 宅男免费午夜| svipshipincom国产片| 大型黄色视频在线免费观看| 日韩欧美免费精品| 久久精品国产99精品国产亚洲性色 | 国产男女超爽视频在线观看| 99热网站在线观看| 搡老熟女国产l中国老女人| videos熟女内射| 亚洲熟女精品中文字幕| 精品一区二区三区视频在线观看免费 | 老熟妇乱子伦视频在线观看| 久久人妻熟女aⅴ| 国产免费视频播放在线视频| 岛国毛片在线播放| 国产精品电影一区二区三区 | 女人久久www免费人成看片| 免费在线观看黄色视频的| 我的亚洲天堂| 久久婷婷成人综合色麻豆| 精品久久久久久久毛片微露脸| 天堂8中文在线网| 日韩欧美一区视频在线观看| 国产激情久久老熟女| 新久久久久国产一级毛片| 日韩有码中文字幕| 女人被躁到高潮嗷嗷叫费观| 久久香蕉激情| 少妇粗大呻吟视频| 香蕉国产在线看| 水蜜桃什么品种好| 一级片免费观看大全| 热99国产精品久久久久久7| 露出奶头的视频| 亚洲 欧美一区二区三区| 1024香蕉在线观看| 久久午夜综合久久蜜桃| a级片在线免费高清观看视频| 午夜激情久久久久久久| 久久狼人影院| 国产精品免费大片| 国产高清激情床上av| 亚洲精品乱久久久久久| 欧美精品高潮呻吟av久久| 国产一区二区在线观看av| 老熟妇仑乱视频hdxx| 久久九九热精品免费| 女同久久另类99精品国产91| 日韩大码丰满熟妇| 国产精品亚洲一级av第二区| 日韩大码丰满熟妇| 一区二区三区激情视频| 国产成人影院久久av| 下体分泌物呈黄色| a级毛片在线看网站| 国产精品熟女久久久久浪| 少妇裸体淫交视频免费看高清 | 成年人黄色毛片网站| 亚洲五月色婷婷综合| 国产成人精品无人区| 免费在线观看视频国产中文字幕亚洲| 香蕉丝袜av| 咕卡用的链子| 久久中文字幕一级| 亚洲自偷自拍图片 自拍| 国产精品av久久久久免费| 亚洲国产欧美网| 岛国毛片在线播放| 男男h啪啪无遮挡| 欧美日本中文国产一区发布| 中文字幕av电影在线播放| 国产欧美日韩精品亚洲av| 国产精品一区二区免费欧美| 久久av网站| av天堂久久9| 国产免费av片在线观看野外av| www.精华液| 国产淫语在线视频| 男女高潮啪啪啪动态图| 欧美老熟妇乱子伦牲交| 丰满迷人的少妇在线观看| 国产精品免费一区二区三区在线 | 亚洲色图综合在线观看| 国产在线免费精品| 亚洲av片天天在线观看| 欧美日韩亚洲综合一区二区三区_| 9191精品国产免费久久| 免费观看人在逋| 又黄又粗又硬又大视频| 国产福利在线免费观看视频| 国产精品98久久久久久宅男小说| 亚洲一区二区三区欧美精品| 一边摸一边做爽爽视频免费| 国产精品 国内视频| 9热在线视频观看99| av网站免费在线观看视频| 欧美日韩av久久| 电影成人av| 成年人午夜在线观看视频| 男女下面插进去视频免费观看| 成年人午夜在线观看视频| cao死你这个sao货| 国产精品免费视频内射| 欧美日韩成人在线一区二区| 美女视频免费永久观看网站| 无人区码免费观看不卡 | 欧美另类亚洲清纯唯美| 在线永久观看黄色视频| 久久影院123| 亚洲精品一二三| 两个人看的免费小视频| 人妻一区二区av| 怎么达到女性高潮| 女人久久www免费人成看片| 精品高清国产在线一区| 日韩欧美免费精品| 老熟妇乱子伦视频在线观看| 日韩大码丰满熟妇| 亚洲国产欧美日韩在线播放| 精品一区二区三区四区五区乱码| 国产日韩欧美亚洲二区| 一级片'在线观看视频| 美女午夜性视频免费| 亚洲精品国产一区二区精华液| 一区福利在线观看| 黄色丝袜av网址大全| 亚洲第一青青草原| 国产精品电影一区二区三区 | 高清av免费在线| 欧美 日韩 精品 国产| 一边摸一边抽搐一进一小说 | 不卡av一区二区三区| 女人久久www免费人成看片| 国产成人精品久久二区二区免费| 免费观看av网站的网址| 亚洲欧美精品综合一区二区三区| 久久九九热精品免费| 国产高清视频在线播放一区| 免费不卡黄色视频| av在线播放免费不卡| 丰满人妻熟妇乱又伦精品不卡| 亚洲视频免费观看视频| 色视频在线一区二区三区| 欧美激情极品国产一区二区三区| 欧美乱码精品一区二区三区| 高清黄色对白视频在线免费看| 亚洲欧美精品综合一区二区三区| 如日韩欧美国产精品一区二区三区| 国产成人欧美| 国产成人免费观看mmmm| h视频一区二区三区| 国产高清videossex| 99精品欧美一区二区三区四区| 窝窝影院91人妻| 日韩视频一区二区在线观看| 91精品国产国语对白视频| 亚洲精品一二三| 精品欧美一区二区三区在线| av线在线观看网站| 中文字幕人妻熟女乱码| 麻豆成人av在线观看| 久久九九热精品免费| 狠狠精品人妻久久久久久综合| 免费不卡黄色视频| 午夜福利乱码中文字幕| 亚洲欧美日韩高清在线视频 | 亚洲视频免费观看视频| 国产国语露脸激情在线看| 日本五十路高清| 久久久久视频综合| 另类亚洲欧美激情| 免费在线观看日本一区| 老司机影院毛片| 欧美精品啪啪一区二区三区| 十八禁人妻一区二区| 国产高清视频在线播放一区| 首页视频小说图片口味搜索| 国产av一区二区精品久久| 淫妇啪啪啪对白视频| 高清视频免费观看一区二区| 高清在线国产一区| 国产主播在线观看一区二区| 一本色道久久久久久精品综合| 亚洲全国av大片| 一区二区三区精品91| 久久国产精品影院| 欧美在线一区亚洲| 亚洲午夜精品一区,二区,三区| 汤姆久久久久久久影院中文字幕| 91字幕亚洲| 老司机午夜福利在线观看视频 | 性少妇av在线| 精品人妻在线不人妻| 久久热在线av| 91成年电影在线观看| 欧美精品啪啪一区二区三区| 国产精品偷伦视频观看了| 狠狠精品人妻久久久久久综合| 国产单亲对白刺激| 国产不卡av网站在线观看| 日日爽夜夜爽网站| 男女免费视频国产| 成人特级黄色片久久久久久久 | 精品国产超薄肉色丝袜足j| 国产成人一区二区三区免费视频网站| 久久狼人影院| 脱女人内裤的视频| 亚洲av第一区精品v没综合| 在线观看免费视频日本深夜| 99热国产这里只有精品6| 国产成人系列免费观看| 久久ye,这里只有精品| 五月开心婷婷网| 免费在线观看视频国产中文字幕亚洲| 亚洲人成电影观看| av天堂在线播放| 法律面前人人平等表现在哪些方面| 精品少妇黑人巨大在线播放| 天天躁日日躁夜夜躁夜夜| 色精品久久人妻99蜜桃| 91字幕亚洲| 2018国产大陆天天弄谢| 日本精品一区二区三区蜜桃| 国产成+人综合+亚洲专区| 免费看十八禁软件| 18禁观看日本| 亚洲精品国产一区二区精华液| www.999成人在线观看| 黑人操中国人逼视频| 最近最新免费中文字幕在线| 久久人人97超碰香蕉20202| 一级片'在线观看视频| 精品少妇黑人巨大在线播放| av不卡在线播放| 久久九九热精品免费| videos熟女内射| 国产高清国产精品国产三级| 国产男女超爽视频在线观看| 侵犯人妻中文字幕一二三四区| 国产av又大| 亚洲欧美色中文字幕在线| 两个人看的免费小视频| 久久中文字幕人妻熟女| 黑人巨大精品欧美一区二区mp4| 日韩有码中文字幕| 成年动漫av网址| 中文字幕av电影在线播放| 99国产精品一区二区蜜桃av | 19禁男女啪啪无遮挡网站| 免费日韩欧美在线观看| 亚洲天堂av无毛| 国产成人系列免费观看| 色视频在线一区二区三区| 中文字幕人妻丝袜制服| 欧美日韩黄片免| 18禁裸乳无遮挡动漫免费视频| 久久av网站| 成人永久免费在线观看视频 | 99国产极品粉嫩在线观看| www.熟女人妻精品国产| 天天添夜夜摸| 亚洲精品成人av观看孕妇| 亚洲伊人久久精品综合| 国产91精品成人一区二区三区 | 中亚洲国语对白在线视频| 日韩免费av在线播放| 午夜久久久在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲专区中文字幕在线| 夫妻午夜视频| 国产一区二区三区在线臀色熟女 | 中文字幕精品免费在线观看视频| 成人手机av| 在线观看一区二区三区激情| 91老司机精品| 国产有黄有色有爽视频| 成年版毛片免费区| 亚洲人成电影观看| 亚洲国产成人一精品久久久| 国产97色在线日韩免费| 一区二区日韩欧美中文字幕| 50天的宝宝边吃奶边哭怎么回事| 91成年电影在线观看| 久久精品亚洲熟妇少妇任你| 久久亚洲真实| 欧美亚洲日本最大视频资源| 国产成人av教育| 波多野结衣av一区二区av| 美女扒开内裤让男人捅视频| 黄片小视频在线播放| 午夜福利视频精品| 亚洲九九香蕉| 九色亚洲精品在线播放| tocl精华| 韩国精品一区二区三区| 国产区一区二久久| av欧美777| 91老司机精品| 亚洲av电影在线进入| 免费日韩欧美在线观看| 一边摸一边抽搐一进一小说 | av天堂在线播放| 久久久久国产一级毛片高清牌| avwww免费| 中文字幕最新亚洲高清| 人人妻人人澡人人爽人人夜夜| 亚洲av第一区精品v没综合| 中文字幕高清在线视频| 蜜桃在线观看..| 国精品久久久久久国模美| 亚洲九九香蕉| 激情在线观看视频在线高清 | 亚洲一码二码三码区别大吗| 变态另类成人亚洲欧美熟女 | 人妻久久中文字幕网| 国产不卡一卡二| 国产精品一区二区在线不卡| 成年人免费黄色播放视频| 99久久99久久久精品蜜桃| 中文字幕人妻丝袜制服| 国产精品98久久久久久宅男小说| 亚洲精品自拍成人| 久久婷婷成人综合色麻豆| 中文字幕人妻熟女乱码| 久久人人97超碰香蕉20202| 亚洲国产av影院在线观看| 国产精品麻豆人妻色哟哟久久| 一本色道久久久久久精品综合| av在线播放免费不卡| 最近最新中文字幕大全电影3 | 国产亚洲午夜精品一区二区久久|