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

    探地雷達(dá)薄層信號(hào)的譜反演算法*

    2011-11-23 06:35:02黃忠來張建中黃吉林
    關(guān)鍵詞:演算法反射面反射系數(shù)

    黃忠來 張建中 黃吉林

    (1)廈門大學(xué)信息科學(xué)與技術(shù)學(xué)院,廈門 361005 2)中國海洋大學(xué)海洋地球科學(xué)學(xué)院,青島 266100)

    探地雷達(dá)薄層信號(hào)的譜反演算法*

    黃忠來1)張建中2)黃吉林1)

    (1)廈門大學(xué)信息科學(xué)與技術(shù)學(xué)院,廈門 361005 2)中國海洋大學(xué)海洋地球科學(xué)學(xué)院,青島 266100)

    為了從多層介質(zhì)的探地雷達(dá)回波數(shù)據(jù)中提取薄層的位置和厚度信息,推導(dǎo)和驗(yàn)證了在頻率域內(nèi)實(shí)現(xiàn)的譜反演算法。該算法給出了各反射面的廣義反射系數(shù),并建立了用于聯(lián)系反射系數(shù)序列頻譜和薄層參數(shù)的代價(jià)函數(shù)并進(jìn)行反演。仿真結(jié)果表明,對(duì)于厚度小于調(diào)諧厚度的薄層,該算法也可以進(jìn)行準(zhǔn)確反演,從而提高了對(duì)于薄層的分辨能力。

    探地雷達(dá);薄層;譜反演算法;廣義反射系數(shù);代價(jià)函數(shù)

    1 引言

    探地雷達(dá) (Ground Penetrating Radar,GPR)通過向地下發(fā)射電磁脈沖,并接收地下回波信號(hào)來探測地下目標(biāo),它具有分辨率高、無損、快速等優(yōu)點(diǎn)。本文研究的對(duì)象是含有薄層的層狀結(jié)構(gòu)地下介質(zhì),目的是通過反演算法得到地下各反射面的位置以及薄層的厚度,這里的薄層是指厚度小于電磁波在其中波長 1/4的層。對(duì)于層狀介質(zhì)參數(shù)的反演,文獻(xiàn)[1]針對(duì)地震勘探共中心點(diǎn)方法 (CMP)中的模型進(jìn)行改進(jìn),并在時(shí)域反演,以求取地下兩層介質(zhì)的厚度和介電常數(shù),但是反演時(shí)需要在時(shí)域得到反射面回波的到達(dá)時(shí)間。文獻(xiàn)[2]研究路面層狀結(jié)構(gòu),利用雷達(dá)波傳播規(guī)律進(jìn)行時(shí)域反演,得到層的復(fù)介電常數(shù)和厚度,但是同樣沒有考慮目標(biāo)層是薄層的情況。本文在文獻(xiàn)[3]提出的地震勘探譜反演算法的基礎(chǔ)上,結(jié)合電磁波在層狀介質(zhì)中傳播的特點(diǎn),提出用廣義反射系數(shù)代替原算法中的反射系數(shù),推導(dǎo)并驗(yàn)證了對(duì)地下薄層 GPR信號(hào)進(jìn)行反演的頻域算法,提高了對(duì)地下薄層的分辨力。

    2 譜反演算法過程

    設(shè) re=(r1+r2)/2,ro=(r1-r2)/2,則:

    其中 Re[g(f)]=2recos(πfT),Im[g(f)]=2rosin(πfT)分別是實(shí)部和虛部。從定義式可以看出 re是函數(shù) f(r1,r2)=r1δ(t-t1)+r2δ(t-t2)做奇偶分解后的偶部分,ro則是奇部分(圖 2)。

    令 G(f)表示 g(f)的幅度譜,則:

    圖1 反射系數(shù)對(duì)Fig.1 Reflection coefficients pairs

    圖2 反射系數(shù)對(duì)奇偶分解Fig.2 Even-odd decomposition of coefficients pairs

    對(duì) G(f)取導(dǎo)數(shù)得:

    將 (3)、(4)兩式相乘得:

    找到使得代價(jià)函數(shù)值為零的參數(shù)和后,就可以根據(jù)

    求出上、下反射系數(shù)的奇偶分量,進(jìn)而得到反射系數(shù)值。

    由于事先并不知道反射系數(shù)譜 g(f),所以在計(jì)算 G(f)時(shí)要利用接收到的雷達(dá)回波和雷達(dá)發(fā)射子波的頻譜。此外,算法建立代價(jià)函數(shù)的前提是分析窗口時(shí)間零點(diǎn)位置位于兩反射系數(shù)的中點(diǎn)。下面,將在此算法的基礎(chǔ)上推導(dǎo)適用于多反射系數(shù)的代價(jià)函數(shù),本文中所推導(dǎo)出的函數(shù)表達(dá)式和文獻(xiàn)[3]中的有些不同。

    以包含 4個(gè)反射系數(shù)的模型為例,圖 3中反射系數(shù) r1~r4分別為 0.1、0.3、-0.2、0.4,t1~t4分別是 4個(gè)反射系數(shù)的距離地面的時(shí)間位置,T1、T2是前后兩對(duì)反射系數(shù)的時(shí)間厚度,則反射系數(shù)序列g(shù)(t)及其頻譜 g(f)可以表示為 g(t)=r1(t)+r2(t) =[r1δ(t-t1)+r2δ(t-t1-T1)]+[r3δ(t-t3)+r4δ (t-t3-T2)],以及 g(f)=r1(f)+r2(f)。

    圖3 反射系數(shù)對(duì)及其奇偶分解Fig.3 Coefficients pairs and their even-odd decompositions

    由前面兩反射系數(shù)的推導(dǎo)可知:

    將 4個(gè)反射系數(shù) (兩個(gè)反射系數(shù)對(duì))推廣到 n個(gè)反射系數(shù)對(duì),則代價(jià)函數(shù)可表示為:

    此代價(jià)函數(shù)含有以下未知量:反射系數(shù)對(duì)中,上層反射系數(shù)的時(shí)間位置 t,兩反射系數(shù)之間時(shí)間厚度T,反射系數(shù)對(duì)的偶分量 re以及反射系數(shù)對(duì)的奇分量 ro。以上就是譜分解算法的大致過程。通過最小化代價(jià)函數(shù)的值,可以得到地下各層的時(shí)間位置、厚度以及反射系數(shù)。但是,算法中討論的是簡單的反射系數(shù)褶積模型,未討論波在地下傳播時(shí)的衰減。對(duì)于探地雷達(dá)波來說,反演出的反射系數(shù)并不等于反射面真正的反射系數(shù),因?yàn)榈叵陆橘|(zhì)實(shí)際上是非理想的有耗介質(zhì),導(dǎo)致電磁波在傳播時(shí)會(huì)逐漸衰減,影響其傳播的除了反射系數(shù)、透射系數(shù)外還有衰減常數(shù)。因此原代價(jià)函數(shù)里的反射系數(shù)奇偶分量需要重新定義。

    3 廣義反射系數(shù)

    同理,第二個(gè)反射界面的廣義反射系數(shù)為:

    以此類推,第 i個(gè)廣義反射系數(shù)為:

    式中 zi為波在第 i層中單向傳播的距離,c為光速,并認(rèn)為α0=0。

    圖4 電磁波在地下傳播路徑示意圖Fig.4 Electromagnetic wave propagation path under ground

    4 反演算法

    對(duì)于包含多參數(shù)的代價(jià)函數(shù)的反演,本文采用在模擬退火算法上改進(jìn)得到的隨機(jī)爬山法[4]。其算法步驟為:1)產(chǎn)生初始參數(shù)向量 X={x1,x2,…, xi,…,xN},每個(gè)參數(shù)的值可以在規(guī)定的范圍內(nèi)隨機(jī)生成,也可以根據(jù)實(shí)際情況下得到的數(shù)據(jù)進(jìn)行設(shè)定。這時(shí),向量對(duì)應(yīng)的代價(jià)函數(shù)為Obj。

    2)在{x1,x2,…,xi,…,xN}中隨機(jī)選取一個(gè)參數(shù)進(jìn)行修改,例如選取的是第 i個(gè)變量 xi,則修改后xi變?yōu)榍?xi+sign*d* rang,sign是在 -1和 1兩者間隨機(jī)選取的符號(hào),d為修改的步長,用以控制反演的速度和精度,rand為 0~1之間的隨機(jī)數(shù)。

    重復(fù) 2)到 4)的操作,直到最后計(jì)算出的代價(jià)函數(shù)值滿足收斂條件,或者反演的次數(shù)超過設(shè)定值為止。每次修改時(shí)的步長 d不是固定的,開始時(shí)為了加快收斂速度,步長可以取得較大,以后逐漸減小,從而提高每次修改的精度,盡量使算法收斂到全局最優(yōu)值。

    5 數(shù)據(jù)仿真

    5.1 楔形模型仿真

    為討論目標(biāo)層厚度對(duì)算法結(jié)果的影響,利用GPRMax軟件[5]進(jìn)行正演仿真,產(chǎn)生作為實(shí)際參考值的數(shù)據(jù)。建立上、下層面反射系數(shù)同號(hào)和異號(hào)兩種情況下的楔形模型A和B(圖 5和圖 6分別是A和B的時(shí)域回波圖)。模型 A中,令楔形上界面反射系數(shù)為 -0.25,下界面為 0.25。模型 B中的楔形,令上界面反射系數(shù) -0.14,下界面為 -0.11。楔形厚度最厚處為發(fā)射波在楔形中波長的 1/2,雷達(dá)發(fā)射波中心頻率設(shè)定為 300 MHz(圖 5~8)。

    圖7和圖 8分別是楔形模型 A和 B的反射系數(shù)反演結(jié)果。從圖7和圖8可以看到,在兩種反射系數(shù)組合下的模型中,1/4波長調(diào)諧厚度均位于記錄道第 25道,當(dāng)楔形厚度逐漸減小到 1/8波長時(shí),算法仍然可以正確反演出反射面的位置和厚度,這說明了算法對(duì)于薄層的分辨效果非常好。

    圖6 模型B回波Fig.6 Received waves ofmodelB

    圖7 模型A反射系數(shù)反演結(jié)果Fig.7 Inversion results of coefficients in modelA

    圖8 模型B反射系數(shù)反演結(jié)果Fig.8 Inversion results of coefficients in modelB

    5.2 多反射系數(shù)層狀模型

    設(shè)計(jì)模型在地下有 4個(gè)反射界面,反射系數(shù)分別是 r1=-0.250 0、r2=0.111 1、r3=0.200 0、r4= 0.500 0。其相應(yīng)的廣義反射系數(shù)分別是 gr1= -0.113 1、gr2=0.043 7、gr3= -0.030 6、gr4= 0.069 9。地下包含了 2個(gè)層,上方層的厚度設(shè)計(jì)成電磁波在層中波長的 1/10,下方層的厚度設(shè)計(jì)成波長的 1/2,即上方是個(gè)相對(duì)薄層,而下方是相對(duì)厚層。

    表1 實(shí)際值與反演結(jié)果比較Tab.1 Comparison between i nversion results and true values

    圖9 回波、發(fā)射波及反射系數(shù)序列頻譜Fig.9 Spectrum of reflected wave,incident wave and reflection coefficients

    圖10 反演結(jié)果與實(shí)際值比較Fig.10 Comparison between inversion results and true values

    圖9是雷達(dá)回波、發(fā)射波的頻譜以及由前兩者計(jì)算出的模型反射系數(shù)序列頻譜,圖 10是由反演出的參數(shù)計(jì)算得到的頻譜和真實(shí)反射系數(shù)序列頻譜之間的對(duì)比,上方是兩個(gè)頻譜的實(shí)部,下方是虛部,其中紅色星劃線是真實(shí)值,藍(lán)色圈劃線是反演值。從表 1的多個(gè)反射面模型反演結(jié)果可以看到,薄層的位置以及厚度的反演結(jié)果準(zhǔn)確,表明當(dāng)?shù)叵聻槎鄬咏Y(jié)構(gòu)時(shí),算法仍然正確有效。其中薄層的廣義反射系數(shù)反演值與真實(shí)值存在一定誤差。

    6 結(jié)語

    通過研究電磁波在地下層狀介質(zhì)中的傳播規(guī)律,重新定義了反射系數(shù)表達(dá)式,推導(dǎo)并驗(yàn)證了適用于 GPR信號(hào)的譜反演算法。正是因?yàn)橛辛诵碌膹V義反射系數(shù)定義式,才可以根據(jù)實(shí)際模型參數(shù)來驗(yàn)證算法反演結(jié)果的正確性。利用譜反演算法可以得到厚度小于 1/8波長的薄層的反射面位置,層厚以及每個(gè)反射面廣義反射系數(shù)的大小。仿真結(jié)果表明該算法對(duì)于提高薄層的分辨有很好的效果,但是在優(yōu)化算法的結(jié)果上存在多解性問題,有進(jìn)一步改進(jìn)的余地。

    1 Kao Chienping,et al.Measurement of layer thickness and permittivity using a new multilayer model from GPR data [J].IEEE Transactions on geoscience and remote sensing, 2007,45(8):2 463-2 470.

    2 張蓓.路面結(jié)構(gòu)層材料介電特性及其厚度反演分析的系統(tǒng)識(shí)別方法——路面雷達(dá)關(guān)鍵技術(shù)研究 [D].重慶大學(xué),2003. (Zhang Bei.System identification method for back-calculating the dielectric property and thickness of pavement structures_study on applied technology of ground penetrating radar[D].ChongqingUniversity,2003)

    3 Charles I Puryear and John P Castagna.Layer-thickness deter mination and stratigraphic interpretation using spectral inversion:Theory and application[J].Geophysics,2008, 73(2):37-48.

    4 張霖斌,紀(jì)晨,姚振興.疊后地震道反演的隨機(jī)爬山法[J].石油地球物理勘探,1997,32(1):75-80.(Zhang Linbin,JiChen and Yao Zhenxing.Stochastic hill-climbing algorithm of post-stack seismic inversion[J].Oil Geophysical Prospecting,1997,32(1):75-80)

    5 Giannopoulos A.Modelling ground penetrating radar by Gpr Max[J].Construction and Building Materials,2005, 19:755-762.

    6 秦德文.基于譜反演的薄層預(yù)測與反演方法研究[D].中國石油大學(xué),2009.(Qin Dewen.Thin bed prediction using spectral inversion and inversion method analysis[D].China University of Petroleum,2009)

    7 Liu C R,et al.New model for estimating the thickness and per mittivity of subsurface layers from GPR data[J]. IEE Proe.-Rodur SonorNaviz,2002,149(6):315-319.

    A SPECTRAL INVERSI ON ALGORITHM FOR GPR SIGNALS OF THIN LAYERS

    Huang Zhonglai1),Zhang Jianzhong2)and Huang Jilin1)

    (1)School of Infor m ation Science and Technology,X iam en University,X iam en 361005 2)College of M arine Geo-science,Ocean University of China,Q ingdao 266100)

    In order to extract infor mation including thin layer’s position and thickness,from reflected wave data of Ground Penetrating Radar(GPR)in multi-layered situation,a spectral inversion algorithm realized in frequency domain is presented.Through analysis of electromagnetic wave propagation in multi-layermedia,generalized reflection coefficients are derived.After dividing the coefficients into pairs followed by decomposing the pairs into even and odd components,a cost function is established to connect reflection coefficients’spectrum and layer’sparameters before inversion is carried out.Some simulations show that accurate inversion results can still be expected when thickness of the layer is less than tuning thickness,which therefore i mproves resolution of thin layers.

    Ground Penetrating Radar;thin layer(GPR);spectral inversion algorithm;general reflection coefficients;cost function

    waves ofmodelA

    1671-5942(2011)04-0154-06

    2011-03-16

    國家自然科學(xué)基金(40774065)

    黃忠來,博士研究生,主要研究方向?yàn)樘降乩走_(dá)及隨機(jī)信號(hào)處理.E-mail:huangzhonglai@163.com

    P631.2

    A

    猜你喜歡
    演算法反射面反射系數(shù)
    智能反射面輔助通信中的信道估計(jì)方法
    超大規(guī)模智能反射面輔助的近場移動(dòng)通信研究
    智能反射面輔助的覆蓋增強(qiáng)技術(shù)綜述
    《四庫全書總目》子部天文演算法、術(shù)數(shù)類提要獻(xiàn)疑
    國學(xué)(2021年0期)2022-01-18 05:59:08
    單多普勒天氣雷達(dá)非對(duì)稱VAP風(fēng)場反演算法
    多道隨機(jī)稀疏反射系數(shù)反演
    石油物探(2020年6期)2020-11-25 02:38:46
    球面波PP反射系數(shù)的頻變特征研究
    運(yùn)動(dòng)平臺(tái)下X波段雷達(dá)海面風(fēng)向反演算法
    一種重新賦形副反射面的環(huán)焦天線設(shè)計(jì)
    沙質(zhì)沉積物反射系數(shù)的寬帶測量方法
    国产精华一区二区三区| 国产一区二区在线观看日韩| 成人国产综合亚洲| 成人精品一区二区免费| 大又大粗又爽又黄少妇毛片口| 成人无遮挡网站| 女生性感内裤真人,穿戴方法视频| 亚洲一级一片aⅴ在线观看| 久久精品久久久久久噜噜老黄 | 欧美日韩黄片免| 国产视频一区二区在线看| h日本视频在线播放| 天天一区二区日本电影三级| 在线观看舔阴道视频| 亚洲最大成人中文| 色综合婷婷激情| 国产探花在线观看一区二区| 免费观看精品视频网站| 伊人久久精品亚洲午夜| 午夜免费成人在线视频| 极品教师在线免费播放| 日韩 亚洲 欧美在线| 黄色丝袜av网址大全| 国产高清三级在线| 亚洲av二区三区四区| 老司机深夜福利视频在线观看| 日本一本二区三区精品| 国产在线精品亚洲第一网站| 成人美女网站在线观看视频| 中文资源天堂在线| 日韩精品青青久久久久久| 丰满人妻一区二区三区视频av| 极品教师在线免费播放| 桃色一区二区三区在线观看| 日日夜夜操网爽| av福利片在线观看| 日韩人妻高清精品专区| 亚洲熟妇熟女久久| 99久久精品国产国产毛片| 校园春色视频在线观看| 天天一区二区日本电影三级| 日韩欧美国产一区二区入口| 亚洲熟妇中文字幕五十中出| 国产真实乱freesex| 在线播放国产精品三级| 少妇熟女aⅴ在线视频| 久久国产乱子免费精品| 色哟哟·www| 午夜福利在线观看吧| 99热这里只有精品一区| 亚洲最大成人中文| 日韩中字成人| 久久久久久国产a免费观看| 99精品久久久久人妻精品| 十八禁国产超污无遮挡网站| 成人高潮视频无遮挡免费网站| 国产av在哪里看| 久久人人爽人人爽人人片va| 精品日产1卡2卡| 免费在线观看日本一区| 久久久久国产精品人妻aⅴ院| 国产女主播在线喷水免费视频网站 | 成人特级黄色片久久久久久久| 天天躁日日操中文字幕| 国产精品国产三级国产av玫瑰| 亚洲美女黄片视频| 日本 av在线| 亚洲精品成人久久久久久| 亚洲成人中文字幕在线播放| 97人妻精品一区二区三区麻豆| 日韩中字成人| 亚洲一级一片aⅴ在线观看| 亚洲精品国产成人久久av| 欧美日韩乱码在线| 午夜福利视频1000在线观看| 国产 一区 欧美 日韩| 看免费成人av毛片| 日日啪夜夜撸| 国产精品自产拍在线观看55亚洲| 一进一出抽搐动态| 国产午夜精品论理片| 又粗又爽又猛毛片免费看| 国产高清视频在线观看网站| 丝袜美腿在线中文| 内地一区二区视频在线| 两性午夜刺激爽爽歪歪视频在线观看| 深夜a级毛片| 亚洲av熟女| 久久精品国产99精品国产亚洲性色| 午夜日韩欧美国产| 久久人人精品亚洲av| 有码 亚洲区| 欧美性猛交黑人性爽| 长腿黑丝高跟| 亚洲欧美日韩东京热| 欧美+日韩+精品| 色综合色国产| 欧美日本视频| 亚洲av五月六月丁香网| 少妇熟女aⅴ在线视频| 精品一区二区三区视频在线观看免费| 1000部很黄的大片| 久久久午夜欧美精品| 搡老岳熟女国产| 不卡一级毛片| 中文字幕av成人在线电影| 成人av一区二区三区在线看| 麻豆成人午夜福利视频| 高清在线国产一区| 国产精品电影一区二区三区| 看免费成人av毛片| 国产精品国产高清国产av| 欧美区成人在线视频| 网址你懂的国产日韩在线| 五月伊人婷婷丁香| 99热只有精品国产| 免费不卡的大黄色大毛片视频在线观看 | 国产av麻豆久久久久久久| 中国美白少妇内射xxxbb| 日本-黄色视频高清免费观看| 99久久久亚洲精品蜜臀av| 婷婷精品国产亚洲av在线| 深夜a级毛片| 亚洲综合色惰| 久久精品国产亚洲av香蕉五月| 日本精品一区二区三区蜜桃| 中文字幕熟女人妻在线| 色精品久久人妻99蜜桃| 性色avwww在线观看| 乱系列少妇在线播放| 久久久国产成人精品二区| 男人的好看免费观看在线视频| 少妇被粗大猛烈的视频| 国产伦精品一区二区三区视频9| 91麻豆精品激情在线观看国产| 国产伦精品一区二区三区视频9| 国内毛片毛片毛片毛片毛片| 国产精品,欧美在线| 综合色av麻豆| 一区二区三区四区激情视频 | 免费无遮挡裸体视频| 欧美区成人在线视频| 日本与韩国留学比较| www日本黄色视频网| 欧美bdsm另类| 少妇丰满av| 精品久久久久久,| 亚洲男人的天堂狠狠| 欧美色视频一区免费| 一本久久中文字幕| 极品教师在线视频| 亚洲人成网站高清观看| 一级av片app| 简卡轻食公司| 乱系列少妇在线播放| 亚洲一级一片aⅴ在线观看| 日韩欧美免费精品| 女人十人毛片免费观看3o分钟| 1024手机看黄色片| 人妻久久中文字幕网| 真实男女啪啪啪动态图| 老司机深夜福利视频在线观看| 国产女主播在线喷水免费视频网站 | 伦理电影大哥的女人| 成年版毛片免费区| 午夜爱爱视频在线播放| 午夜精品久久久久久毛片777| 99精品久久久久人妻精品| 如何舔出高潮| 午夜视频国产福利| 国产乱人伦免费视频| 欧美又色又爽又黄视频| 一区福利在线观看| 床上黄色一级片| 一夜夜www| 一级黄片播放器| 国产 一区 欧美 日韩| 中文字幕高清在线视频| av福利片在线观看| 91午夜精品亚洲一区二区三区 | 观看免费一级毛片| 中文字幕av成人在线电影| 最近视频中文字幕2019在线8| 免费黄网站久久成人精品| 在线观看66精品国产| 欧美黑人欧美精品刺激| 乱系列少妇在线播放| 亚洲国产日韩欧美精品在线观看| 搡老岳熟女国产| 99久久精品国产国产毛片| 国产精品电影一区二区三区| 国产精品无大码| 久久久久久久精品吃奶| 免费人成在线观看视频色| 精品久久久久久久久亚洲 | 成人高潮视频无遮挡免费网站| 日韩欧美三级三区| 色视频www国产| 又黄又爽又免费观看的视频| 国产成年人精品一区二区| 成年免费大片在线观看| 日韩精品中文字幕看吧| 日韩欧美三级三区| 性欧美人与动物交配| 免费观看在线日韩| 偷拍熟女少妇极品色| 亚洲av不卡在线观看| 午夜福利在线观看免费完整高清在 | 免费黄网站久久成人精品| 99视频精品全部免费 在线| 哪里可以看免费的av片| 日本在线视频免费播放| 日本黄大片高清| 极品教师在线免费播放| 一本精品99久久精品77| av.在线天堂| 99在线人妻在线中文字幕| 99热这里只有是精品在线观看| 日韩精品中文字幕看吧| 十八禁网站免费在线| 国产69精品久久久久777片| 国产乱人伦免费视频| 久久久久久久久久黄片| 亚洲一区高清亚洲精品| 黄色日韩在线| 亚洲精品国产成人久久av| 午夜福利成人在线免费观看| 成年女人看的毛片在线观看| 国产精品,欧美在线| 无人区码免费观看不卡| 亚洲精品一卡2卡三卡4卡5卡| 精品午夜福利在线看| 欧美3d第一页| 久久香蕉精品热| 变态另类丝袜制服| 男人狂女人下面高潮的视频| 黄色配什么色好看| 国产精品日韩av在线免费观看| 中出人妻视频一区二区| 老熟妇仑乱视频hdxx| 国内毛片毛片毛片毛片毛片| 国产激情偷乱视频一区二区| 黄片wwwwww| 国产亚洲精品久久久com| 国产亚洲精品综合一区在线观看| 婷婷六月久久综合丁香| 久99久视频精品免费| 亚洲最大成人中文| 国产精品三级大全| 又爽又黄无遮挡网站| 熟女电影av网| 亚洲狠狠婷婷综合久久图片| 久久精品国产亚洲av天美| 欧美一区二区精品小视频在线| av在线天堂中文字幕| 免费黄网站久久成人精品| 国产女主播在线喷水免费视频网站 | 五月伊人婷婷丁香| 搡老妇女老女人老熟妇| 亚洲欧美日韩高清在线视频| 一级a爱片免费观看的视频| 午夜老司机福利剧场| 麻豆一二三区av精品| 欧美日韩精品成人综合77777| 俺也久久电影网| 国产精品不卡视频一区二区| 久久精品夜夜夜夜夜久久蜜豆| 免费观看人在逋| 国产真实乱freesex| 久久久久精品国产欧美久久久| 日本熟妇午夜| av中文乱码字幕在线| 亚洲精品粉嫩美女一区| 97热精品久久久久久| 亚洲精品456在线播放app | 国产三级在线视频| 亚洲精品粉嫩美女一区| av女优亚洲男人天堂| 亚洲美女搞黄在线观看 | 亚洲一区高清亚洲精品| 真人做人爱边吃奶动态| 欧美三级亚洲精品| 日日啪夜夜撸| a级毛片a级免费在线| 国产美女午夜福利| 精品久久久久久久久久免费视频| 波野结衣二区三区在线| 日本 av在线| 欧美日韩乱码在线| 一个人观看的视频www高清免费观看| 欧美丝袜亚洲另类 | 一a级毛片在线观看| 亚洲狠狠婷婷综合久久图片| 男人和女人高潮做爰伦理| 在线看三级毛片| 少妇的逼水好多| 国产精品一区二区性色av| 88av欧美| 久久九九热精品免费| 精品一区二区三区av网在线观看| 国产精品一区二区三区四区免费观看 | 男插女下体视频免费在线播放| 日韩精品中文字幕看吧| 国产一区二区亚洲精品在线观看| 久久精品国产鲁丝片午夜精品 | 精品午夜福利视频在线观看一区| 国产爱豆传媒在线观看| a在线观看视频网站| 国产伦在线观看视频一区| 亚洲乱码一区二区免费版| 别揉我奶头~嗯~啊~动态视频| 伦理电影大哥的女人| 成人特级av手机在线观看| 国产69精品久久久久777片| 99久久久亚洲精品蜜臀av| 中亚洲国语对白在线视频| netflix在线观看网站| 亚洲熟妇熟女久久| 亚洲国产精品久久男人天堂| 午夜影院日韩av| 国产精品嫩草影院av在线观看 | 搡老岳熟女国产| 亚洲国产色片| 免费av不卡在线播放| 久久久久久大精品| 精品久久久噜噜| 美女黄网站色视频| 久久精品国产亚洲av天美| 日本五十路高清| 久久午夜亚洲精品久久| 国产欧美日韩精品一区二区| 久久久久久久久大av| 成年女人毛片免费观看观看9| 两个人的视频大全免费| 久久午夜福利片| 欧美潮喷喷水| 他把我摸到了高潮在线观看| 免费看av在线观看网站| 久久精品人妻少妇| 99久久精品国产国产毛片| 午夜日韩欧美国产| 色哟哟·www| 免费在线观看日本一区| 国产免费av片在线观看野外av| 欧美3d第一页| 97碰自拍视频| 男插女下体视频免费在线播放| 一级黄色大片毛片| 国内毛片毛片毛片毛片毛片| 午夜免费男女啪啪视频观看 | 搡老妇女老女人老熟妇| 色综合婷婷激情| 亚洲av免费在线观看| 精品午夜福利视频在线观看一区| 亚洲成人久久性| 久久天躁狠狠躁夜夜2o2o| 尾随美女入室| 欧美精品国产亚洲| 亚洲精品乱码久久久v下载方式| 夜夜看夜夜爽夜夜摸| 成人午夜高清在线视频| 真人做人爱边吃奶动态| 国产欧美日韩精品亚洲av| 99久久成人亚洲精品观看| 久久人人精品亚洲av| 桃色一区二区三区在线观看| 国产高潮美女av| 久久久久九九精品影院| 亚洲av五月六月丁香网| av视频在线观看入口| 久久久久久久久中文| 热99在线观看视频| 美女高潮的动态| 久久精品人妻少妇| 能在线免费观看的黄片| 免费看光身美女| 亚洲七黄色美女视频| 亚洲成av人片在线播放无| 亚洲熟妇中文字幕五十中出| 国产成人a区在线观看| 搡女人真爽免费视频火全软件 | 亚洲成人久久爱视频| 在线观看免费视频日本深夜| 欧美zozozo另类| 伦精品一区二区三区| 亚洲性夜色夜夜综合| 日日夜夜操网爽| 久久精品国产清高在天天线| 久久精品国产99精品国产亚洲性色| 国产精品久久久久久亚洲av鲁大| 久99久视频精品免费| 不卡一级毛片| 久久久久久大精品| 亚洲三级黄色毛片| 国产欧美日韩精品亚洲av| 国产大屁股一区二区在线视频| 美女免费视频网站| 亚洲成a人片在线一区二区| 99久久久亚洲精品蜜臀av| 午夜福利视频1000在线观看| 全区人妻精品视频| 一个人免费在线观看电影| 国产免费男女视频| 一级毛片久久久久久久久女| 亚洲精华国产精华液的使用体验 | 亚洲自拍偷在线| 在线免费观看不下载黄p国产 | aaaaa片日本免费| 身体一侧抽搐| 男人和女人高潮做爰伦理| 亚洲人成网站在线播放欧美日韩| 欧美激情在线99| .国产精品久久| 久久6这里有精品| 欧美人与善性xxx| 国产精品自产拍在线观看55亚洲| 窝窝影院91人妻| 淫秽高清视频在线观看| 欧美高清成人免费视频www| 变态另类丝袜制服| 日日夜夜操网爽| 国产精品爽爽va在线观看网站| 啦啦啦韩国在线观看视频| 美女高潮喷水抽搐中文字幕| av在线老鸭窝| 国产黄a三级三级三级人| 免费在线观看日本一区| 国产成人福利小说| 国产精品野战在线观看| 嫩草影院精品99| 看黄色毛片网站| 成人特级av手机在线观看| 国产伦一二天堂av在线观看| 午夜亚洲福利在线播放| 国产精品久久久久久av不卡| 日韩中文字幕欧美一区二区| 熟女人妻精品中文字幕| 国产真实伦视频高清在线观看 | 丰满人妻一区二区三区视频av| 大又大粗又爽又黄少妇毛片口| 色哟哟·www| 在线观看舔阴道视频| 两个人视频免费观看高清| 国产精品1区2区在线观看.| 俺也久久电影网| 成人毛片a级毛片在线播放| 久久久国产成人免费| 日日啪夜夜撸| 美女高潮的动态| 日本黄大片高清| 国产真实乱freesex| 91久久精品电影网| 亚洲五月天丁香| 国产精品久久久久久久电影| 成年人黄色毛片网站| 夜夜看夜夜爽夜夜摸| 色播亚洲综合网| 老熟妇乱子伦视频在线观看| 中国美女看黄片| 国产精品98久久久久久宅男小说| 日本欧美国产在线视频| 国产精品一区二区性色av| 精品一区二区三区视频在线| 91午夜精品亚洲一区二区三区 | 日本色播在线视频| 欧美丝袜亚洲另类 | 一本精品99久久精品77| 久久精品综合一区二区三区| 亚洲va日本ⅴa欧美va伊人久久| 国模一区二区三区四区视频| 一进一出抽搐gif免费好疼| 国产成人av教育| 亚洲欧美日韩无卡精品| 老师上课跳d突然被开到最大视频| 欧美精品国产亚洲| 久久久久精品国产欧美久久久| 亚洲美女黄片视频| 97热精品久久久久久| 久久久久久久久久成人| 天堂动漫精品| 精品久久久久久成人av| 此物有八面人人有两片| 欧美黑人巨大hd| 成人永久免费在线观看视频| 精品乱码久久久久久99久播| 美女高潮的动态| 在线天堂最新版资源| 免费人成视频x8x8入口观看| 国产在线男女| 99热网站在线观看| 性欧美人与动物交配| 午夜福利高清视频| 18禁裸乳无遮挡免费网站照片| 日韩国内少妇激情av| 国产av不卡久久| 黄色丝袜av网址大全| 亚洲国产精品sss在线观看| 成年人黄色毛片网站| 看十八女毛片水多多多| 欧美性感艳星| a级一级毛片免费在线观看| 波多野结衣高清作品| 啦啦啦韩国在线观看视频| 亚洲无线在线观看| 久久精品国产亚洲av涩爱 | 白带黄色成豆腐渣| www.www免费av| 三级国产精品欧美在线观看| 久久久国产成人精品二区| 日本 av在线| 国产大屁股一区二区在线视频| 国产一区二区亚洲精品在线观看| 伦精品一区二区三区| 免费人成视频x8x8入口观看| 午夜福利高清视频| 麻豆av噜噜一区二区三区| 国产不卡一卡二| 免费在线观看日本一区| 乱人视频在线观看| 欧美绝顶高潮抽搐喷水| 亚洲四区av| 欧美+亚洲+日韩+国产| 国产视频内射| 久久婷婷人人爽人人干人人爱| 国产 一区 欧美 日韩| 99在线人妻在线中文字幕| 国产精品人妻久久久久久| 欧美一区二区亚洲| 网址你懂的国产日韩在线| 亚洲欧美清纯卡通| 22中文网久久字幕| 女人被狂操c到高潮| 在线观看一区二区三区| 内地一区二区视频在线| 99久久精品热视频| 国产精品电影一区二区三区| 俄罗斯特黄特色一大片| 成人高潮视频无遮挡免费网站| 一边摸一边抽搐一进一小说| 免费在线观看影片大全网站| 很黄的视频免费| 黄片wwwwww| 九九爱精品视频在线观看| 深夜精品福利| 免费一级毛片在线播放高清视频| 日韩强制内射视频| 欧美精品啪啪一区二区三区| 99riav亚洲国产免费| 国产高清视频在线播放一区| aaaaa片日本免费| 在线观看一区二区三区| 久久久久久久亚洲中文字幕| 久久精品综合一区二区三区| 在线国产一区二区在线| 国产精品久久久久久av不卡| 我的老师免费观看完整版| 国产视频内射| 精品国产三级普通话版| 国产精品一区二区免费欧美| 国产欧美日韩精品亚洲av| 麻豆一二三区av精品| 国产黄色小视频在线观看| 精品福利观看| 美女cb高潮喷水在线观看| 国产乱人视频| 91麻豆av在线| 国产精品福利在线免费观看| 久久久国产成人精品二区| 久久人人爽人人爽人人片va| 国产伦一二天堂av在线观看| 久久久久久久午夜电影| 亚洲无线观看免费| 国产一区二区三区在线臀色熟女| 中文资源天堂在线| 男人和女人高潮做爰伦理| 免费看a级黄色片| 久久精品国产亚洲av天美| 免费高清视频大片| 久久久久免费精品人妻一区二区| 啪啪无遮挡十八禁网站| 亚洲成人久久性| 国产成人aa在线观看| 国产不卡一卡二| 天堂网av新在线| 国产精品爽爽va在线观看网站| 99久久九九国产精品国产免费| 亚洲av免费高清在线观看| 村上凉子中文字幕在线| 一本精品99久久精品77| 美女高潮喷水抽搐中文字幕| 婷婷亚洲欧美| 亚洲欧美日韩高清专用| 色在线成人网| 在线播放国产精品三级| 成人特级av手机在线观看| 老熟妇乱子伦视频在线观看| 真人一进一出gif抽搐免费| 国内揄拍国产精品人妻在线| 久久久精品大字幕| 欧美日韩乱码在线| av.在线天堂| 97超级碰碰碰精品色视频在线观看| 一区福利在线观看| 99热精品在线国产| 国产高清有码在线观看视频| 亚洲国产精品合色在线| 最后的刺客免费高清国语| 99国产极品粉嫩在线观看| 在线播放国产精品三级| 色综合站精品国产| 欧美激情久久久久久爽电影| 国产精华一区二区三区| 最好的美女福利视频网| 一本精品99久久精品77| 欧美一区二区国产精品久久精品|