• <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
    久久午夜综合久久蜜桃| 别揉我奶头~嗯~啊~动态视频| 午夜福利欧美成人| 中文字幕色久视频| 国产黄a三级三级三级人| 妹子高潮喷水视频| 老鸭窝网址在线观看| 88av欧美| 欧美激情 高清一区二区三区| 亚洲成人精品中文字幕电影| 亚洲精品一区av在线观看| 亚洲成国产人片在线观看| 在线观看免费视频日本深夜| bbb黄色大片| 日本 av在线| 一本综合久久免费| 亚洲七黄色美女视频| 他把我摸到了高潮在线观看| 侵犯人妻中文字幕一二三四区| 18禁美女被吸乳视频| 97碰自拍视频| 淫妇啪啪啪对白视频| 亚洲精品久久国产高清桃花| 成人三级黄色视频| 国产免费男女视频| 亚洲成国产人片在线观看| 亚洲成av片中文字幕在线观看| 国产精品99久久99久久久不卡| 亚洲人成伊人成综合网2020| 91av网站免费观看| 日韩精品中文字幕看吧| 亚洲成国产人片在线观看| 最近最新中文字幕大全免费视频| 好男人电影高清在线观看| 国产成人av激情在线播放| 性欧美人与动物交配| 老鸭窝网址在线观看| 丰满的人妻完整版| 日韩免费av在线播放| 一级a爱视频在线免费观看| 麻豆一二三区av精品| av片东京热男人的天堂| 母亲3免费完整高清在线观看| 天天添夜夜摸| 天堂√8在线中文| 女人高潮潮喷娇喘18禁视频| 日韩大码丰满熟妇| 岛国视频午夜一区免费看| 琪琪午夜伦伦电影理论片6080| 在线观看午夜福利视频| 精品熟女少妇八av免费久了| 一级,二级,三级黄色视频| 欧美丝袜亚洲另类 | 国产激情欧美一区二区| 好看av亚洲va欧美ⅴa在| 欧美日韩黄片免| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲男人天堂网一区| 亚洲av第一区精品v没综合| 免费看美女性在线毛片视频| 日韩欧美免费精品| 色播亚洲综合网| 黄片大片在线免费观看| 51午夜福利影视在线观看| 人人澡人人妻人| 色在线成人网| 欧美日韩黄片免| 欧美成人午夜精品| 亚洲一区二区三区不卡视频| 法律面前人人平等表现在哪些方面| 在线观看免费日韩欧美大片| 精品久久久久久久毛片微露脸| 午夜福利视频1000在线观看 | 午夜亚洲福利在线播放| 国产精品久久久av美女十八| 黄色 视频免费看| 国产成人欧美在线观看| 久久精品国产亚洲av高清一级| 一级a爱片免费观看的视频| 美女大奶头视频| 亚洲国产欧美日韩在线播放| 曰老女人黄片| 欧美 亚洲 国产 日韩一| 看黄色毛片网站| av天堂在线播放| 成人三级做爰电影| 欧美激情久久久久久爽电影 | 国产精品电影一区二区三区| 妹子高潮喷水视频| 丝袜美足系列| 我的亚洲天堂| 夜夜夜夜夜久久久久| 最近最新免费中文字幕在线| 极品教师在线免费播放| 天堂动漫精品| 男女做爰动态图高潮gif福利片 | 国产亚洲精品久久久久久毛片| 精品国产国语对白av| 最近最新中文字幕大全电影3 | 国产欧美日韩精品亚洲av| 久久久久久免费高清国产稀缺| 亚洲av片天天在线观看| 免费观看精品视频网站| 看免费av毛片| 制服人妻中文乱码| 国产av一区在线观看免费| 成人国语在线视频| 桃红色精品国产亚洲av| 99国产精品一区二区蜜桃av| 免费搜索国产男女视频| 亚洲av第一区精品v没综合| 成在线人永久免费视频| 91精品三级在线观看| 久久中文看片网| 日韩欧美国产一区二区入口| 亚洲精品在线观看二区| 女人高潮潮喷娇喘18禁视频| 午夜影院日韩av| av天堂在线播放| av天堂在线播放| 日本免费一区二区三区高清不卡 | 中文字幕人妻丝袜一区二区| 男女做爰动态图高潮gif福利片 | 久久久国产精品麻豆| 国产精品二区激情视频| 欧美老熟妇乱子伦牲交| 欧美黄色片欧美黄色片| 免费在线观看完整版高清| 亚洲男人天堂网一区| 国产成年人精品一区二区| 99精品欧美一区二区三区四区| 欧美激情 高清一区二区三区| 久久欧美精品欧美久久欧美| 成年版毛片免费区| 成人18禁高潮啪啪吃奶动态图| 色av中文字幕| www.精华液| 亚洲男人的天堂狠狠| 色av中文字幕| 午夜成年电影在线免费观看| 性欧美人与动物交配| aaaaa片日本免费| 两人在一起打扑克的视频| 99久久国产精品久久久| 无人区码免费观看不卡| 后天国语完整版免费观看| 国产精品 国内视频| 香蕉丝袜av| 久久久久九九精品影院| 久久久久九九精品影院| 黄片大片在线免费观看| 国产欧美日韩一区二区精品| 男人操女人黄网站| 香蕉国产在线看| 国产精品美女特级片免费视频播放器 | 亚洲一码二码三码区别大吗| 亚洲成人免费电影在线观看| 老汉色∧v一级毛片| 色哟哟哟哟哟哟| 97碰自拍视频| 人人妻,人人澡人人爽秒播| 亚洲全国av大片| 日韩欧美国产一区二区入口| 悠悠久久av| 丝袜在线中文字幕| 国产一级毛片七仙女欲春2 | netflix在线观看网站| 亚洲电影在线观看av| 免费少妇av软件| 男女下面进入的视频免费午夜 | 亚洲精品av麻豆狂野| 亚洲五月天丁香| 中文字幕av电影在线播放| 多毛熟女@视频| 成熟少妇高潮喷水视频| 夜夜躁狠狠躁天天躁| 免费久久久久久久精品成人欧美视频| 亚洲第一欧美日韩一区二区三区| 在线播放国产精品三级| 亚洲精品av麻豆狂野| 免费无遮挡裸体视频| 国产一区二区三区视频了| 国语自产精品视频在线第100页| 熟妇人妻久久中文字幕3abv| 一级毛片女人18水好多| 国产片内射在线| 中出人妻视频一区二区| 大香蕉久久成人网| 免费人成视频x8x8入口观看| 在线永久观看黄色视频| 国产精品国产高清国产av| 啦啦啦韩国在线观看视频| 老司机福利观看| 久久久久久亚洲精品国产蜜桃av| 精品熟女少妇八av免费久了| 成年人黄色毛片网站| cao死你这个sao货| 51午夜福利影视在线观看| 免费女性裸体啪啪无遮挡网站| 亚洲国产看品久久| 国产蜜桃级精品一区二区三区| 国产三级在线视频| 最好的美女福利视频网| 激情视频va一区二区三区| 此物有八面人人有两片| 黄色女人牲交| 99香蕉大伊视频| 日本 欧美在线| 国产精品久久久久久精品电影 | 国产人伦9x9x在线观看| 极品人妻少妇av视频| 大码成人一级视频| 91麻豆av在线| 亚洲av电影不卡..在线观看| 日韩高清综合在线| 变态另类丝袜制服| 婷婷六月久久综合丁香| 色综合亚洲欧美另类图片| 国产精品日韩av在线免费观看 | xxx96com| 在线国产一区二区在线| 午夜激情av网站| 波多野结衣巨乳人妻| 欧洲精品卡2卡3卡4卡5卡区| 亚洲第一欧美日韩一区二区三区| 国内精品久久久久久久电影| 97人妻天天添夜夜摸| 国产成人一区二区三区免费视频网站| 视频区欧美日本亚洲| 黄色女人牲交| 在线观看www视频免费| 国产精品亚洲美女久久久| 精品久久久精品久久久| 波多野结衣巨乳人妻| 免费在线观看日本一区| 精品第一国产精品| 黄色女人牲交| 国产蜜桃级精品一区二区三区| 欧美黑人精品巨大| 久久久久久大精品| a在线观看视频网站| 亚洲av成人av| 成人免费观看视频高清| 久久人妻熟女aⅴ| 国产成人精品久久二区二区91| 国产高清有码在线观看视频 | 色播亚洲综合网| 亚洲五月婷婷丁香| 久久草成人影院| 天堂影院成人在线观看| 97人妻精品一区二区三区麻豆 | 黄网站色视频无遮挡免费观看| 一边摸一边抽搐一进一出视频| 老司机在亚洲福利影院| 99re在线观看精品视频| 欧美激情 高清一区二区三区| 欧美日本中文国产一区发布| 叶爱在线成人免费视频播放| a级毛片在线看网站| 精品一区二区三区四区五区乱码| 日本三级黄在线观看| 在线国产一区二区在线| 麻豆成人av在线观看| 12—13女人毛片做爰片一| 村上凉子中文字幕在线| 午夜福利欧美成人| 国产精品久久久人人做人人爽| 国产91精品成人一区二区三区| 岛国视频午夜一区免费看| 精品久久久久久久毛片微露脸| 国产成人影院久久av| 欧美乱妇无乱码| 亚洲人成电影观看| 中亚洲国语对白在线视频| 日韩av在线大香蕉| av在线天堂中文字幕| 18禁美女被吸乳视频| 嫁个100分男人电影在线观看| 色精品久久人妻99蜜桃| 精品一区二区三区四区五区乱码| 亚洲男人的天堂狠狠| 伦理电影免费视频| 色播在线永久视频| 国产精品国产高清国产av| 欧美一区二区精品小视频在线| 亚洲第一欧美日韩一区二区三区| 非洲黑人性xxxx精品又粗又长| 日韩欧美国产在线观看| 国产成人啪精品午夜网站| 91大片在线观看| 日日摸夜夜添夜夜添小说| 色老头精品视频在线观看| 自线自在国产av| 欧美黄色片欧美黄色片| 国产精品久久久久久人妻精品电影| 十八禁人妻一区二区| 97超级碰碰碰精品色视频在线观看| 女人精品久久久久毛片| 国内精品久久久久久久电影| 亚洲在线自拍视频| 在线观看舔阴道视频| 日本免费一区二区三区高清不卡 | 大香蕉久久成人网| √禁漫天堂资源中文www| 亚洲专区字幕在线| 欧美一区二区精品小视频在线| 宅男免费午夜| 啪啪无遮挡十八禁网站| aaaaa片日本免费| 99精品在免费线老司机午夜| 亚洲五月色婷婷综合| av中文乱码字幕在线| 免费高清视频大片| 久久久久国内视频| 日韩欧美在线二视频| 搞女人的毛片| 97碰自拍视频| 亚洲精品国产精品久久久不卡| 久久人妻福利社区极品人妻图片| tocl精华| 久久午夜综合久久蜜桃| 97人妻天天添夜夜摸| 国产亚洲精品第一综合不卡| 天堂√8在线中文| 亚洲人成电影观看| 午夜成年电影在线免费观看| 精品高清国产在线一区| 身体一侧抽搐| 精品国产美女av久久久久小说| 国产成人精品无人区| 免费不卡黄色视频| 可以在线观看毛片的网站| 亚洲中文字幕一区二区三区有码在线看 | 黄片大片在线免费观看| 日本黄色视频三级网站网址| 手机成人av网站| 免费人成视频x8x8入口观看| 久久久久九九精品影院| 亚洲国产中文字幕在线视频| 无限看片的www在线观看| 亚洲伊人色综图| 亚洲男人天堂网一区| 色精品久久人妻99蜜桃| 美女扒开内裤让男人捅视频| 熟妇人妻久久中文字幕3abv| 国产野战对白在线观看| 成人亚洲精品一区在线观看| 国产精品美女特级片免费视频播放器 | 国产精品免费一区二区三区在线| 国产成人精品久久二区二区91| 亚洲欧美精品综合久久99| 精品福利观看| 精品第一国产精品| 搡老熟女国产l中国老女人| 亚洲欧美精品综合一区二区三区| 一进一出抽搐动态| 九色国产91popny在线| 一边摸一边做爽爽视频免费| 国产xxxxx性猛交| 变态另类丝袜制服| 18美女黄网站色大片免费观看| 亚洲天堂国产精品一区在线| 国产伦一二天堂av在线观看| 韩国精品一区二区三区| 免费av毛片视频| 男女下面插进去视频免费观看| 国产亚洲欧美98| 香蕉久久夜色| 最好的美女福利视频网| 日本免费a在线| 十分钟在线观看高清视频www| 亚洲色图av天堂| 午夜影院日韩av| 制服丝袜大香蕉在线| 国产单亲对白刺激| 好看av亚洲va欧美ⅴa在| 99久久久亚洲精品蜜臀av| av有码第一页| 国产麻豆成人av免费视频| 亚洲午夜理论影院| 午夜福利高清视频| 午夜福利免费观看在线| 别揉我奶头~嗯~啊~动态视频| 成人18禁高潮啪啪吃奶动态图| 亚洲av五月六月丁香网| 丰满的人妻完整版| 午夜福利视频1000在线观看 | av免费在线观看网站| 婷婷六月久久综合丁香| 免费看十八禁软件| 国产激情欧美一区二区| 天天躁夜夜躁狠狠躁躁| videosex国产| 国内精品久久久久久久电影| netflix在线观看网站| 成人三级黄色视频| 国产精品av久久久久免费| 欧美一级毛片孕妇| 国产精品1区2区在线观看.| 国产一区在线观看成人免费| 久久精品人人爽人人爽视色| 国产麻豆69| 国产伦一二天堂av在线观看| 亚洲国产看品久久| 一区二区三区精品91| 高潮久久久久久久久久久不卡| 久久婷婷成人综合色麻豆| 国产精品自产拍在线观看55亚洲| 国产精品 欧美亚洲| 国产日韩一区二区三区精品不卡| 欧美一级a爱片免费观看看 | 国产成人影院久久av| 我的亚洲天堂| 此物有八面人人有两片| 亚洲五月天丁香| 黄色毛片三级朝国网站| 嫩草影视91久久| 久久精品91无色码中文字幕| 可以免费在线观看a视频的电影网站| 天天一区二区日本电影三级 | 老司机福利观看| 欧美成人午夜精品| 国产一区二区在线av高清观看| 免费不卡黄色视频| 手机成人av网站| 国产91精品成人一区二区三区| 后天国语完整版免费观看| 国产成人精品久久二区二区免费| 欧美av亚洲av综合av国产av| 麻豆av在线久日| 啦啦啦韩国在线观看视频| 国产亚洲精品第一综合不卡| 色精品久久人妻99蜜桃| 精品久久久久久久久久免费视频| 91字幕亚洲| 国产高清有码在线观看视频 | 欧美激情久久久久久爽电影 | 亚洲国产中文字幕在线视频| 日韩中文字幕欧美一区二区| av有码第一页| 国产成人啪精品午夜网站| 不卡一级毛片| 丝袜美腿诱惑在线| 91老司机精品| 精品国产一区二区三区四区第35| 国产av又大| 免费无遮挡裸体视频| 一个人免费在线观看的高清视频| 欧美大码av| 美女大奶头视频| 两个人视频免费观看高清| 999久久久国产精品视频| 欧美午夜高清在线| av福利片在线| 久久久久国产精品人妻aⅴ院| 亚洲男人天堂网一区| 亚洲成a人片在线一区二区| 亚洲国产精品合色在线| 亚洲成av片中文字幕在线观看| 日韩一卡2卡3卡4卡2021年| 亚洲色图av天堂| 成人亚洲精品一区在线观看| 国产精品永久免费网站| 亚洲国产毛片av蜜桃av| 操美女的视频在线观看| 午夜日韩欧美国产| 成人国产一区最新在线观看| 午夜精品在线福利| av电影中文网址| 午夜福利免费观看在线| 人人澡人人妻人| 他把我摸到了高潮在线观看| 高清在线国产一区| 丰满的人妻完整版| 国产精品久久电影中文字幕| 精品久久久久久,| 无人区码免费观看不卡| 999精品在线视频| 日韩欧美三级三区| а√天堂www在线а√下载| 成人精品一区二区免费| 亚洲精品久久成人aⅴ小说| 午夜老司机福利片| 午夜激情av网站| 90打野战视频偷拍视频| 黄色丝袜av网址大全| 日本五十路高清| a在线观看视频网站| 又黄又粗又硬又大视频| 国产欧美日韩精品亚洲av| 婷婷丁香在线五月| 少妇粗大呻吟视频| cao死你这个sao货| 男男h啪啪无遮挡| 亚洲第一青青草原| 大码成人一级视频| av中文乱码字幕在线| 如日韩欧美国产精品一区二区三区| 在线十欧美十亚洲十日本专区| 国产一区二区三区综合在线观看| 欧美中文综合在线视频| 午夜影院日韩av| 亚洲人成电影免费在线| 国产精品国产高清国产av| 琪琪午夜伦伦电影理论片6080| 欧美日韩亚洲国产一区二区在线观看| 大型黄色视频在线免费观看| a级毛片在线看网站| 9191精品国产免费久久| 午夜福利,免费看| 精品日产1卡2卡| 999久久久国产精品视频| 日韩大码丰满熟妇| 欧美日韩瑟瑟在线播放| 国产精品自产拍在线观看55亚洲| 黄片播放在线免费| 日韩精品免费视频一区二区三区| 99久久国产精品久久久| 欧美日韩黄片免| 麻豆av在线久日| 欧美 亚洲 国产 日韩一| 亚洲av电影在线进入| 黄片播放在线免费| 国产91精品成人一区二区三区| 欧美一级毛片孕妇| 国产精品久久久久久亚洲av鲁大| 欧美乱妇无乱码| 精品不卡国产一区二区三区| 国产精品精品国产色婷婷| 日日摸夜夜添夜夜添小说| 免费在线观看完整版高清| 午夜久久久在线观看| 一级a爱片免费观看的视频| 国产麻豆69| 一进一出抽搐gif免费好疼| 日韩大码丰满熟妇| 曰老女人黄片| 久久精品国产亚洲av高清一级| 99久久99久久久精品蜜桃| 变态另类丝袜制服| 久久国产精品人妻蜜桃| 黄色女人牲交| 欧美中文综合在线视频| 日本欧美视频一区| 最新在线观看一区二区三区| www.自偷自拍.com| 一级作爱视频免费观看| 不卡一级毛片| 91国产中文字幕| 久久香蕉精品热| 日本 欧美在线| 亚洲国产精品合色在线| 久久人妻福利社区极品人妻图片| 黄色毛片三级朝国网站| 国产精品二区激情视频| 国产人伦9x9x在线观看| 国产主播在线观看一区二区| 亚洲avbb在线观看| 国产精品久久久av美女十八| 十分钟在线观看高清视频www| 韩国精品一区二区三区| 别揉我奶头~嗯~啊~动态视频| 老汉色∧v一级毛片| 国产欧美日韩一区二区精品| 久久久久久免费高清国产稀缺| 大香蕉久久成人网| 久久精品91蜜桃| 给我免费播放毛片高清在线观看| 成人特级黄色片久久久久久久| 又大又爽又粗| 国产精品98久久久久久宅男小说| 啪啪无遮挡十八禁网站| 欧美不卡视频在线免费观看 | 法律面前人人平等表现在哪些方面| 国产一区在线观看成人免费| 亚洲欧美精品综合久久99| 日本欧美视频一区| 午夜福利成人在线免费观看| av欧美777| 亚洲第一电影网av| 日韩欧美一区二区三区在线观看| 国产人伦9x9x在线观看| 精品国产美女av久久久久小说| 久久午夜综合久久蜜桃| 亚洲男人天堂网一区| 怎么达到女性高潮| 久久久水蜜桃国产精品网| 高潮久久久久久久久久久不卡| 欧美日本中文国产一区发布| 久久久久久久久中文| 国产一区二区三区在线臀色熟女| 91九色精品人成在线观看| 久久人妻熟女aⅴ| 美女国产高潮福利片在线看| 人成视频在线观看免费观看| 国产精品国产高清国产av| cao死你这个sao货| 亚洲午夜精品一区,二区,三区| 久久天躁狠狠躁夜夜2o2o| 麻豆成人av在线观看| 大香蕉久久成人网| 激情视频va一区二区三区| 久9热在线精品视频| 亚洲av成人av| 精品国内亚洲2022精品成人| 亚洲熟妇熟女久久| 一区二区三区国产精品乱码| 一个人免费在线观看的高清视频| 国产亚洲精品av在线| 极品教师在线免费播放| 午夜福利影视在线免费观看| 男女下面插进去视频免费观看| 变态另类成人亚洲欧美熟女 | 国产成人av教育| 乱人伦中国视频|