• <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ù)的寬帶測量方法
    纯流量卡能插随身wifi吗| 少妇 在线观看| 中文字幕制服av| 欧美黄色片欧美黄色片| 老鸭窝网址在线观看| 国产精品秋霞免费鲁丝片| 午夜成年电影在线免费观看| 亚洲va日本ⅴa欧美va伊人久久| 天天躁夜夜躁狠狠躁躁| 啦啦啦 在线观看视频| 国产精华一区二区三区| 亚洲欧美日韩高清在线视频| 久久久久久人人人人人| 女性生殖器流出的白浆| 操出白浆在线播放| 亚洲国产欧美网| 建设人人有责人人尽责人人享有的| 亚洲精品久久成人aⅴ小说| 亚洲人成77777在线视频| 日韩有码中文字幕| 日日摸夜夜添夜夜添小说| 精品国产超薄肉色丝袜足j| 国产精品国产av在线观看| 男女下面插进去视频免费观看| 十分钟在线观看高清视频www| 亚洲午夜理论影院| 自线自在国产av| 久久久水蜜桃国产精品网| 青草久久国产| 精品国内亚洲2022精品成人 | 日韩成人在线观看一区二区三区| 人人妻人人澡人人爽人人夜夜| 久久人妻福利社区极品人妻图片| 国产精品亚洲一级av第二区| 国产精品亚洲av一区麻豆| 在线播放国产精品三级| 91大片在线观看| 最近最新中文字幕大全电影3 | 国产xxxxx性猛交| 91麻豆精品激情在线观看国产 | 国产三级黄色录像| 国产午夜精品久久久久久| 精品久久久久久久毛片微露脸| 久久天躁狠狠躁夜夜2o2o| 午夜福利在线观看吧| 国产av又大| 嫁个100分男人电影在线观看| 国产精品国产av在线观看| 国产亚洲欧美在线一区二区| 免费女性裸体啪啪无遮挡网站| 国产精华一区二区三区| 一进一出抽搐动态| 91大片在线观看| 精品高清国产在线一区| 超碰成人久久| 水蜜桃什么品种好| 国产欧美日韩一区二区三区在线| videos熟女内射| 日本黄色日本黄色录像| 欧美在线一区亚洲| 日本撒尿小便嘘嘘汇集6| 亚洲五月天丁香| 国产成人精品久久二区二区91| 久久精品亚洲av国产电影网| 国产熟女午夜一区二区三区| 免费久久久久久久精品成人欧美视频| 无限看片的www在线观看| 好男人电影高清在线观看| 中文字幕人妻丝袜一区二区| 亚洲黑人精品在线| 女警被强在线播放| 国产激情久久老熟女| 丁香六月欧美| 好看av亚洲va欧美ⅴa在| 欧美精品啪啪一区二区三区| 9热在线视频观看99| 亚洲欧美精品综合一区二区三区| 国产视频一区二区在线看| 一边摸一边做爽爽视频免费| 妹子高潮喷水视频| 丝袜在线中文字幕| 中文字幕制服av| 亚洲 国产 在线| 久久久国产欧美日韩av| 免费少妇av软件| 亚洲精品一二三| 久久中文字幕人妻熟女| 俄罗斯特黄特色一大片| 久久人妻福利社区极品人妻图片| 叶爱在线成人免费视频播放| 成人亚洲精品一区在线观看| 视频区欧美日本亚洲| 巨乳人妻的诱惑在线观看| 免费在线观看视频国产中文字幕亚洲| 黄色毛片三级朝国网站| 精品国产国语对白av| 国产成+人综合+亚洲专区| 欧美日韩瑟瑟在线播放| 麻豆乱淫一区二区| av免费在线观看网站| 建设人人有责人人尽责人人享有的| 亚洲中文字幕日韩| 自线自在国产av| 人人妻,人人澡人人爽秒播| 午夜影院日韩av| 国产男靠女视频免费网站| 国产精品久久久av美女十八| 久久亚洲精品不卡| 免费在线观看完整版高清| 免费在线观看完整版高清| 亚洲久久久国产精品| 狠狠婷婷综合久久久久久88av| 午夜亚洲福利在线播放| 啪啪无遮挡十八禁网站| 啪啪无遮挡十八禁网站| 母亲3免费完整高清在线观看| 亚洲性夜色夜夜综合| 人人妻人人澡人人爽人人夜夜| 大香蕉久久网| 久久久久久久午夜电影 | 丰满饥渴人妻一区二区三| 久久婷婷成人综合色麻豆| 伊人久久大香线蕉亚洲五| 国产精品综合久久久久久久免费 | 美女高潮喷水抽搐中文字幕| 天堂√8在线中文| 亚洲熟妇熟女久久| 成年人免费黄色播放视频| avwww免费| av视频免费观看在线观看| 99热只有精品国产| 久久久久精品人妻al黑| 欧美激情高清一区二区三区| 亚洲国产精品一区二区三区在线| 中文欧美无线码| 成人三级做爰电影| 久久久久精品国产欧美久久久| 国产成人av激情在线播放| 丝袜美足系列| 在线观看免费视频日本深夜| 91字幕亚洲| 天天躁狠狠躁夜夜躁狠狠躁| 日韩熟女老妇一区二区性免费视频| 免费人成视频x8x8入口观看| 欧美激情高清一区二区三区| 超色免费av| 99精品久久久久人妻精品| 首页视频小说图片口味搜索| www.精华液| 日日夜夜操网爽| 高清欧美精品videossex| 国产视频一区二区在线看| 大型黄色视频在线免费观看| 国产精品久久久人人做人人爽| 欧美午夜高清在线| 美女扒开内裤让男人捅视频| 国产精品亚洲av一区麻豆| 大片电影免费在线观看免费| 免费观看精品视频网站| 在线观看66精品国产| 激情视频va一区二区三区| av天堂在线播放| 18禁国产床啪视频网站| 亚洲av电影在线进入| 在线看a的网站| 欧美乱妇无乱码| 成年女人毛片免费观看观看9 | 无限看片的www在线观看| 美女国产高潮福利片在线看| 999久久久精品免费观看国产| 人人妻,人人澡人人爽秒播| 精品视频人人做人人爽| 天堂动漫精品| 男女之事视频高清在线观看| 欧美亚洲 丝袜 人妻 在线| 69精品国产乱码久久久| 91麻豆av在线| 精品少妇久久久久久888优播| 老司机亚洲免费影院| 在线观看www视频免费| 国产亚洲精品第一综合不卡| 91麻豆精品激情在线观看国产 | 99国产精品免费福利视频| 日韩欧美免费精品| 欧美日韩福利视频一区二区| 久久久国产一区二区| 最近最新中文字幕大全电影3 | 欧美黄色淫秽网站| 久久国产精品人妻蜜桃| 久久精品亚洲精品国产色婷小说| 黄色 视频免费看| 久久精品aⅴ一区二区三区四区| 黄频高清免费视频| 亚洲五月婷婷丁香| 亚洲第一欧美日韩一区二区三区| 国产日韩一区二区三区精品不卡| 黑人操中国人逼视频| 在线观看免费视频日本深夜| 久久精品国产a三级三级三级| 国产精品国产av在线观看| 国产成人精品在线电影| 国产精品成人在线| 亚洲成a人片在线一区二区| 久久人人爽av亚洲精品天堂| 免费日韩欧美在线观看| 亚洲少妇的诱惑av| www.999成人在线观看| 免费人成视频x8x8入口观看| 99久久国产精品久久久| 国产人伦9x9x在线观看| 精品久久久精品久久久| 国产精品久久久久久精品古装| 免费看a级黄色片| 在线免费观看的www视频| 国产成人精品无人区| 久久久国产成人免费| 国内久久婷婷六月综合欲色啪| 在线国产一区二区在线| 国产精品亚洲一级av第二区| 国产在线观看jvid| 国产精品九九99| 亚洲五月天丁香| 婷婷丁香在线五月| 免费一级毛片在线播放高清视频 | 久久中文字幕一级| 久久性视频一级片| 大片电影免费在线观看免费| 成人精品一区二区免费| 建设人人有责人人尽责人人享有的| 丁香欧美五月| 成人免费观看视频高清| 午夜福利免费观看在线| videosex国产| 亚洲黑人精品在线| 欧美另类亚洲清纯唯美| 在线观看舔阴道视频| 久久久久久免费高清国产稀缺| 久久人人97超碰香蕉20202| 婷婷精品国产亚洲av在线 | 色综合欧美亚洲国产小说| 久久 成人 亚洲| 三上悠亚av全集在线观看| 亚洲中文av在线| 老司机影院毛片| 久久精品国产99精品国产亚洲性色 | 亚洲精品美女久久久久99蜜臀| а√天堂www在线а√下载 | 久久婷婷成人综合色麻豆| 女性生殖器流出的白浆| 久久草成人影院| 国产视频一区二区在线看| 国产午夜精品久久久久久| x7x7x7水蜜桃| 手机成人av网站| 精品国产美女av久久久久小说| 精品熟女少妇八av免费久了| 午夜免费鲁丝| 黄色 视频免费看| 露出奶头的视频| 99久久综合精品五月天人人| 一区在线观看完整版| 中文字幕高清在线视频| 建设人人有责人人尽责人人享有的| 国产不卡av网站在线观看| 黄色片一级片一级黄色片| av电影中文网址| 国产精品香港三级国产av潘金莲| 亚洲五月天丁香| 美女福利国产在线| 日本欧美视频一区| 欧美亚洲 丝袜 人妻 在线| 夜夜躁狠狠躁天天躁| 涩涩av久久男人的天堂| 一级片免费观看大全| 国产精品电影一区二区三区 | 国产精品成人在线| 欧美成人免费av一区二区三区 | 99国产极品粉嫩在线观看| 国产精品国产av在线观看| a在线观看视频网站| xxxhd国产人妻xxx| 亚洲精品在线观看二区| 黄色视频,在线免费观看| 伦理电影免费视频| 亚洲一区中文字幕在线| 精品一区二区三卡| 99精国产麻豆久久婷婷| 十八禁人妻一区二区| 亚洲精品国产一区二区精华液| 王馨瑶露胸无遮挡在线观看| 自线自在国产av| 在线免费观看的www视频| av网站免费在线观看视频| 精品欧美一区二区三区在线| 精品久久久精品久久久| 国产一区二区三区视频了| 亚洲国产看品久久| 日本精品一区二区三区蜜桃| 亚洲熟女精品中文字幕| 成人影院久久| 水蜜桃什么品种好| 另类亚洲欧美激情| 如日韩欧美国产精品一区二区三区| 91精品国产国语对白视频| 亚洲精品中文字幕在线视频| 99国产精品免费福利视频| 国产av又大| 中文字幕高清在线视频| 天天影视国产精品| 日韩人妻精品一区2区三区| 亚洲精品国产色婷婷电影| 黄色 视频免费看| 国内毛片毛片毛片毛片毛片| 色综合欧美亚洲国产小说| 黄网站色视频无遮挡免费观看| 五月开心婷婷网| 免费不卡黄色视频| 曰老女人黄片| 国产精品二区激情视频| 国产麻豆69| 老司机在亚洲福利影院| 男女床上黄色一级片免费看| 免费不卡黄色视频| 国产精品99久久99久久久不卡| 一二三四社区在线视频社区8| 女警被强在线播放| www日本在线高清视频| 日本黄色日本黄色录像| 国产99久久九九免费精品| 亚洲第一青青草原| 黄片播放在线免费| 99精品欧美一区二区三区四区| 黄色丝袜av网址大全| 一区福利在线观看| 欧洲精品卡2卡3卡4卡5卡区| 80岁老熟妇乱子伦牲交| 在线播放国产精品三级| 国产91精品成人一区二区三区| 多毛熟女@视频| 国产又色又爽无遮挡免费看| 欧美日韩亚洲国产一区二区在线观看 | 国产不卡一卡二| 少妇的丰满在线观看| 91成人精品电影| 性少妇av在线| 十八禁高潮呻吟视频| 热99re8久久精品国产| 中文字幕制服av| 国产熟女午夜一区二区三区| 三上悠亚av全集在线观看| 久久香蕉激情| 成人国产一区最新在线观看| 免费不卡黄色视频| 在线视频色国产色| 久久久久视频综合| 一a级毛片在线观看| 自拍欧美九色日韩亚洲蝌蚪91| 午夜影院日韩av| 免费不卡黄色视频| 国产成人系列免费观看| 亚洲综合色网址| 精品久久久久久,| 女人被狂操c到高潮| 亚洲美女黄片视频| 亚洲综合色网址| 国产一区在线观看成人免费| 日韩欧美在线二视频 | 一级毛片女人18水好多| 欧美激情极品国产一区二区三区| 久久 成人 亚洲| 国产精品九九99| 日日摸夜夜添夜夜添小说| 黄色视频不卡| 国产成人精品久久二区二区免费| 成在线人永久免费视频| av电影中文网址| 日韩欧美免费精品| 亚洲一区中文字幕在线| 久久久久久久久久久久大奶| 一边摸一边做爽爽视频免费| 嫁个100分男人电影在线观看| 女同久久另类99精品国产91| 两性夫妻黄色片| 精品国产亚洲在线| 欧美人与性动交α欧美精品济南到| 欧美 日韩 精品 国产| 中文字幕人妻熟女乱码| 久久精品国产a三级三级三级| x7x7x7水蜜桃| 国产在线观看jvid| 夫妻午夜视频| 在线免费观看的www视频| 久久久久久久久免费视频了| 欧美激情极品国产一区二区三区| 精品少妇久久久久久888优播| 亚洲精品一二三| 熟女少妇亚洲综合色aaa.| 亚洲专区中文字幕在线| 亚洲一区二区三区欧美精品| 成人永久免费在线观看视频| 亚洲av成人av| 大码成人一级视频| 精品一区二区三区视频在线观看免费 | 91大片在线观看| 国产精品影院久久| 青草久久国产| 久久精品亚洲精品国产色婷小说| 亚洲片人在线观看| 久9热在线精品视频| 国产精品一区二区精品视频观看| 亚洲精品国产精品久久久不卡| 国产精品亚洲av一区麻豆| 天堂俺去俺来也www色官网| 黄频高清免费视频| 人人澡人人妻人| 国产国语露脸激情在线看| 久久精品亚洲熟妇少妇任你| 国产午夜精品久久久久久| 久久人人爽av亚洲精品天堂| 亚洲中文av在线| 亚洲欧美一区二区三区久久| 亚洲av日韩在线播放| 水蜜桃什么品种好| 午夜精品在线福利| 99久久人妻综合| 亚洲免费av在线视频| 丁香欧美五月| 正在播放国产对白刺激| 捣出白浆h1v1| 搡老熟女国产l中国老女人| 日本a在线网址| 国产精品99久久99久久久不卡| 丝袜人妻中文字幕| 纯流量卡能插随身wifi吗| 性少妇av在线| 亚洲精品久久成人aⅴ小说| av天堂久久9| 国产欧美日韩精品亚洲av| 男男h啪啪无遮挡| 精品国产乱码久久久久久男人| 欧美日韩瑟瑟在线播放| 国产黄色免费在线视频| 久久久久久人人人人人| 在线av久久热| 久热这里只有精品99| 欧美久久黑人一区二区| 午夜福利免费观看在线| 亚洲精品美女久久av网站| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美三级三区| 成人精品一区二区免费| 亚洲av日韩在线播放| 在线观看免费视频日本深夜| 99热国产这里只有精品6| 亚洲伊人色综图| 婷婷丁香在线五月| 99热只有精品国产| 王馨瑶露胸无遮挡在线观看| 99国产精品一区二区蜜桃av | 亚洲一卡2卡3卡4卡5卡精品中文| videos熟女内射| 午夜久久久在线观看| 欧美黑人欧美精品刺激| 精品国产国语对白av| 国产成人精品久久二区二区91| 久久香蕉精品热| 午夜精品久久久久久毛片777| 男人舔女人的私密视频| 人妻久久中文字幕网| 亚洲视频免费观看视频| 一区二区三区国产精品乱码| 最近最新中文字幕大全免费视频| 在线观看www视频免费| 美女 人体艺术 gogo| 国产一区有黄有色的免费视频| 黄色丝袜av网址大全| 在线天堂中文资源库| 一本大道久久a久久精品| 一级,二级,三级黄色视频| 悠悠久久av| 首页视频小说图片口味搜索| 91麻豆精品激情在线观看国产 | 嫩草影视91久久| 国产单亲对白刺激| 18禁国产床啪视频网站| 亚洲av成人一区二区三| 亚洲av欧美aⅴ国产| 久久久精品区二区三区| 99久久精品国产亚洲精品| 亚洲精品久久成人aⅴ小说| 亚洲欧美日韩高清在线视频| 日韩欧美免费精品| 亚洲中文字幕日韩| 极品人妻少妇av视频| 欧美黑人精品巨大| 中文字幕av电影在线播放| 国产亚洲精品久久久久5区| 一二三四在线观看免费中文在| 中文字幕人妻熟女乱码| 国产一区在线观看成人免费| 我的亚洲天堂| 精品久久久久久久毛片微露脸| 成年版毛片免费区| 中文字幕人妻熟女乱码| 国产又色又爽无遮挡免费看| 精品国内亚洲2022精品成人 | 日日夜夜操网爽| 性色av乱码一区二区三区2| 日韩中文字幕欧美一区二区| bbb黄色大片| 国产又爽黄色视频| 怎么达到女性高潮| 黄频高清免费视频| 欧美日韩黄片免| e午夜精品久久久久久久| 另类亚洲欧美激情| 精品卡一卡二卡四卡免费| 亚洲精华国产精华精| 天天躁夜夜躁狠狠躁躁| www.精华液| 国产一区二区激情短视频| a级毛片黄视频| 国产av又大| 在线观看一区二区三区激情| 91老司机精品| 久久人妻熟女aⅴ| 激情视频va一区二区三区| 国产成人免费观看mmmm| 国产精品免费一区二区三区在线 | 91精品三级在线观看| 欧美人与性动交α欧美软件| 制服诱惑二区| 国产一区有黄有色的免费视频| 黑人欧美特级aaaaaa片| 国产伦人伦偷精品视频| 老汉色∧v一级毛片| 中文字幕色久视频| 亚洲av成人av| 亚洲一区中文字幕在线| 日韩有码中文字幕| 久热这里只有精品99| 国产激情欧美一区二区| 亚洲黑人精品在线| 日韩免费av在线播放| www.精华液| 久久久国产精品麻豆| 九色亚洲精品在线播放| 很黄的视频免费| 少妇被粗大的猛进出69影院| 欧美日韩国产mv在线观看视频| 十八禁人妻一区二区| 亚洲美女黄片视频| 人人妻人人添人人爽欧美一区卜| 18禁美女被吸乳视频| 亚洲国产欧美日韩在线播放| 美女扒开内裤让男人捅视频| 五月开心婷婷网| 少妇 在线观看| 精品久久久久久电影网| 亚洲视频免费观看视频| 日韩免费高清中文字幕av| 午夜精品在线福利| 99re在线观看精品视频| 看片在线看免费视频| 夜夜爽天天搞| 日本五十路高清| 久久久久久久久免费视频了| 满18在线观看网站| 亚洲全国av大片| 国产成人啪精品午夜网站| 人人妻人人澡人人爽人人夜夜| 搡老熟女国产l中国老女人| 两性夫妻黄色片| 午夜福利欧美成人| 免费观看精品视频网站| 电影成人av| 国产在线一区二区三区精| 成人18禁高潮啪啪吃奶动态图| 国产精品一区二区在线观看99| 9热在线视频观看99| 精品久久久久久电影网| 深夜精品福利| 一边摸一边抽搐一进一出视频| 美女 人体艺术 gogo| 一进一出好大好爽视频| 欧美人与性动交α欧美软件| 在线观看日韩欧美| 女人高潮潮喷娇喘18禁视频| 校园春色视频在线观看| 欧美日韩中文字幕国产精品一区二区三区 | 亚洲av日韩在线播放| 男人舔女人的私密视频| 制服人妻中文乱码| 天堂√8在线中文| 丝袜美足系列| 久久中文字幕一级| 精品国产乱子伦一区二区三区| 国产精品.久久久| 一本综合久久免费| 天天操日日干夜夜撸| 亚洲专区国产一区二区| 在线观看免费日韩欧美大片| 亚洲全国av大片| 欧美黄色淫秽网站| 色在线成人网| 正在播放国产对白刺激| 亚洲精品国产精品久久久不卡| 中文字幕精品免费在线观看视频| 日韩人妻精品一区2区三区| 女人被躁到高潮嗷嗷叫费观| 久久久久久亚洲精品国产蜜桃av| 成年动漫av网址| 人人妻人人添人人爽欧美一区卜| 亚洲性夜色夜夜综合| 手机成人av网站|