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

    基于PS-InSAR的曹妃甸填海區(qū)域形變監(jiān)測(cè)

    2020-09-03 03:56:18陳星彤
    礦山測(cè)量 2020年4期
    關(guān)鍵詞:曹妃甸海區(qū)差分

    李 星,陳星彤

    (華北理工大學(xué),河北 唐山 063000)

    從2008年起,曹妃甸填海區(qū)為了響應(yīng)‘一港雙城’的政策,開始了大規(guī)模的加速發(fā)展。尤其是華北理工大學(xué)等一系列高校的搬入組成了曹妃甸大學(xué)城,京唐首鋼等各種鋼廠、鐵廠的搬入所引起的大規(guī)模的人流擾動(dòng)、基礎(chǔ)設(shè)施的建設(shè),勢(shì)必會(huì)造成地表的擾動(dòng)。為了更好地建設(shè)曹妃甸填海區(qū)以及保護(hù)人民財(cái)產(chǎn)安全,對(duì)曹妃甸填海區(qū)進(jìn)行地表形變監(jiān)測(cè)是極其有必要的。

    以往的水準(zhǔn)測(cè)量等監(jiān)測(cè)地表形變受環(huán)境天氣影響,而合成孔徑雷達(dá)(SAR)獲取數(shù)據(jù)有全天時(shí)、全天候的特點(diǎn)。因此學(xué)者們提出了用合成孔徑雷達(dá)干涉測(cè)量(InSAR)技術(shù)、合成孔徑雷達(dá)差分干涉測(cè)量技術(shù)(D-InSAR)監(jiān)測(cè)地表形變。但該技術(shù)受時(shí)空失相干和大氣延遲問題的影響較為嚴(yán)重,學(xué)者們又提出了永久散射體差分干涉技術(shù)(PS-InSAR)。Kiseleva對(duì)Caucasus的黑海海岸進(jìn)行地表沉降監(jiān)測(cè),得到了很好的效果[1]。Gokhan對(duì)Istanbul地區(qū)采用PS-InSAR技術(shù)進(jìn)行地面運(yùn)動(dòng)監(jiān)測(cè),并對(duì)沉降原因進(jìn)行分析[2]。裴媛媛等、劉鵬等、Peng H等采用永久散射體時(shí)序分析(PS-InSAR)技術(shù)分別對(duì)上海市臨港新城主城區(qū)、黃河三角洲東營地區(qū)、杭州灣南部沿海地區(qū)這三個(gè)地區(qū)進(jìn)行地表形變監(jiān)測(cè),反演了各個(gè)地表形變過程,分析其沉降原因[3-5]。

    本文詳細(xì)介紹了 PS-InSAR技術(shù)用于形變監(jiān)測(cè)的技術(shù)原理,在此基礎(chǔ)上對(duì)20幅覆蓋唐山市曹妃甸區(qū)的Sentinel_1A影像進(jìn)行了沉降監(jiān)測(cè)實(shí)驗(yàn),通過取得的沉降規(guī)律,加以分析其沉降原因。

    1 PS-InSAR的原理與關(guān)鍵技術(shù)

    1.1 PS-InSAR的原理及技術(shù)流程

    PS-InSAR技術(shù)是從地區(qū)一致時(shí)間不同的多幅SAR影像中提取出散射特性穩(wěn)定的永久散射體(persistent scatterers,PS)點(diǎn),如人工建筑體、裸露的巖石等硬物體,用提取出的PS點(diǎn)構(gòu)建網(wǎng)絡(luò),并提取其相位信息,從而建立關(guān)于相位的函數(shù)模型,采用合理的算法,精確地反演地表形變。

    在K幅SAR影像中,首先對(duì)主影像進(jìn)行選取,將其與剩余影像配準(zhǔn)及干涉處理,生成K-1幅干涉條紋圖;去除用外部參考DEM數(shù)據(jù)模擬的高程相位,生成去平后的干涉圖;在差分干涉處理后得到的平均強(qiáng)度圖像上提取PS點(diǎn),對(duì)PS點(diǎn)進(jìn)行線性形變相位和DEM誤差的估計(jì),從差分干涉相位中將其去除;通過加入在時(shí)間域和空間域的高通、低通濾波去除獲得的殘余相位信息,包括非線性形變相位、大氣相位、噪聲相位;最后解纏形變相位,對(duì)解纏結(jié)果時(shí)序分析,得到可靠的地表形變信息。

    1.2 關(guān)鍵技術(shù)

    (1)公共主影像的選取

    影響干涉相位的精度以及準(zhǔn)確性的因素有很多,時(shí)間基線、空間基線和Doppler質(zhì)心頻率基線是最直接的,而這三者會(huì)隨著主影像的不同而變化。因此,主影像的選取是PS-InSAR技術(shù)中的一個(gè)關(guān)鍵問題。多個(gè)影響因素綜合考慮,讓其成為最佳組合,即綜合相關(guān)系數(shù)γm最大:

    (1)

    (2)PS點(diǎn)的選取

    選取永久散射體(PS)點(diǎn)也是PS-InSAR技術(shù)中的關(guān)鍵步驟。該技術(shù)通過分析處理這些高相干點(diǎn)目標(biāo)上的干涉相位,從而得到時(shí)間序列上的某地區(qū)形變信息。在處理過程中,建模的條件是影像上要有足夠數(shù)量和密度的PS點(diǎn)。選取PS點(diǎn)的方法主要有相位離差閾值法、相干系數(shù)閾值法和振幅離差閾值法的單種閾值方法以及兩者結(jié)合或三者結(jié)合的復(fù)合閾值方法。PS點(diǎn)選取的方法中最為經(jīng)典的是振幅離差閾值法。該方法的基本原理是,信噪比高的點(diǎn)像元的振幅離差指數(shù)可以和其相位的噪聲水平等價(jià)衡量。對(duì)于標(biāo)準(zhǔn)差為δ的高斯噪聲的SAR影像而言,其振幅信息和 Rice 分布相似,若點(diǎn)具有高斯信噪比的特征,則該點(diǎn)的分布與高斯分布相近,即振幅離差指數(shù)DA:

    (2)

    式中,δA為平均強(qiáng)度圖像的振幅標(biāo)準(zhǔn)差,mA為平均強(qiáng)度圖像的振幅平均值[7]。

    (3)干涉相位模型及形變相位反演

    PS-InSAR技術(shù)中最主要、最關(guān)鍵的一步就是相位模型的建立以及形變相位的反演。每個(gè)提取的PS點(diǎn)的解纏干涉相位φunw:

    φunw=φtopo+φdef+φatm+φnoise

    (3)

    式中,φtopo為DEM高程誤差引起的地形相位;φdef為沿LOS方向的地表形變相位;φatm為大氣影響相位;φnoise為噪聲相位。

    在時(shí)間域與空間域上DEM高程誤差引起的地形相位φtopo、大氣影響相位φatm、噪聲相位φnoise有不同的表現(xiàn)形式,可通過采用高通濾波和低通濾波的方式來去除,最后,形變相位φdef:

    (4)

    式中,φlin是線性的形變相位;φnon-lin是非線性的形變相位,V是線性的形變速率;T是時(shí)間間隔;Dnon-lin是非線性形變。

    由此LOS方向線性形變速率得以獲取[8]。

    2 曹妃甸填海區(qū)形變監(jiān)測(cè)與分析

    2.1 研究區(qū)概況及數(shù)據(jù)選取

    唐山市曹妃甸區(qū)位于唐山南部、環(huán)渤海中心地帶,地理坐標(biāo)為北緯39°07′~39°27′,東經(jīng)118°12′~118°43′。曹妃甸區(qū)的地質(zhì)構(gòu)造為填海造陸,土質(zhì)松軟,近幾年曹妃甸區(qū)進(jìn)行了大規(guī)模建設(shè),引進(jìn)工業(yè)、商業(yè)和學(xué)校。而地表形變會(huì)對(duì)各種基礎(chǔ)設(shè)施的安全造成威脅,這是曹妃甸填海區(qū)地表變形監(jiān)測(cè)的原因。

    本文選取20景2018年8月~2019年4月,覆蓋曹妃甸地區(qū)的Sentinel_1A影像,其分辨率為5 m×20 m。對(duì)曹妃甸填海區(qū)進(jìn)行裁剪,裁剪區(qū)域如圖1黑色線框所示,外部的DEM數(shù)據(jù)采用的是分辨率為30 m的STRM數(shù)據(jù)。

    圖1 研究區(qū)域

    2.2 形變信息提取與分析

    通過公共主影像的選取方法,最終確定2018年11月30日的圖像為公共主影像,并與其他輔影像進(jìn)行配準(zhǔn)、干涉,而干涉數(shù)據(jù)的時(shí)間基線與空間基線的連接圖如圖2所示。從圖中可以看出,空間基線的范圍是在-76.246 1 ~ 96.068 9 m,時(shí)間基線范圍是在-120 ~ 131 d。最后干涉處理完成后得到了19對(duì)干涉條紋圖,如圖3所示。圖中大部分干涉條紋圖的條紋不是很明顯,對(duì)后面PS點(diǎn)的提取有極大的影響。去除干涉圖中利用外部DEM數(shù)據(jù)模擬的高程相位后,得到的去平后的差分干涉圖效果更好,如圖4所示。

    圖2 數(shù)據(jù)的連接圖

    在生成差分干涉圖之后,開始選取PS點(diǎn),由于其地形條件,將振幅離差指數(shù)閾值DA設(shè)置為0.8,PS點(diǎn)的密度為25個(gè)/km2,總共選取了1 107 106個(gè)PS點(diǎn),結(jié)果如圖5所示。PS點(diǎn)大都落在人工建筑,如房屋、道路等硬面地物上,而不在草地、養(yǎng)殖池等散射特性較弱的地方。

    圖3 干涉條紋圖

    圖4 去平后的差分干涉圖

    圖5 永久散射體(PS)分布

    PS點(diǎn)選取后,對(duì)其進(jìn)行線性形變相位和DEM誤差的估計(jì),從差分干涉相位中去除得到殘余相位。由于大氣引起的相位延遲和噪聲屬于高頻信號(hào),而非線性形變相位屬于低頻信號(hào),則可以分別采用時(shí)間域的高通濾波和空間域的低通濾波去除。2019年2月22日的影像能更直觀的表達(dá)濾波前后的效果,如圖6所示,可以看出濾波后的影像相較于濾波前降低了相干噪聲,有效地去除了一些由大氣引起的噪聲,效果明顯。

    圖6 差分干涉圖(a為濾波前,b為濾波后)

    由此,對(duì)濾波后的影像進(jìn)行線性形變相位解纏,得到真實(shí)的相位信息以及各PS點(diǎn)在LOS方向的形變速率信息,如圖7所示。

    圖7 曹妃甸填海區(qū)的年平均沉降速率圖

    由圖7可以看出,整體的平均沉降速率約為-26.56~22.45 mm/a,可見從2018年8月起該地區(qū)仍然存在著較為嚴(yán)重的沉降。紅圈范圍內(nèi)的A區(qū)和B區(qū)是沉降嚴(yán)重的區(qū)域,C區(qū)為抬升區(qū)。A區(qū)位于曹妃甸區(qū)填海區(qū),緊靠海堤,該區(qū)域的沉降速率達(dá)到了-23 mm/a以上,沉降尤為顯著。B區(qū)位于曹妃甸區(qū)的柳贊鎮(zhèn)附近,沉降速率最大為-26.532 5 mm/a,是整個(gè)研究區(qū)域沉降最為嚴(yán)重的地區(qū)。C區(qū)為抬升區(qū)最明顯的地區(qū),最大抬升速率為22.440 8 mm/a。

    為了更好的探討其沉降原因,首先對(duì)大部分沉降文獻(xiàn)進(jìn)行了研究,發(fā)現(xiàn)除了地震等大型自然災(zāi)害的直接因素外,地表沉降的主要原因是地下水的開采和地上人工建筑物建造這兩個(gè)間接因素。其次,查找該區(qū)域的水文地質(zhì)資料,該區(qū)域的含水層形成原因是河流沖積和海積,屬第四系松散巖類,其空隙水一般分為Ⅰ、Ⅱ、Ⅲ組,第Ⅰ、第Ⅱ含水組屬潛水及微承壓水,一般不作為開采層。第Ⅲ含水組為深層承壓水,底界深度一般為160~500 m,水位埋深約為l5 m。20世紀(jì)80年代由于超采深層地下水形成以南堡開發(fā)區(qū)和唐??h城為中心的兩個(gè)沉降漏斗,受這兩個(gè)沉降漏斗的影響,該區(qū)地下水位呈下降趨勢(shì)。這也是該區(qū)地面整體下降的主要原因。就A區(qū)而言,其位于曹妃甸工業(yè)區(qū),我們對(duì)該區(qū)域進(jìn)行了實(shí)地考察,發(fā)現(xiàn)其附近有原油碼頭和鋼廠(首鋼京唐等)。2018年首鋼京唐共生產(chǎn)815萬t,每天至少需要近10萬t淡水(首鋼京唐相關(guān)報(bào)道)。由此可見,該區(qū)域地表沉降的原因主要是地下水的開采。B區(qū)位于柳贊鎮(zhèn)—一個(gè)以養(yǎng)殖海產(chǎn)品為主要經(jīng)濟(jì)產(chǎn)業(yè)的地區(qū),其附近有許多養(yǎng)殖池,如蝦池、蟹池。根據(jù)資料顯示,2019年3月柳贊鎮(zhèn)南部5萬多畝露天海水養(yǎng)殖區(qū)域已全數(shù)建成并進(jìn)入海水養(yǎng)殖準(zhǔn)備階段。由此可見,B區(qū)沉降嚴(yán)重的主要原因?yàn)轲B(yǎng)殖池的建造與開挖。C區(qū)位于唐山灣生態(tài)城,曹妃甸東站附近。生態(tài)城復(fù)合功能用地規(guī)劃用地面積387.90 ha,占城市建設(shè)用地總面積的20.28%左右[9],曹妃甸東站于2017年開始建設(shè),2019年運(yùn)營。唐山灣生態(tài)城即將被打造為臨海宜居城市,建設(shè)面積在不斷地增加,可見養(yǎng)殖池的平填與地上建設(shè)用地的建造為C區(qū)地表提升的主要原因。

    為了更加直觀的了解PS-InSAR技術(shù)所處理的沉降結(jié)果,分別在沉降尤為嚴(yán)重的A區(qū)、B區(qū)選定了兩個(gè)PS點(diǎn)——Dot1(38°56′10.88′′N、118°31′44.12′′E)、Dot2(39°13′39.32′′N、118°37′25.82′′E),其位置在圖7中標(biāo)記顯示。這兩個(gè)點(diǎn)的沉降結(jié)果如圖8所示,在時(shí)間區(qū)域內(nèi),兩者的沉降過程以非線性的方式沉降,且其沉降坡度逐漸增大,由此可見,A、B兩區(qū)域的沉降極為嚴(yán)重,也間接印證了圖7沉降速率的準(zhǔn)確性。Dot1的最大累積沉降量為21.856 mm,Dot2的最大累積沉降量為15.625 mm。從這張折線圖可以看出,在沉降過程中,雖然總體是呈現(xiàn)沉降趨勢(shì),但是某一時(shí)間區(qū)域內(nèi)會(huì)有些許抬升,這主要是與該區(qū)域的地質(zhì)條件、地下水開采及地上建設(shè)有極大關(guān)系。

    圖8 A、B區(qū)域中Dot1、Dot2的累積沉降量

    3 結(jié) 語

    本文采用PS-InSAR技術(shù)對(duì)曹妃甸填海區(qū)進(jìn)行地表形變監(jiān)測(cè),得到了很好的監(jiān)測(cè)效果。從結(jié)果看出該地區(qū)的確有著極大的地表變形,最大年平均沉降率達(dá)到了-26.558 76 mm/a,沉降區(qū)域最為嚴(yán)重的地方在鋼廠附近,年平均沉降率達(dá)到了-23 mm/a。該區(qū)域的沉降與本身填海造陸且松軟的地質(zhì)條件有一定的關(guān)聯(lián),同時(shí),工業(yè)區(qū)建設(shè)運(yùn)營、長期的地下水開采、道路的開通運(yùn)營會(huì)對(duì)這個(gè)本來就是填海造陸的地區(qū)產(chǎn)生了一定的負(fù)擔(dān)。本文也證實(shí)了采用PS-InSAR技術(shù)進(jìn)行地表沉降監(jiān)測(cè)的有效性,為曹妃甸的開發(fā)提供了參考。

    猜你喜歡
    曹妃甸海區(qū)差分
    不正規(guī)半日潮海區(qū)高(低)潮選取的探討
    中國銀行股份有限公司曹妃甸自貿(mào)區(qū)分行
    數(shù)列與差分
    中國銀行股份有限公司曹妃甸分行
    地鐵長大過海區(qū)間隧道人員疏散模擬分析
    3個(gè)不同海區(qū)鳶烏賊漁業(yè)生物學(xué)的初步比較
    崛起的曹妃甸
    靜海區(qū)林業(yè)結(jié)構(gòu)調(diào)整與經(jīng)濟(jì)效益研究
    基于差分隱私的大數(shù)據(jù)隱私保護(hù)
    相對(duì)差分單項(xiàng)測(cè)距△DOR
    太空探索(2014年1期)2014-07-10 13:41:50
    老女人水多毛片| 舔av片在线| 成人午夜精彩视频在线观看| 国产精品久久电影中文字幕| 久久精品久久久久久久性| 亚洲欧美中文字幕日韩二区| 国产人妻一区二区三区在| 尤物成人国产欧美一区二区三区| 免费观看精品视频网站| 久久久久性生活片| 在线天堂最新版资源| 两个人的视频大全免费| 三级国产精品欧美在线观看| 久久久午夜欧美精品| 国产三级中文精品| 欧美bdsm另类| 精品久久久久久久久亚洲| 国产成人91sexporn| 少妇丰满av| 久久99热这里只有精品18| 最好的美女福利视频网| 国产极品精品免费视频能看的| 久久精品91蜜桃| 欧美精品一区二区大全| 色哟哟哟哟哟哟| 日产精品乱码卡一卡2卡三| 中文欧美无线码| 一区福利在线观看| 亚洲自偷自拍三级| 九九久久精品国产亚洲av麻豆| h日本视频在线播放| 成年av动漫网址| 色尼玛亚洲综合影院| 国产精品久久电影中文字幕| 麻豆国产av国片精品| 亚洲18禁久久av| 成人国产麻豆网| 久久精品久久久久久噜噜老黄 | 国产极品精品免费视频能看的| 亚洲中文字幕一区二区三区有码在线看| 国产免费一级a男人的天堂| ponron亚洲| 日产精品乱码卡一卡2卡三| 亚洲欧洲国产日韩| 国产精品国产三级国产av玫瑰| 久久精品影院6| 美女xxoo啪啪120秒动态图| av在线播放精品| 亚洲五月天丁香| 69人妻影院| 99热这里只有是精品在线观看| 亚洲欧美日韩高清在线视频| 国产一级毛片七仙女欲春2| 99久久中文字幕三级久久日本| 午夜免费激情av| 亚洲美女搞黄在线观看| 小说图片视频综合网站| 美女黄网站色视频| 国产亚洲精品久久久com| 国产又黄又爽又无遮挡在线| 欧美日韩一区二区视频在线观看视频在线 | 亚洲精品日韩av片在线观看| 亚洲内射少妇av| 国内精品宾馆在线| 免费人成在线观看视频色| 国产极品精品免费视频能看的| 99久久中文字幕三级久久日本| 亚洲欧美清纯卡通| 黄色欧美视频在线观看| 国产三级中文精品| 欧美日韩乱码在线| 国产91av在线免费观看| 看十八女毛片水多多多| 免费不卡的大黄色大毛片视频在线观看 | avwww免费| а√天堂www在线а√下载| 女人被狂操c到高潮| 国产精品野战在线观看| 99久久精品一区二区三区| or卡值多少钱| 美女高潮的动态| 久久九九热精品免费| 深爱激情五月婷婷| 校园人妻丝袜中文字幕| 国产成人影院久久av| 99在线视频只有这里精品首页| 麻豆国产av国片精品| 少妇人妻精品综合一区二区 | 99久久精品一区二区三区| 久久久久国产网址| 人妻制服诱惑在线中文字幕| 在线a可以看的网站| 亚洲激情五月婷婷啪啪| 久久精品国产自在天天线| 亚洲最大成人av| 青春草亚洲视频在线观看| 麻豆成人午夜福利视频| 国产av在哪里看| 此物有八面人人有两片| 男的添女的下面高潮视频| 亚洲精品日韩av片在线观看| 长腿黑丝高跟| 亚洲av一区综合| 男人的好看免费观看在线视频| 在线免费观看不下载黄p国产| 晚上一个人看的免费电影| 亚洲在线自拍视频| 欧美高清成人免费视频www| 特级一级黄色大片| 久久久久国产网址| 成人毛片a级毛片在线播放| 国产亚洲5aaaaa淫片| 国产精品,欧美在线| 青青草视频在线视频观看| 久久九九热精品免费| 国产av麻豆久久久久久久| 午夜福利高清视频| 国产黄色小视频在线观看| 国产69精品久久久久777片| 精品国内亚洲2022精品成人| 偷拍熟女少妇极品色| 日本黄大片高清| 日本黄色视频三级网站网址| 亚洲精品自拍成人| www.色视频.com| 变态另类丝袜制服| 狠狠狠狠99中文字幕| 能在线免费观看的黄片| 亚洲第一电影网av| 91狼人影院| 国产老妇伦熟女老妇高清| 国产淫片久久久久久久久| 色尼玛亚洲综合影院| av在线亚洲专区| 久久久久久久久久久免费av| www.av在线官网国产| 欧美日韩乱码在线| 中文字幕熟女人妻在线| 成人亚洲欧美一区二区av| 久久久久久久久中文| 免费av观看视频| 在线国产一区二区在线| 久久久久免费精品人妻一区二区| 国产白丝娇喘喷水9色精品| 人体艺术视频欧美日本| 国产精品久久久久久亚洲av鲁大| 亚洲国产欧美在线一区| 观看美女的网站| 一级毛片我不卡| 特大巨黑吊av在线直播| 天堂av国产一区二区熟女人妻| 人人妻人人澡人人爽人人夜夜 | 18禁在线播放成人免费| 亚洲欧美清纯卡通| 成熟少妇高潮喷水视频| 日日摸夜夜添夜夜添av毛片| 12—13女人毛片做爰片一| 日韩精品有码人妻一区| 美女国产视频在线观看| 波多野结衣高清作品| 又粗又硬又长又爽又黄的视频 | 精品欧美国产一区二区三| 成人亚洲精品av一区二区| 久久亚洲国产成人精品v| 欧美一区二区精品小视频在线| 色吧在线观看| 欧美zozozo另类| 午夜激情欧美在线| 日韩大尺度精品在线看网址| 免费无遮挡裸体视频| 精品人妻视频免费看| or卡值多少钱| 精品一区二区三区视频在线| 国产久久久一区二区三区| 91狼人影院| 国产三级在线视频| 精品久久国产蜜桃| 啦啦啦啦在线视频资源| 夜夜看夜夜爽夜夜摸| 久久久久久久午夜电影| 精品久久久久久久久久久久久| 日韩av不卡免费在线播放| 国产黄片视频在线免费观看| 亚洲第一电影网av| 一个人看的www免费观看视频| 欧美日韩综合久久久久久| 好男人视频免费观看在线| 国产人妻一区二区三区在| 精品人妻一区二区三区麻豆| 日韩国内少妇激情av| h日本视频在线播放| 在线观看一区二区三区| 免费大片18禁| 一边摸一边抽搐一进一小说| 性欧美人与动物交配| 最近最新中文字幕大全电影3| 大香蕉久久网| 亚洲精品自拍成人| 99热网站在线观看| a级毛片免费高清观看在线播放| 欧美xxxx黑人xx丫x性爽| 国产精品蜜桃在线观看 | 日本在线视频免费播放| 国产成人a区在线观看| 日韩高清综合在线| 99riav亚洲国产免费| 国产人妻一区二区三区在| 亚洲精品乱码久久久v下载方式| 天堂av国产一区二区熟女人妻| 中文亚洲av片在线观看爽| 国产蜜桃级精品一区二区三区| 69av精品久久久久久| 青春草国产在线视频 | 日本黄大片高清| 色综合亚洲欧美另类图片| 国产高清有码在线观看视频| 精品日产1卡2卡| 免费看光身美女| av天堂在线播放| 日韩强制内射视频| 韩国av在线不卡| 99热精品在线国产| 午夜福利高清视频| 免费观看人在逋| 啦啦啦啦在线视频资源| 日韩一区二区视频免费看| 不卡视频在线观看欧美| 亚洲丝袜综合中文字幕| 最近手机中文字幕大全| 国产精品人妻久久久久久| 国产午夜精品久久久久久一区二区三区| 最近视频中文字幕2019在线8| 热99re8久久精品国产| 人人妻人人看人人澡| 国产三级中文精品| а√天堂www在线а√下载| 国产久久久一区二区三区| 99视频精品全部免费 在线| www日本黄色视频网| 久久韩国三级中文字幕| 欧美成人免费av一区二区三区| 日日摸夜夜添夜夜添av毛片| 国产精品av视频在线免费观看| 国产老妇伦熟女老妇高清| 国产伦在线观看视频一区| 少妇人妻一区二区三区视频| 精品久久久久久久末码| 久久精品国产鲁丝片午夜精品| 一区二区三区高清视频在线| 色吧在线观看| 能在线免费观看的黄片| 欧美变态另类bdsm刘玥| 国产精品无大码| 久久精品综合一区二区三区| 老司机影院成人| 啦啦啦韩国在线观看视频| 1000部很黄的大片| 少妇的逼好多水| 色视频www国产| 尾随美女入室| 男的添女的下面高潮视频| 久久久精品94久久精品| 最近的中文字幕免费完整| 免费一级毛片在线播放高清视频| 日韩人妻高清精品专区| 青青草视频在线视频观看| 亚洲av成人精品一区久久| 午夜福利在线在线| av视频在线观看入口| 男人和女人高潮做爰伦理| 国产免费男女视频| 亚洲精品成人久久久久久| 欧美成人一区二区免费高清观看| 亚洲欧美精品综合久久99| 成人毛片60女人毛片免费| 久久婷婷人人爽人人干人人爱| 国产亚洲av嫩草精品影院| av专区在线播放| 欧美bdsm另类| 最近手机中文字幕大全| 少妇的逼水好多| 欧美bdsm另类| 午夜福利视频1000在线观看| 国产精品人妻久久久久久| 亚洲激情五月婷婷啪啪| 亚洲最大成人av| 老司机福利观看| 亚洲一区二区三区色噜噜| 哪里可以看免费的av片| 麻豆乱淫一区二区| 国产伦理片在线播放av一区 | 在线免费观看的www视频| 毛片一级片免费看久久久久| 国产乱人偷精品视频| 国产免费一级a男人的天堂| 激情 狠狠 欧美| 国内久久婷婷六月综合欲色啪| 在线国产一区二区在线| 91av网一区二区| 插阴视频在线观看视频| 在线国产一区二区在线| avwww免费| 99热全是精品| 欧美激情久久久久久爽电影| 日韩三级伦理在线观看| 精品欧美国产一区二区三| 日韩大尺度精品在线看网址| 少妇人妻一区二区三区视频| 亚洲婷婷狠狠爱综合网| 日本在线视频免费播放| 99久久久亚洲精品蜜臀av| 欧美性感艳星| 观看免费一级毛片| 97人妻精品一区二区三区麻豆| 男人和女人高潮做爰伦理| 97超碰精品成人国产| 日韩欧美精品v在线| 男插女下体视频免费在线播放| 一本精品99久久精品77| 国产女主播在线喷水免费视频网站 | 人妻久久中文字幕网| 12—13女人毛片做爰片一| 青青草视频在线视频观看| 99热这里只有是精品50| 国产精品综合久久久久久久免费| 亚洲欧洲日产国产| 免费黄网站久久成人精品| 男女啪啪激烈高潮av片| 99riav亚洲国产免费| 伦精品一区二区三区| 亚洲国产欧美在线一区| 国产精品久久电影中文字幕| 99国产极品粉嫩在线观看| 卡戴珊不雅视频在线播放| 国产在线男女| 色综合亚洲欧美另类图片| 免费黄网站久久成人精品| 国产爱豆传媒在线观看| 欧美日韩乱码在线| 99久久人妻综合| 悠悠久久av| 99国产精品一区二区蜜桃av| 联通29元200g的流量卡| 美女脱内裤让男人舔精品视频 | 久久人人爽人人片av| 少妇高潮的动态图| 99热只有精品国产| 精品无人区乱码1区二区| 黄色配什么色好看| 高清日韩中文字幕在线| 男人和女人高潮做爰伦理| 赤兔流量卡办理| 国产三级在线视频| 欧美激情在线99| 嘟嘟电影网在线观看| 99久久九九国产精品国产免费| 精品国产三级普通话版| 麻豆成人午夜福利视频| 亚洲在线观看片| 一个人看的www免费观看视频| 九九久久精品国产亚洲av麻豆| 中国美女看黄片| 欧美最黄视频在线播放免费| 国产亚洲5aaaaa淫片| 午夜福利在线观看吧| 少妇的逼水好多| 搡女人真爽免费视频火全软件| 久久久国产成人免费| 精品久久久久久久人妻蜜臀av| 国产午夜福利久久久久久| 国产精品一区二区三区四区久久| 成年免费大片在线观看| 亚洲经典国产精华液单| av免费在线看不卡| 日本黄色片子视频| 成人特级黄色片久久久久久久| 欧美日韩乱码在线| 啦啦啦观看免费观看视频高清| 国产一区二区在线观看日韩| 免费看日本二区| 只有这里有精品99| 如何舔出高潮| 亚洲欧洲国产日韩| 国产av一区在线观看免费| 国产伦精品一区二区三区视频9| 国产91av在线免费观看| 九九久久精品国产亚洲av麻豆| 免费黄网站久久成人精品| 欧美精品一区二区大全| а√天堂www在线а√下载| 色视频www国产| 久久久久性生活片| 国产精品99久久久久久久久| 国产私拍福利视频在线观看| 国产成人a∨麻豆精品| 午夜视频国产福利| 草草在线视频免费看| 简卡轻食公司| 美女脱内裤让男人舔精品视频 | 狂野欧美激情性xxxx在线观看| 国产精品永久免费网站| 神马国产精品三级电影在线观看| av视频在线观看入口| 一本精品99久久精品77| 三级国产精品欧美在线观看| 日日干狠狠操夜夜爽| 美女被艹到高潮喷水动态| 色综合亚洲欧美另类图片| 久久精品国产亚洲网站| 天堂中文最新版在线下载 | 激情 狠狠 欧美| 老女人水多毛片| 亚洲成人中文字幕在线播放| 亚洲欧美中文字幕日韩二区| av福利片在线观看| a级毛片a级免费在线| 国内久久婷婷六月综合欲色啪| 久久久久久九九精品二区国产| 国产成人a∨麻豆精品| 亚洲,欧美,日韩| 国产精品爽爽va在线观看网站| 亚洲精华国产精华液的使用体验 | 成年版毛片免费区| 18禁裸乳无遮挡免费网站照片| 国内精品美女久久久久久| 亚洲国产色片| 免费看a级黄色片| 婷婷色av中文字幕| 菩萨蛮人人尽说江南好唐韦庄 | 国产亚洲精品av在线| 青青草视频在线视频观看| 一个人看视频在线观看www免费| 一个人观看的视频www高清免费观看| .国产精品久久| 一本一本综合久久| 男的添女的下面高潮视频| 欧美色欧美亚洲另类二区| 国产一区二区三区av在线 | 亚洲av成人av| 麻豆一二三区av精品| 我要搜黄色片| 岛国毛片在线播放| 毛片女人毛片| 亚洲精品亚洲一区二区| 干丝袜人妻中文字幕| 亚洲精品日韩av片在线观看| 级片在线观看| 九色成人免费人妻av| 国产精品人妻久久久久久| 亚洲av一区综合| h日本视频在线播放| 午夜免费激情av| 久久久久久九九精品二区国产| 秋霞在线观看毛片| 3wmmmm亚洲av在线观看| 99九九线精品视频在线观看视频| 午夜精品一区二区三区免费看| АⅤ资源中文在线天堂| 日韩成人伦理影院| 六月丁香七月| 麻豆久久精品国产亚洲av| 久久精品国产亚洲网站| 国产一区二区在线观看日韩| 成人美女网站在线观看视频| 久久久久久大精品| 欧美+日韩+精品| 久久韩国三级中文字幕| 悠悠久久av| 亚洲成a人片在线一区二区| 日韩中字成人| 国产精品麻豆人妻色哟哟久久 | 丰满人妻一区二区三区视频av| 青春草视频在线免费观看| 成人漫画全彩无遮挡| 伊人久久精品亚洲午夜| 能在线免费看毛片的网站| ponron亚洲| 精品一区二区三区视频在线| 亚洲中文字幕日韩| 国产精品女同一区二区软件| 免费看光身美女| 精品午夜福利在线看| av黄色大香蕉| 少妇人妻一区二区三区视频| 欧美xxxx性猛交bbbb| 国产亚洲精品久久久久久毛片| 91在线精品国自产拍蜜月| 国产亚洲精品av在线| 噜噜噜噜噜久久久久久91| 亚洲在线观看片| 国产精品一及| 精品久久久久久久久久久久久| 亚洲美女视频黄频| 久久精品国产亚洲av天美| av免费观看日本| 亚洲欧美精品综合久久99| 免费看a级黄色片| 精品人妻熟女av久视频| 亚洲成人av在线免费| 99国产精品一区二区蜜桃av| 亚洲最大成人中文| 狠狠狠狠99中文字幕| 两性午夜刺激爽爽歪歪视频在线观看| 免费观看人在逋| 91狼人影院| 国产一区二区三区在线臀色熟女| 亚洲成av人片在线播放无| av在线播放精品| 人体艺术视频欧美日本| 少妇人妻精品综合一区二区 | 高清毛片免费观看视频网站| 一进一出抽搐gif免费好疼| 国产v大片淫在线免费观看| 日日撸夜夜添| 在线国产一区二区在线| 国产精品久久视频播放| 12—13女人毛片做爰片一| 亚洲欧美精品综合久久99| 精品午夜福利在线看| 91狼人影院| 亚洲精品国产av成人精品| 日本免费一区二区三区高清不卡| 亚洲av不卡在线观看| 国产一区亚洲一区在线观看| 国产成人一区二区在线| 久久九九热精品免费| 欧美日本亚洲视频在线播放| 老女人水多毛片| 激情 狠狠 欧美| 在线观看免费视频日本深夜| 国产成人aa在线观看| 国产 一区精品| 丝袜喷水一区| 日韩强制内射视频| 九九在线视频观看精品| 日本撒尿小便嘘嘘汇集6| 高清日韩中文字幕在线| 少妇的逼水好多| 好男人视频免费观看在线| 欧美精品一区二区大全| 九九爱精品视频在线观看| 精品久久久久久久久av| 国产精品国产高清国产av| 麻豆国产av国片精品| 欧美日韩精品成人综合77777| 亚洲在线观看片| 18禁在线播放成人免费| 天堂中文最新版在线下载 | 亚洲精品日韩在线中文字幕 | 成年女人永久免费观看视频| 伊人久久精品亚洲午夜| 成人av在线播放网站| 日本色播在线视频| 真实男女啪啪啪动态图| 嫩草影院新地址| 欧美日韩一区二区视频在线观看视频在线 | 国产激情偷乱视频一区二区| 卡戴珊不雅视频在线播放| 久久久久久九九精品二区国产| 内地一区二区视频在线| 国产成人精品婷婷| 国产白丝娇喘喷水9色精品| 可以在线观看的亚洲视频| 日韩中字成人| 小蜜桃在线观看免费完整版高清| 日韩av在线大香蕉| 久久欧美精品欧美久久欧美| 天堂影院成人在线观看| 一个人免费在线观看电影| 成人特级黄色片久久久久久久| 婷婷色av中文字幕| 又粗又硬又长又爽又黄的视频 | 综合色丁香网| 色哟哟哟哟哟哟| 成熟少妇高潮喷水视频| a级毛片免费高清观看在线播放| 人人妻人人澡人人爽人人夜夜 | 日本免费a在线| 在线天堂最新版资源| 国产精品久久久久久久久免| 国产一区亚洲一区在线观看| 99久久九九国产精品国产免费| 99热6这里只有精品| 亚洲av男天堂| 天天一区二区日本电影三级| 国产精品一区www在线观看| 黄片wwwwww| 波多野结衣巨乳人妻| 午夜免费激情av| 精品久久久久久久久亚洲| 中文字幕久久专区| 亚洲高清免费不卡视频| 亚洲精品日韩在线中文字幕 | 国产精品女同一区二区软件| 寂寞人妻少妇视频99o| 亚洲av电影不卡..在线观看| 噜噜噜噜噜久久久久久91| 欧美日韩一区二区视频在线观看视频在线 | 亚洲欧美日韩东京热| 欧美色视频一区免费| 久久精品夜夜夜夜夜久久蜜豆| 中国美白少妇内射xxxbb| 免费观看的影片在线观看| 精品久久久久久久久亚洲| 国产成人午夜福利电影在线观看| 一本久久中文字幕| 国产一区二区在线观看日韩| 亚洲人成网站在线观看播放| 亚洲最大成人手机在线| 亚洲欧美成人综合另类久久久 | 男人舔奶头视频| 久久久a久久爽久久v久久| 两个人视频免费观看高清| 日韩人妻高清精品专区| 天堂影院成人在线观看|