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

    差分干涉測量(D-InSAR)技術(shù)在礦區(qū)地面沉降監(jiān)測中的應(yīng)用

    2010-08-29 01:39:10余景波劉國林曹振坦
    全球定位系統(tǒng) 2010年5期
    關(guān)鍵詞:數(shù)據(jù)處理差分礦區(qū)

    余景波,劉國林,曹振坦

    (山東科技大學(xué)測繪科學(xué)與工程學(xué)院,山東青島266510)

    0 引 言

    采礦引起的地面沉降和塌陷是煤礦開采地區(qū)經(jīng)常發(fā)生的一種破壞性災(zāi)難,這不僅僅是個環(huán)境問題,而且會影響到煤礦附近地區(qū)社會和諧與可持續(xù)發(fā)展[1]。因此,有必要利用先進(jìn)的技術(shù)進(jìn)行監(jiān)測和控制礦區(qū)地表沉降引起的破壞,從而保證礦區(qū)及礦區(qū)附近地區(qū)的安全和穩(wěn)定發(fā)展。以InSAR技術(shù)為基礎(chǔ)進(jìn)行地面沉降監(jiān)測,是國內(nèi)外研究的熱點之一。Gabriel等[2]首次論證了差分干涉測量技術(shù)可用于探測厘米級的地表形變。Wegm等[3]利用DInSAR對德國 Ruhrgebiet地區(qū)地面沉降進(jìn)行研究,獲取了因采礦而導(dǎo)致的地表沉降信息。Ge等[4]將GPS和InSAR技術(shù)結(jié)合在一起進(jìn)行礦區(qū)地表形變測量,準(zhǔn)確獲取了澳大利亞Appin地區(qū)煤礦開采造成的地面沉降信息。吳立新等[5]利用多時相D-InSAR結(jié)合角反射器方法進(jìn)行了工礦區(qū)地表沉降監(jiān)測。

    采用雙軌法和三軌法兩種差分干涉測量方法對6景覆蓋濟(jì)寧某礦區(qū)的ENVISAT ASAR數(shù)據(jù)進(jìn)行差分干涉處理,并對實驗處理結(jié)果進(jìn)行討論和分析。

    1 基本原理

    1.1 InSAR及差分干涉測量(D-InSAR)的基本原理

    InSAR技術(shù)是以合成孔徑雷達(dá)復(fù)數(shù)據(jù)提取的相位信息為信息源獲取地表三維信息和變化信息的一項技術(shù),其基本思想是:利用兩副天線同時成像或一副天線相隔一定時間重復(fù)成像,獲取同一區(qū)域的復(fù)雷達(dá)圖像對,由于兩副天線與地面某一目標(biāo)之間的距離不等,形成干涉圖,干涉圖中的相位值即為兩次成像的相位差測量值,根據(jù)兩次成像的相位差與地面目標(biāo)的三維空間位置之間存在的幾何關(guān)系,利用飛行軌道的參數(shù),即可測定地面目標(biāo)的三維坐標(biāo),它可以用來提供大范圍的高精度數(shù)字高程模型(DEM),并用于探測地表形變[6]。

    差分干涉測量(D-InSAR)技術(shù)是以雷達(dá)干涉測量為出發(fā)點,通過星載或機(jī)載的雷達(dá)天線(單天線或雙天線)主動向地面目標(biāo)發(fā)射電磁波,并接收地面目標(biāo)的后向散射回波,同時記錄各自回波的相位信息。由于在所獲取的相位信息中,含有大氣延遲、軌道誤差、平地效應(yīng)、地形起伏、目標(biāo)兩次成像過程中的微小形變及噪聲信息,因此,D-InSAR的處理過程就是通過外部DEM(“雙軌法”)或其它干涉影像對獲取的相位(“三軌法”或“四軌法”)進(jìn)行差分處理,以便去除無用相位信息而獲取形變相位的一個過程[7]。

    1.2 雙軌法和三軌法數(shù)據(jù)處理流程

    差分干涉測量數(shù)據(jù)處理,經(jīng)過配準(zhǔn)、重采樣、濾波、平地效應(yīng)消除、干涉條紋濾波和相位解纏,相高轉(zhuǎn)換,得到差分干涉圖、增強(qiáng)干涉圖、相干圖、以及地面形變圖,通過進(jìn)一步后處理得到地面沉降分布圖。InSAR數(shù)據(jù)處理的基本流程如圖1所示[8]:

    圖1 InSAR數(shù)據(jù)處理基本流程

    雙軌法是Massonnet等[9]提出的,該方法是利用監(jiān)測區(qū)域地表變化前后兩幅影像生成干涉紋圖,再利用事先獲取的DEM數(shù)據(jù)模擬紋圖,從干涉紋圖中去除地形信息以得到地表變化信息。無需進(jìn)行相位解纏,處理工作量少是該方法的優(yōu)點。但是,對于無DEM 數(shù)據(jù)的地區(qū),無法使用該方法,同時引入DEM數(shù)據(jù)的同時有可能帶入新誤差,這些是該方法的缺點。雙軌法數(shù)據(jù)處理基本流程如圖2所示:

    圖2 雙軌法差分干涉數(shù)據(jù)處理流程

    三軌法是Zebker等[10]人提出的一種數(shù)據(jù)處理方法,該方法利用三幅影像生成兩幅干涉紋圖,一副作為地形干涉對反映地形信息,另一幅作為形變干涉對反映地表形變信息,進(jìn)行平地效應(yīng)消除后,分別進(jìn)行相位解纏,最后利用差分干涉原理計算得到地形信息。三軌法的主要優(yōu)點是無需地面信息、數(shù)據(jù)間的配準(zhǔn)比較容易實現(xiàn);其主要缺點是相位解纏質(zhì)量的好壞將對最終結(jié)果產(chǎn)生影響,但在相干性較高地區(qū)可以獲得較好的結(jié)果。三軌法數(shù)據(jù)處理的基本流程如圖3所示。

    圖3 三軌法差分?jǐn)?shù)據(jù)處理簡單流程圖

    2 實驗研究數(shù)據(jù)選取

    研究的區(qū)域位于濟(jì)寧市的某煤礦,該煤礦位于濟(jì)寧煤田西北部,礦井范圍南北長4~4.5 km,東西寬4~6 km。礦區(qū)內(nèi)地形平坦,地勢東北略高,西南稍低,自然坡度為萬分之七。氣候溫和,四季分明,交通方便。

    利用ENVISAT衛(wèi)星獲取的6景ASAR數(shù)據(jù)進(jìn)行差分干涉測量處理,ASAR數(shù)據(jù)參數(shù)見表1所示。

    表1 所選用的 ASAR數(shù)據(jù)列表

    根據(jù)研究的需要,這6景ASAR數(shù)據(jù)可以組成不同的干涉像對,所需的干涉像對垂直基線和時間間隔參數(shù)情況見表2和表3。

    對表2中的干涉像對進(jìn)行雙軌法差分干涉處理,對表3的干涉像對進(jìn)行三軌法差分干涉處理。

    表2 干涉像對參數(shù)情況列表

    表3 干涉像對參數(shù)情況列表

    雙軌法利用外部DEM數(shù)據(jù)來消除干涉相位圖中的地形因素影響,在實驗中采用分辨率為30 m Aster數(shù)據(jù)作為雙軌法差分干涉處理數(shù)據(jù)的外部DEM。

    3 差分干涉處理和結(jié)果分析

    采用雙軌法處理表1中的干涉像對,把2009年4月3日和2009年6月12日獲取的影像為主圖像,2009年5月8日和2009年7月17日獲取的影像為輔圖像,同時引入分辨率為30 m的Aster數(shù)據(jù)作為外部DEM,進(jìn)行雙軌法差分干涉處理,處理后分別生成差分干涉圖、增強(qiáng)干涉圖、相干圖、地面形變圖,如圖4、5所示。

    在圖4和圖5的地面形變圖中,可以看到盡管所選取的實驗數(shù)據(jù)時間跨度不同,但是差分干涉處理生成地面形變圖中沉降分布范圍基本上大體一致,但是由于干涉選用不同時期的數(shù)據(jù)進(jìn)行差分處理,數(shù)據(jù)相干性存在差異,得出結(jié)果與實際情況會有誤差。這說明雙軌法可以從整體上確定礦區(qū)的沉降分布情況。

    在圖4和圖5的差分干涉圖中,可以看到在圖4的差分干涉圖中沒有明顯的條紋,在圖5差分干涉圖的中央位置有一明顯的條紋,在其他地方都沒有表現(xiàn)出來。從表2中,可以看出生成圖4結(jié)果的干涉像對垂直基線較長,基線長度是相干一個重要的限制因素,垂直基線過長可能導(dǎo)致失相干。造成圖4中差分干涉圖沒明顯干涉條紋的一個重要原因可能是該干涉像對垂直基線過長,造成該干涉像對相干性不強(qiáng)。圖4的增強(qiáng)干涉圖,相干效果不好,表現(xiàn)出雜亂無章的彩色斑點;圖5增強(qiáng)干涉圖中,有一相干性較好的地方,可能這個地方是該時間段,地面沉降最大的所在地,在其它地方基本上沒有干涉出來,只是一些噪聲表現(xiàn)出的雜亂無章的彩色斑點。上面是對雙軌法差分干涉處理數(shù)據(jù)生成干涉圖的分析。由于雙軌法需要引入外部DEM,外部數(shù)據(jù)帶來誤差可能也是導(dǎo)致上面干涉圖現(xiàn)象的其中一個原因。

    將2008年12月19日獲取的影像作為主圖像,分別以2009年 2月 27日和 2009年5月 8日獲取的影像為輔影像,其中,把2008年12月19日和2009年2月27日獲取的影像作為地表形變前獲取的兩幅影像,作為地形干涉對進(jìn)行處理;同樣,把2009年5月8日獲取的影像作為地表形變后獲取的影像,將其與主圖像作為形變干涉對進(jìn)行處理;同樣,以2009年5月8日獲取的影像作為主影像,將2009年5月8日和 2009年6月12日獲取的影像作為地形干涉對進(jìn)行處理;把2009年5月8日和2009年7月17日獲取的影像作為形變干涉對進(jìn)行處理;最后,進(jìn)行三軌法差分干涉處理,獲取差分干涉圖、增強(qiáng)干涉圖相干圖和地面形變圖,如圖6、7所示。

    在圖6和圖7中的地面形變圖中,也可看到沉降分布的相同之處,在形變圖中形成了一些“漏斗區(qū)”。這可以說明,三軌法可以檢測到區(qū)域地表形變發(fā)生的大致位置和范圍,并且從得到的結(jié)果,基本上可以確定地表形變大致位置和范圍。這表明三軌法可以整體上確定礦區(qū)沉降分布情況。

    在圖6和圖7中,圖6中差分干涉圖比圖7中差分干涉圖干涉條紋多些,除了垂直基線長度的原因,季節(jié)可能也是一個原因,即時間去相干。生成圖6結(jié)果的干涉像對所處的時間段正是研究區(qū)域植被稀疏、雨水較少、天氣干燥的時期,這有助于提高獲取的ASAR數(shù)據(jù)差分干涉處理精度,提高干涉效果;而生成圖7結(jié)果的干涉像對,正處于研究區(qū)的夏季,植被茂盛、雨水較多,這可能對干涉效果有影響;再者,生成圖6和圖7結(jié)果,需要先進(jìn)行地形干涉像對差分干涉處理,從表3中,可以看出兩個作為地形干涉像對的時間間隔相差不是很大,前者干涉像對時間處于研究區(qū)冬季,后則處于夏季;并且,由于季節(jié)不同,生成含地形信息干涉圖效果不同,這就造成了最終結(jié)果的差異;再者,作為地形干涉像對的垂直基線長度對DEM的生成具有影響,DEM的質(zhì)量又會影響到差分干涉處理的結(jié)果。在圖6和圖7的增強(qiáng)干涉圖中,干涉條紋明顯程度也不一樣,前者較多、較明顯;后者則較差一些。除了上面的分析外,相位解纏的質(zhì)量也有差別,這也是造成上面現(xiàn)象出現(xiàn)的一個原因。以上是對三軌法數(shù)據(jù)處理生成干涉圖的分析。

    從差分干涉圖到增強(qiáng)干涉圖有時需要除去殘余地形相位,主要有梯度擬合法和基線向量調(diào)整法進(jìn)行殘余地形相位去除兩種方法。在數(shù)據(jù)處理中因干涉圖條紋不明顯,可以看作沒有殘余地形相位,沒有進(jìn)行殘余地形相位去除操作,因而實際存在的殘余地形相位會對增強(qiáng)干涉圖質(zhì)量帶來影響,從而影響最終差分干涉處理的結(jié)果。

    相干圖表示的是兩幅影像的相關(guān)性,在相干性低的地方,可以通過相干圖的對比檢測出礦區(qū)地表沉降的大致位置和范圍。垂直基線長度,對相干圖質(zhì)量影響很大。從圖4~圖7的相干圖中可以看出圖4相干圖質(zhì)量很差,其他三幅相干圖質(zhì)量較好,能夠看出明顯地物的特征。

    對圖4~圖7中地面形變圖進(jìn)一步處理,可以獲取該研究區(qū)域在各時間段的沉降分布,從而可以從宏觀尺度分析該區(qū)域在不同時間段的沉降量和沉降分布。

    從圖8~圖9中,可以看出選取的研究區(qū)域都存在著不同程度的沉降,一些不明顯的沉降可能是開采地下水等原因引起的,也可能是數(shù)據(jù)處理過程不可避免的誤差引起的。在數(shù)據(jù)處理中,選用的差分干涉處理方法不同,最后產(chǎn)生的結(jié)果會有差別。這說明差分干涉測量技術(shù)可以得到礦區(qū)的沉降分布。

    以雙軌法處理的2009年6月12日-2009年7月17日期間獲取的干涉像對和三軌法處理2009年5月8日-2009年6月12日-2009年7月17日期間或取得干涉像對為例,進(jìn)行沉降面積估算。這兩個時間段的沉降面積估算如表4、5所示。

    表4 2009年6月12日-2009年7月17日沉降面積估算(雙軌法)

    表5 2008年12月19日-2009年2月27日-2009年5月8日沉降面積估算(三軌法)

    從表4和表5,可以看出雙軌差分干涉測量和三軌差分干涉測量都可對不同沉降值的沉降面積進(jìn)行估算。這說明差分干涉測量技術(shù)可以對沉降面積進(jìn)行估算。

    4 結(jié) 論

    隨著社會經(jīng)濟(jì)的快速發(fā)展,我國對煤炭資源的開采需求將進(jìn)一步擴(kuò)大,開采深度和范圍也將進(jìn)一步拓展,勢必會造成礦區(qū)地表沉降災(zāi)害的發(fā)生,因此必須對礦區(qū)地表沉陷進(jìn)行監(jiān)測。差分干涉測量技術(shù)具有其自身的優(yōu)勢,越來越受到人們的重視。

    采用雙軌法和三軌法差分干涉處理方法處理了6景覆蓋濟(jì)寧某礦區(qū)ENVISAT ASAR數(shù)據(jù),并對數(shù)據(jù)處理結(jié)果進(jìn)行了分析和討論,從而證明了差分干涉測量(D-InSAR)技術(shù)在礦區(qū)地面沉降監(jiān)測應(yīng)用中的可行性。但是,應(yīng)用差分干涉測量技術(shù)獲取礦區(qū)地表形變信息受多種因素影響,因而需要進(jìn)一步將其測量結(jié)果與水準(zhǔn)測量或GPS測量結(jié)果進(jìn)行比較,以評價其測量精度,這是我們今后進(jìn)一步的研究工作。

    [1]侯建國,楊成生,王亞男,王宏宇.煤礦開采沉降監(jiān)測-基于合成孔徑雷達(dá)差分干涉測量的研究[J].自然災(zāi)害學(xué)報,2009,18(6):71-72.

    [2]Gabriel A K,Goldstein R M,Zebker H A.Mapping Small Elevation Changes Over Large Areas:Differential Radar Interferometry[J].Journal of Geophysiacal Research,1989,94(B7):9183-9191.

    [3]Wegm U,Strozzi T,Werner C,et al.Monitoring of mining-Induced Surface Deformation in the Ruhrgebiet(Germany)with SAR Interferometry[C]∥IGARSS'00,Honolulu,USA,2000.

    [4]Ge L L,Chen H Y,Han S,et al.Intergrated GPS and Interferometric SAR Techniques for Highly Dense Crustal Deformation Monitoring[C]∥The 14th Internal Technology meeting of the Satellite Division of the U.S.Inst.Of Navigation,Salt Lake City,Utah,2001.

    [5]吳立新,高均海,葛大慶,等.基于D-InSAR的煤礦區(qū)開采沉降遙感監(jiān)測技術(shù)分析[J].地理學(xué)與國土研究,2004,20(2):22-25.

    [6]陶秋香,劉國林,馬飛虎,等.軌道差分干涉測量技術(shù)在地面沉降監(jiān)測中的應(yīng)用[J].測量信息與工程,2008,33(4):1-2.

    [7]李晶晶,郭增長.基于D-InSAR技術(shù)的煤礦區(qū)開采沉降監(jiān)測[J].河南理工大學(xué)學(xué)報,2006,25(4):306-309.

    [8]李平湘,楊 杰.雷達(dá)干涉測量原理與應(yīng)用[M].北京:測繪出版社,2006(44):107-108.

    [9]Assonnnet D M.The displacement field of the Lander earth-quake mapped by radar interferometry[J].Nature,1993(364):138-142.

    [10]Zerker H A,Rosen P A,Goldstein R M,etal.On the derivation of coseismic displacement fields using Differential Radar Interferometry of the Landers[J].J Geophys Res,1994,99(B10):19617-19634.

    猜你喜歡
    數(shù)據(jù)處理差分礦區(qū)
    認(rèn)知診斷缺失數(shù)據(jù)處理方法的比較:零替換、多重插補(bǔ)與極大似然估計法*
    ILWT-EEMD數(shù)據(jù)處理的ELM滾動軸承故障診斷
    數(shù)列與差分
    加納Amanforom礦區(qū)Ⅲ號隱伏金礦帶的發(fā)現(xiàn)與評價
    加納Amanforom礦區(qū)Ⅲ號隱伏金礦帶的發(fā)現(xiàn)與評價
    湖北省保康縣堰邊上礦區(qū)發(fā)現(xiàn)超大型磷礦
    廣東省蕉嶺縣作壁坑礦區(qū)探明超大型銣礦
    基于希爾伯特- 黃變換的去噪法在外測數(shù)據(jù)處理中的應(yīng)用
    基于差分隱私的大數(shù)據(jù)隱私保護(hù)
    相對差分單項測距△DOR
    太空探索(2014年1期)2014-07-10 13:41:50
    精品一区二区三区视频在线| 亚洲熟妇中文字幕五十中出| 赤兔流量卡办理| 看黄色毛片网站| 国产精品乱码一区二三区的特点| 网址你懂的国产日韩在线| 欧美日韩瑟瑟在线播放| 亚洲中文字幕一区二区三区有码在线看| 亚洲乱码一区二区免费版| 精品乱码久久久久久99久播| 国产精品一区二区三区四区久久| 精品午夜福利在线看| 亚洲狠狠婷婷综合久久图片| 欧美在线一区亚洲| 国产真实乱freesex| 国产成人av教育| 成熟少妇高潮喷水视频| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 国产成人aa在线观看| 欧美色欧美亚洲另类二区| 国内久久婷婷六月综合欲色啪| 国产伦精品一区二区三区四那| 黄色配什么色好看| 亚洲专区国产一区二区| 亚洲中文字幕日韩| 2021天堂中文幕一二区在线观| 成人三级黄色视频| 天堂动漫精品| 午夜福利免费观看在线| a级毛片a级免费在线| 麻豆国产97在线/欧美| 色综合站精品国产| 男女之事视频高清在线观看| 十八禁人妻一区二区| 亚洲人成伊人成综合网2020| 精品99又大又爽又粗少妇毛片 | 又爽又黄a免费视频| 免费一级毛片在线播放高清视频| 中文字幕av在线有码专区| 99久久九九国产精品国产免费| 熟女电影av网| 久久欧美精品欧美久久欧美| 九九久久精品国产亚洲av麻豆| 国产av不卡久久| 成人无遮挡网站| 国产精品久久电影中文字幕| 黄色日韩在线| 亚洲自拍偷在线| av在线天堂中文字幕| 日韩有码中文字幕| 午夜视频国产福利| 午夜福利18| 国产精品自产拍在线观看55亚洲| 又爽又黄无遮挡网站| 欧美高清成人免费视频www| 大型黄色视频在线免费观看| 国产精品女同一区二区软件 | 国产精品人妻久久久久久| 欧美一区二区精品小视频在线| 久久精品国产清高在天天线| 1024手机看黄色片| 18+在线观看网站| 成人国产一区最新在线观看| 舔av片在线| 在线看三级毛片| 中文资源天堂在线| 日韩中字成人| 日日干狠狠操夜夜爽| 婷婷丁香在线五月| 波野结衣二区三区在线| 宅男免费午夜| 日韩有码中文字幕| 欧美激情国产日韩精品一区| 亚洲片人在线观看| 国产亚洲av嫩草精品影院| 99热这里只有是精品50| 国内少妇人妻偷人精品xxx网站| 久久精品影院6| av女优亚洲男人天堂| 啦啦啦韩国在线观看视频| 国产精品综合久久久久久久免费| 国产一区二区亚洲精品在线观看| 国产视频内射| 老熟妇乱子伦视频在线观看| 深夜精品福利| 美女 人体艺术 gogo| 韩国av一区二区三区四区| av在线观看视频网站免费| 亚洲性夜色夜夜综合| 欧美色欧美亚洲另类二区| 国产伦一二天堂av在线观看| eeuss影院久久| 国产精品久久久久久人妻精品电影| 一卡2卡三卡四卡精品乱码亚洲| 欧美激情久久久久久爽电影| av欧美777| 欧美黑人巨大hd| 别揉我奶头 嗯啊视频| 麻豆成人av在线观看| www.色视频.com| 97热精品久久久久久| 国产又黄又爽又无遮挡在线| 亚洲专区中文字幕在线| 欧美zozozo另类| 欧美高清性xxxxhd video| .国产精品久久| 国产精品久久久久久精品电影| 久久久久国内视频| 99久久精品一区二区三区| 91午夜精品亚洲一区二区三区 | 夜夜夜夜夜久久久久| 欧美+日韩+精品| 国产野战对白在线观看| 欧美xxxx黑人xx丫x性爽| 国产av一区在线观看免费| 好看av亚洲va欧美ⅴa在| 日本在线视频免费播放| 在线观看美女被高潮喷水网站 | 国产中年淑女户外野战色| 精品无人区乱码1区二区| 特大巨黑吊av在线直播| 免费观看精品视频网站| 午夜免费男女啪啪视频观看 | 欧美xxxx黑人xx丫x性爽| 国产亚洲av嫩草精品影院| 简卡轻食公司| 成年人黄色毛片网站| 内地一区二区视频在线| 特级一级黄色大片| 永久网站在线| 亚州av有码| 女人十人毛片免费观看3o分钟| 听说在线观看完整版免费高清| 国产黄色小视频在线观看| 18禁黄网站禁片午夜丰满| 免费人成在线观看视频色| 小说图片视频综合网站| 91av网一区二区| 久久精品国产清高在天天线| 日日干狠狠操夜夜爽| 如何舔出高潮| 免费av观看视频| 亚洲成人久久性| 九九久久精品国产亚洲av麻豆| 两个人视频免费观看高清| 两人在一起打扑克的视频| 亚洲在线观看片| 女同久久另类99精品国产91| 日本与韩国留学比较| 亚洲国产欧洲综合997久久,| av视频在线观看入口| 99热6这里只有精品| 成人一区二区视频在线观看| 12—13女人毛片做爰片一| 欧美黑人欧美精品刺激| a在线观看视频网站| 又黄又爽又刺激的免费视频.| av天堂中文字幕网| 99久久精品热视频| 性色avwww在线观看| 级片在线观看| 黄色一级大片看看| 日韩中字成人| 国产精品爽爽va在线观看网站| 日韩亚洲欧美综合| 亚洲三级黄色毛片| 国产爱豆传媒在线观看| 久久精品影院6| 18禁黄网站禁片午夜丰满| 久久人人精品亚洲av| 宅男免费午夜| 中亚洲国语对白在线视频| 国产高清三级在线| 亚洲中文字幕日韩| 日本精品一区二区三区蜜桃| 1024手机看黄色片| 看片在线看免费视频| 变态另类成人亚洲欧美熟女| 婷婷六月久久综合丁香| 99在线人妻在线中文字幕| 一区二区三区免费毛片| 欧美色欧美亚洲另类二区| 国产精品人妻久久久久久| 国产黄片美女视频| av中文乱码字幕在线| 免费高清视频大片| 午夜精品在线福利| 国内久久婷婷六月综合欲色啪| 精品福利观看| 亚洲国产精品sss在线观看| 88av欧美| 色精品久久人妻99蜜桃| 91麻豆av在线| 天天一区二区日本电影三级| 精品人妻视频免费看| 我要看日韩黄色一级片| 免费在线观看成人毛片| 天堂av国产一区二区熟女人妻| 永久网站在线| 亚洲精品日韩av片在线观看| 观看美女的网站| 麻豆av噜噜一区二区三区| 99久久99久久久精品蜜桃| 亚洲精品亚洲一区二区| 我要看日韩黄色一级片| aaaaa片日本免费| 老司机深夜福利视频在线观看| 午夜a级毛片| 夜夜爽天天搞| 久久久精品大字幕| 美女大奶头视频| 久久性视频一级片| 黄色一级大片看看| 黄色视频,在线免费观看| 国产在视频线在精品| 又粗又爽又猛毛片免费看| 性插视频无遮挡在线免费观看| 国产探花在线观看一区二区| 久久久成人免费电影| 国产中年淑女户外野战色| 三级男女做爰猛烈吃奶摸视频| a级毛片a级免费在线| 黄色丝袜av网址大全| 婷婷精品国产亚洲av在线| 在线观看午夜福利视频| 久久久久免费精品人妻一区二区| 哪里可以看免费的av片| 人妻久久中文字幕网| 亚洲成人久久性| 日韩av在线大香蕉| netflix在线观看网站| 久久久久久久久大av| 日本撒尿小便嘘嘘汇集6| 99视频精品全部免费 在线| 国内毛片毛片毛片毛片毛片| 美女被艹到高潮喷水动态| 一a级毛片在线观看| 麻豆国产av国片精品| 黄色一级大片看看| 一级毛片久久久久久久久女| 免费在线观看日本一区| 国产精品一及| 精品久久久久久成人av| 亚洲真实伦在线观看| 18禁裸乳无遮挡免费网站照片| 中文资源天堂在线| 国产一区二区亚洲精品在线观看| 亚洲人成网站高清观看| av欧美777| 免费在线观看亚洲国产| 99久久久亚洲精品蜜臀av| 丝袜美腿在线中文| 少妇的逼水好多| 国产精品久久久久久精品电影| 精品无人区乱码1区二区| 成人特级黄色片久久久久久久| 免费观看人在逋| 免费观看精品视频网站| 在线免费观看的www视频| 欧美中文日本在线观看视频| 波多野结衣高清无吗| 免费观看的影片在线观看| 日本熟妇午夜| 九九热线精品视视频播放| 久久国产精品人妻蜜桃| 少妇被粗大猛烈的视频| aaaaa片日本免费| 亚洲一区二区三区色噜噜| 欧美激情在线99| 亚洲av中文字字幕乱码综合| АⅤ资源中文在线天堂| 日韩欧美国产在线观看| 久久国产乱子免费精品| 亚洲成a人片在线一区二区| 午夜福利在线观看免费完整高清在 | 午夜福利18| 麻豆成人午夜福利视频| av天堂中文字幕网| 一个人观看的视频www高清免费观看| 国产精品电影一区二区三区| 欧美黑人欧美精品刺激| 老女人水多毛片| 一a级毛片在线观看| 欧美成人a在线观看| 中文字幕高清在线视频| АⅤ资源中文在线天堂| 国产成人a区在线观看| av欧美777| 老司机午夜十八禁免费视频| 中文资源天堂在线| av欧美777| 国产精华一区二区三区| 可以在线观看毛片的网站| 久久久久久九九精品二区国产| 国产一区二区亚洲精品在线观看| 男人和女人高潮做爰伦理| 国产成人a区在线观看| 99国产极品粉嫩在线观看| 亚洲五月天丁香| 国内久久婷婷六月综合欲色啪| 午夜两性在线视频| 久久婷婷人人爽人人干人人爱| 亚洲最大成人手机在线| 精品一区二区三区av网在线观看| 国产乱人视频| 国产精品日韩av在线免费观看| 免费在线观看亚洲国产| 91狼人影院| 最新中文字幕久久久久| 欧美成人性av电影在线观看| 中文字幕免费在线视频6| 日韩中文字幕欧美一区二区| 国产91精品成人一区二区三区| 午夜免费激情av| 在线国产一区二区在线| 男人舔女人下体高潮全视频| 日韩高清综合在线| 人妻丰满熟妇av一区二区三区| 人人妻,人人澡人人爽秒播| 亚洲aⅴ乱码一区二区在线播放| 乱码一卡2卡4卡精品| 久久久久久九九精品二区国产| 亚洲av一区综合| 直男gayav资源| 此物有八面人人有两片| 国产一区二区三区在线臀色熟女| 精品99又大又爽又粗少妇毛片 | 九九在线视频观看精品| 成人特级av手机在线观看| 特大巨黑吊av在线直播| ponron亚洲| 欧美成狂野欧美在线观看| 亚洲精品一卡2卡三卡4卡5卡| 亚洲av一区综合| 国产精品人妻久久久久久| 国产三级黄色录像| 天堂网av新在线| 国产91精品成人一区二区三区| 国产久久久一区二区三区| 自拍偷自拍亚洲精品老妇| 日韩 亚洲 欧美在线| 1000部很黄的大片| 亚洲人成电影免费在线| 每晚都被弄得嗷嗷叫到高潮| 精品熟女少妇八av免费久了| 亚洲成人久久性| 男人和女人高潮做爰伦理| 欧美区成人在线视频| 久久亚洲真实| 又紧又爽又黄一区二区| 99riav亚洲国产免费| 国产精品98久久久久久宅男小说| 精品久久久久久久末码| 小说图片视频综合网站| 欧美黑人巨大hd| 国产毛片a区久久久久| 观看美女的网站| 成人永久免费在线观看视频| 18禁在线播放成人免费| 成人永久免费在线观看视频| 最近中文字幕高清免费大全6 | 在线观看免费视频日本深夜| 91久久精品国产一区二区成人| 欧美乱色亚洲激情| 国产伦人伦偷精品视频| 99国产精品一区二区三区| 亚洲av一区综合| 亚洲无线观看免费| 在线观看美女被高潮喷水网站 | 国产精品日韩av在线免费观看| 99热精品在线国产| 精品久久久久久久久av| a级一级毛片免费在线观看| 久久精品影院6| 人妻丰满熟妇av一区二区三区| 国产高清视频在线观看网站| 啦啦啦韩国在线观看视频| 欧美精品国产亚洲| 国产色爽女视频免费观看| 色吧在线观看| 欧美一级a爱片免费观看看| 男女床上黄色一级片免费看| 久久久久久久亚洲中文字幕 | 在线播放无遮挡| 嫩草影院精品99| 在线观看舔阴道视频| 国产精品精品国产色婷婷| 久久久久免费精品人妻一区二区| 免费观看人在逋| 老司机深夜福利视频在线观看| 国产伦人伦偷精品视频| 成年版毛片免费区| 成人无遮挡网站| 久久香蕉精品热| 嫁个100分男人电影在线观看| 九九在线视频观看精品| 亚洲欧美日韩高清专用| 好男人电影高清在线观看| 日韩 亚洲 欧美在线| 嫩草影视91久久| x7x7x7水蜜桃| 无遮挡黄片免费观看| 蜜桃亚洲精品一区二区三区| 97人妻精品一区二区三区麻豆| 蜜桃亚洲精品一区二区三区| 久久亚洲精品不卡| 十八禁国产超污无遮挡网站| 变态另类丝袜制服| 午夜免费男女啪啪视频观看 | 少妇的逼好多水| 国产精品野战在线观看| 久久久精品欧美日韩精品| 亚洲国产日韩欧美精品在线观看| 欧美性感艳星| 久久久久免费精品人妻一区二区| 亚洲av熟女| 国产精品国产高清国产av| 男女视频在线观看网站免费| 国产老妇女一区| 久久性视频一级片| 91久久精品电影网| 亚洲欧美激情综合另类| 亚洲国产精品999在线| 亚洲avbb在线观看| 国产精品免费一区二区三区在线| 淫妇啪啪啪对白视频| 久久精品影院6| 精品久久久久久久人妻蜜臀av| 乱人视频在线观看| 欧美激情久久久久久爽电影| 高清毛片免费观看视频网站| 悠悠久久av| 国产麻豆成人av免费视频| 波多野结衣高清无吗| 国产成人欧美在线观看| avwww免费| 乱码一卡2卡4卡精品| 国产高清视频在线观看网站| 99久久99久久久精品蜜桃| 天美传媒精品一区二区| av福利片在线观看| 一本一本综合久久| 久久欧美精品欧美久久欧美| 人妻久久中文字幕网| 精品一区二区免费观看| 99精品在免费线老司机午夜| 一级毛片久久久久久久久女| 国语自产精品视频在线第100页| 日韩高清综合在线| 久久热精品热| 国产真实乱freesex| 成人av一区二区三区在线看| 99久久精品国产亚洲精品| 深夜a级毛片| 欧美另类亚洲清纯唯美| 亚洲av美国av| 波多野结衣巨乳人妻| 少妇裸体淫交视频免费看高清| 婷婷亚洲欧美| 男人和女人高潮做爰伦理| 亚洲 国产 在线| 男人狂女人下面高潮的视频| 久久久成人免费电影| 怎么达到女性高潮| 舔av片在线| 激情在线观看视频在线高清| 亚洲五月婷婷丁香| 精品人妻偷拍中文字幕| 久久久国产成人精品二区| 男女做爰动态图高潮gif福利片| 中文在线观看免费www的网站| 少妇的逼好多水| 搡老熟女国产l中国老女人| 免费搜索国产男女视频| .国产精品久久| 99在线视频只有这里精品首页| 国产三级中文精品| 色av中文字幕| 欧美黄色片欧美黄色片| 国产亚洲精品综合一区在线观看| 91久久精品国产一区二区成人| 精品一区二区三区视频在线观看免费| 国产成人av教育| 欧美在线黄色| 精品一区二区三区av网在线观看| 极品教师在线免费播放| 国语自产精品视频在线第100页| 久久久色成人| 午夜福利18| 99久久成人亚洲精品观看| 亚洲精品粉嫩美女一区| 欧美激情国产日韩精品一区| 亚洲熟妇中文字幕五十中出| 91九色精品人成在线观看| 老鸭窝网址在线观看| 免费在线观看影片大全网站| 国产精品伦人一区二区| 免费电影在线观看免费观看| 在线国产一区二区在线| 亚洲黑人精品在线| 国产一区二区三区视频了| 精品无人区乱码1区二区| 一a级毛片在线观看| 亚洲在线自拍视频| 日韩欧美免费精品| 欧美区成人在线视频| 午夜激情欧美在线| av欧美777| 99热精品在线国产| 搡老妇女老女人老熟妇| 亚洲第一区二区三区不卡| 99热6这里只有精品| 亚洲片人在线观看| 人妻丰满熟妇av一区二区三区| 一级毛片久久久久久久久女| 中文亚洲av片在线观看爽| 在线免费观看不下载黄p国产 | 亚洲国产精品成人综合色| 夜夜看夜夜爽夜夜摸| xxxwww97欧美| 香蕉av资源在线| 免费av观看视频| 性插视频无遮挡在线免费观看| 久久久久久九九精品二区国产| 亚洲第一欧美日韩一区二区三区| 五月玫瑰六月丁香| 特级一级黄色大片| 午夜精品在线福利| 搡老岳熟女国产| 国产成人a区在线观看| 一级毛片久久久久久久久女| 九九在线视频观看精品| 99热这里只有是精品在线观看 | 人人妻人人看人人澡| 久久久久久久精品吃奶| 一区二区三区高清视频在线| 亚洲欧美精品综合久久99| 一二三四社区在线视频社区8| 久久久久免费精品人妻一区二区| 桃红色精品国产亚洲av| 日本与韩国留学比较| 九色国产91popny在线| 舔av片在线| 国产亚洲精品久久久久久毛片| 免费在线观看亚洲国产| 我的女老师完整版在线观看| 亚洲一区二区三区不卡视频| 一本久久中文字幕| 日韩欧美精品v在线| 久久久久久久亚洲中文字幕 | 亚洲国产高清在线一区二区三| 亚洲欧美激情综合另类| 亚洲av.av天堂| 亚洲成av人片在线播放无| 午夜福利在线观看吧| 波多野结衣高清无吗| 天堂影院成人在线观看| 噜噜噜噜噜久久久久久91| av天堂中文字幕网| 老司机午夜十八禁免费视频| 精品乱码久久久久久99久播| 亚洲一区二区三区不卡视频| 男人和女人高潮做爰伦理| 3wmmmm亚洲av在线观看| 可以在线观看毛片的网站| 老鸭窝网址在线观看| 嫩草影视91久久| 一本精品99久久精品77| 俄罗斯特黄特色一大片| 亚洲第一区二区三区不卡| 欧美成人免费av一区二区三区| 九色成人免费人妻av| 国产黄片美女视频| 很黄的视频免费| 免费电影在线观看免费观看| 两性午夜刺激爽爽歪歪视频在线观看| 亚洲最大成人手机在线| 色噜噜av男人的天堂激情| 麻豆一二三区av精品| 午夜老司机福利剧场| 村上凉子中文字幕在线| 成人美女网站在线观看视频| 欧美激情在线99| 亚洲 国产 在线| 国产大屁股一区二区在线视频| 一级毛片久久久久久久久女| 1000部很黄的大片| 别揉我奶头 嗯啊视频| 亚洲性夜色夜夜综合| 永久网站在线| 欧美极品一区二区三区四区| 欧美性猛交黑人性爽| 熟女电影av网| 女人被狂操c到高潮| 精品久久久久久久久av| 亚洲专区国产一区二区| 久久久国产成人精品二区| 99在线人妻在线中文字幕| 国产精品98久久久久久宅男小说| 色视频www国产| 在线免费观看的www视频| 亚洲国产精品合色在线| 亚洲成人久久爱视频| 久久久久久久久大av| 免费看光身美女| 18禁黄网站禁片午夜丰满| 久久中文看片网| 日本熟妇午夜| 美女 人体艺术 gogo| 88av欧美| 老熟妇乱子伦视频在线观看| 97超级碰碰碰精品色视频在线观看| 国产成年人精品一区二区| 欧美精品啪啪一区二区三区| 亚洲人成伊人成综合网2020|